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Este livro é direcionado, principalmente, para os cursos de tecnologia que tenham em seu 
campo de abrangência a noção e a preparação de métodos científicos para a análise de situa- 
ções complexas (inseridas no contexto empresarial) e a tomada de decisões sobre elas. 


Você verá, no transcorrer deste estudo, bem como de suas atividades profissionais, que a pes- 
quisa operacional (PO) é cada vez mais útil no âmbito dos processos e dos métodos destinados 
à customização e à otimização em contextos reais, como produção, armazenagem, distribuição 
e planejamento. 


A PO é amplamente utilizada pelas diversas áreas do conhecimento, como agricultura, finan- 
ças, marketing, meio ambiente, serviços públicos e, em especial, engenharia da produção, que 
emprega muito esse recurso para a compreensão de características de sistemas complexos de 
modelos matemáticos para desenvolvimento de simulações, as quais consistem em diversos 
métodos de resolução, na busca da solução ótima. 


Neste livro, propomo-nos a, de maneira simples, iniciar você, tecnólogo, na compreensão dos 
resultados mostrados pela PO por meio de métodos específicos, para instruí-lo ou subsidiá- 
-lo para a tomada de decisões. 


Deixamos para você, como tarefa de pesquisa, a abordagem histórica da criação da PO, suas 


diversas vertentes e os cenários de sua origem. 


No primeiro capítulo, mostraremos a influência da PO nos processos de tomada de decisão 
nas organizações por meio do uso dos modelos matemáticos, para que você possa ter uma 
ideia da importância dessa ferramenta nos âmbitos estratégico, tático e operacional. No se- 
gundo capítulo, apresentaremos o contexto algébrico da programação linear (PL), com as ca- 
racterizações de funções e de restrições na busca pela otimização de resultados, bem como 
examinaremos as restrições de problemas de PL por meio da interpretação geométrica com a 
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respectiva demonstração. Daremos continuidade a essa tratativa no terceiro capítulo, com a 
exposição do método simplex e dos seus rearranjos para as diversas situações, destacando a 
necessidade de observar nesse processo fatores específicos de cada circunstância e aspectos 


próprios que sejam inerentes à sua aplicação. 


Na sequência, no quarto capítulo, enfocaremos a resolução de problemas de PL por meio do 
software WinQSB. Finalizando, no quinto capítulo, mostraremos o uso da PO na resolução 
dos problemas de transportes e de redes, e a respectiva resolução de problemas por meio do 
WinQSB será abordada no sexto capítulo. As concepções sobre análise de sensibilidade e so- 
bre os processos de simulação serão tratados no sétimo capítulo. Como adendo, apresenta- 
remos no oitavo capítulo uma abordagem sobre o conhecimento matemático matricial (por 
se tratar de prerrequisito à abordagem de PO). Isso será feito a título de rememoração, pois 
esse conteúdo serve de fundamentação, de base, para a resolução algébrica de problemas de 


PO no contexto da modelagem matemática. 


Esta obra, além de conceitos essenciais, traz diversos exercícios resolvidos e para resolver. Essa 

foi a metodologia pela qual optamos (com exemplos e contextualizações), com o objetivo de 

apresentar o passo a passo inicial da PO para você que não é da área das ciências matemáti- 
cas. Além disso, com essa abordagem, a intenção é facilitar a visualização dos procedimentos 

operacionais e a compreensão da importância do uso de tal ferramenta. Destacamos tratar-se 

aqui de uma abordagem introdutória. Você poderá, com base neste estudo, pesquisar e am- 
pliar o seu entendimento, principalmente com conhecimentos advindos do manuseio de pro- 
gramas como o Microsoft” Office Excel, o Linear Interactive and Discret Optimizer (Lindo), 
o Language for Interactive General Optimizer (Lingo) e o Quantitative System for Business 

(WinQSB). 
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E ) Como aproveitar ao máximo este livro 


Conteúdos do capítulo 


* Contexto de aplicação da pesquisa operacional. 
e Passos para estruturar uma pesquisa operacional. 


e Ferramentas matemáticas de modelagem de uma pesquisa operacional. 


Após o estudo deste capítulo, você será capaz de: 


1. reconhecer uma situação-problema; 


2. aplicar os passos de elaboração de uma pesquisa operacional; 


3. tomar decisões com o uso dessa ferramenta. 


Logo na abertura do capítulo, você fica conhecendo os conteúdos que serão nele abordados. 
Você também é informado a respeito das competências que irá desenvolver e dos conheci- 
mentos que irá adquirir com o estudo do capítulo. 


TU “Exercício resolvido 


Vamos formular um exemplo de fácil percepção para a formulação do mo- 
delo matemático chamado de forma-padrão. 


Vamos considerar na formulação matemática do modelo-padrão uma em- 
presa de artefatos de madeira que fabrica: 


* seis porta-trecos por hora, se fizer somente os porta-trecos; 


e quatro potes de madeira para armazenar mantimentos, se fizer so- 
mente potes. 


' Ocorre que nessa produção a empresa gasta: 


* duas unidades de madeira para fabricar uma unidade de 
porta-trecos; 


Ao longo do livro, há uma série de exercícios resolvidos para facilitar a compreensáo dos capítulos. 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal. 


Perguntas & respostas 
O que é método? 


Iremos recorrer ao Dicionário Houaiss da língua portuguesa (Houaiss; 
Villar, 2009), o qual apresenta várias definições com diferentes abran- 
gências e estabelece áreas de referência para o significado do termo. 
Encontramos, por exemplo, que método é “procedimento, técnica 
ou meio de fazer alguma coisa, especialmente de acordo com um pla- 
no”. Na área filosófica, quando falamos em método, referimo-nos a 
um “conjunto sistemático de regras e procedimentos que, se respeita- 
dos em uma investigação cognitiva, conduzem-na à verdade”. E você, 


o que entende sobre método? 


Nesta seção, os autores respondem a dúvidas frequentes relacionadas aos conteúdos do capítulo. 


7l "Síntese 


O uso de ferramentas quantitativas que ofereçam subsídios objetivos 
para a atividade de gerenciamento foi o enfoque que privilegiamos neste 
capítulo. Nesse contexto, a PO apresenta-se como uma ciência prática. Ela 
estabelece parâmetros decisórios confiáveis, nos quais considera fatores 
e cenários dos problemas e estabelece, através de modelos matemáticos, 
a possibilidade de visualizarmos as possíveis soluções para situações 
cujas variáveis, restrições e funções objetivo são lidas a partir de cálculos. 
Cálculos estes modelados a partir das fases de estruturação do problema. 
Por isso também é necessário que você tenha clareza sobre o que seja 
método, sistema e solução ótima (conceitos que desenvolvemos neste 
estudo). Assim, poderá aproveitar esse instrumental nos processos de 


tomada de decisão. 


Você dispõe, ao final do capítulo, de uma síntese que traz os principais conceitos nele abordados. 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal 
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Al Questões para revisão 


1. Em relação às fases da pesquisa operacional, assinale (V) para verdadei- 
ro e (F) para falso: 


) Construção e alteração do modelo. 

) Estabelecimento e controle das soluções. 
) Cálculo do modelo. 

) Consulta aos operários. 

) Teste do modelo e da solução. 

) Compra de software. 


) Formulação do problema. 


Com estas atividades, vocé tem a possibilidade de rever os principais conceitos analisados. 


Ao final do livro, os autores disponibilizam as respostas às questóes, a fim de que vocé 
possa verificar como está sua aprendizagem. 


Para saber mais 


GOMIDE, F. Algoritmos numéricos de busca em otimização. 
Conteúdo para a aula de Planejamento e Análise de Sistemas de 
Produção da Faculdade de Engenharia Elétrica e de Computação da 
Unicamp. Disponível em: <http://www.dca.fee.unicamp.br/-gomide/ 
courses/EA044/transp/EA. 044 BuscaNumericaOtimizacao.pdf>. 
Acesso em: 20 nov. 2009. 

Nesse material, vocé encontra um excelente resumo sobre algoritmos 
de busca, aplicações, problemas de maximização e minimização, bem 


como sobre o método BIG-M. 


Vocé pode consultar as obras indicadas nesta secáo para aprofundar sua aprendizagem. 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal 


capitulo 


Os subsídios da pesquisa operacional para 
o processo de tomada de decisões 
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Conteúdos do capítulo 


e Contexto de aplicação da pesquisa operacional. 
e Passos para estruturar uma pesquisa operacional. 


e Ferramentas matemáticas de modelagem de uma pesquisa operacional. 


Após o estudo deste capítulo, você será capaz de: 


1. reconhecer uma situação-problema; 
2. aplicar os passos de elaboração de uma pesquisa operacional; 


3. tomar decisões com o uso dessa ferramenta. 


Ao nos acompanhar neste estudo, você irá verificar que a programação linear (PL) é uma 
das mais importantes técnicas de otimização da pesquisa operacional (PO), contexto no qual 
muitos problemas podem ser descritos como problemas de otimização sujeitos a determina- 


das restrições. 


Portanto, para que você tenha uma melhor compreensão do ambiente de uso da PL, é impor- 
tante que tenha claro o que é a PO. Muitos autores a definem, de modo geral, como um mé- 
todo para a tomada de decisões. 


Para Loesch e Hein (2009, p. 1), a PO se classifica como uma ciência do conhecimento. Eles 
afirmam que ela “como ciência estrutura processos, propondo um conjunto de alternativas e 
ações, fazendo a previsão e a comparação de valores, de eficiência e custos.” 


Já de acordo com a concepção de Silva et al. (2008, p. 11), a PO se constitui em um sistema: 
“consiste na descrição de um sistema organizado com o auxílio de um modelo, e através da 
experimentação com o modelo, na descoberta da melhor maneira de operar o sistema”. 


Preferimos dizer que a PO é uma ferramenta matemática que auxilia no processo de tomadas 


de decisão em situações reais. Para isso, utilizamos modelos matemáticos estruturados em fases. 
Perguntas & respostas 
O que é método? 


Iremos recorrer ao Dicionário Houaiss da língua portuguesa (Houaiss; Villar, 2009), o qual apre- 
senta várias definições com diferentes abrangências e estabelece áreas de referência para 
o significado do termo. Encontramos, por exemplo, que método é “procedimento, técnica 
oumeio de fazer alguma coisa, especialmente de acordo com um plano”. Na área filosófica, 
quando falamos em método, referimo-nos a um “conjunto sistemático de regras e procedi- 
mentos que, serespeitadosem uma investigação cognitiva, conduzem-naà verdade”. Evocê, 


o que entende sobre método? 


Estamos falando em modelos e métodos. O que isso significa? 
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Podemos dizer que a PO é uma ciéncia aplicada. Ela nos oferece, de acordo com Lisboa (2002), 
instrumentos para: 

e resolver problemas reais; 

e tomar decisões com base em dados e correlações quantitativos; 


e conceber, planejar e operar sistemas fazendo uso de tecnologia e métodos de outras 


esferas do conhecimento; 
e diminuir os custos e aumentar o lucro; 


e encontrar a solução ótima. 
Perguntas & respostas 
O que é uma solução ótima? 
Podemos dizer que é aquela que melhor serve aos objetivos das pessoas e das organizações. 


Como exemplo, podemos citar o propósito de encontrar o lucro máximo ou o custo mínimo. 


A PO é bastante empregada como ferramenta nos processos de tomada de decisões no en- 
frentamento de problemas dos ambientes de negócios (é este o seu aspecto que destacamos 
neste estudo). Por sua vez, os principais instrumentos que a PO utiliza são oriundos dos co- 
nhecimentos matemáticos, estatísticos e de informática. 


Você irá constatar, no transcorrer deste nosso estudo, que se trata da utilização da modela- 
gem matemática aplicada à área de realização de negócios e à área empresarial, consideran- 
do-se tanto o setor privado como o setor püblico. Nesse sentido, Lachtermacher (2009, p. 2) 


relaciona uma série de situações na qual a PO pode ajudar no processo de decisão. São elas: 
e problemas de otimização de recursos; 
e problemas de localização; 


e problemas de roteirização; 
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e problemas de carteiras de investimentos; 


e problemas de alocação de pessoas; 


e problemas de previsão de planejamento; 


e problemas de alocação de verbas de mídia. 
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Nessa atividade!, é necessário realizarmos procedimentos no sentido de organizarmos os 
dados coletados em sistemas computacionais (SIG) para usuários de outras áreas, portanto, 
não técnicos, a fim de que sejam úteis no apoio ao processo gerencial de tomada de decisões. 


Perguntas & respostas 
O que significa SIG? 


No nosso caso, quando falamos em SIG, estamos nos referindo aos sistemas de infor- 
mações gerenciais. Essa sigla é também um acrônimo de SIG - Sistemas de Informação 
Geográfica; não devemos confundir uma com a outra. Os sistemas de informações ge- 
renciais interligam e consolidam informações operacionais, bem como dados anteriores 
dos procedimentos similares, fornecendo, dessa forma, um suporte quantitativo e qua- 
litativo para as tomadas de decisões. Para resumir, podemos dizer, conforme definição 
bastante objetiva encontrada no site do Cepromat (2009), que se constitui no “processo de 
transformação de dados em informações. E, quando esse processo está voltado para a ge- 
ração de informações que são necessárias e utilizadas no processo decisório da empresa, 


diz-se que esse é um sistema de informações gerenciais.” 


Os procedimentos de PO em uma organização abrangem, interligam e beneficiam todos os 
níveis, ou seja, auxiliam as decisões nos âmbitos estratégico, tático e operacional. Esses três 
níveis de uma organização correspondem respectivamente, conforme Sertek, Guindani e 
Martins (2009), à atuação daqueles que definem o rumo da empresa (diretores, assessores), às 
ações de quem atua na liderança e no comando da execução das atividades (coordenadores, 


supervisores) e, por fim, à atuação dos que executam propriamente as tarefas. 


s por que precisamos tomar deci 


Obviamente, essa é uma condição da vida humana. Constantemente somos chamados a to- 
madas de decisões, desde a escolha do supermercado até a compra de um carro ou a cor de 


uma roupa, entre tantas outras. 


No entanto, nosso estudo é direcionado para questões do ambiente organizacional. Esse é 
o marco que delimita nossa abordagem. Mas, como em qualquer situação da vida, esse tema 


apresenta variáveis, as quais caracterizam as situações-problema. 


1 O enfoque de PO que apresentamos nesta obra está inserido na Management Sciences (MS). Conforme explica Lachtermacher 
(2009, p. 2), esta é a “área de estudos que utiliza computadores, estatística e matemática para resolver problemas de negócios”. 
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Nenhuma parte desta 


Afinal, quais sáo os problemas? Vocé verá, na continuidade, por meio dos exemplos, que 
os problemas são situações que a organização precisa resolver para dar seguimento aos seus 
propósitos ou para atingir seus objetivos. 


Diante desses problemas, encontramos fatores que interferem nos processos de tomada de 
decisão e que, quando estamos na condição de gerenciar a situação, precisam ser considera- 
dos. Veja no Quadro 1.1. 


Quadro 1.1 — Fatores e cenários que interferem na tomada de decisão 


Fatores Cenários 
"as Está relacionada ao impacto que a decisão pode provocar na organização 
Importância 2,4 
(ganhos ou prejuízos). 
O nümero de decisores, individual ou em grupo, simplifica ou torna mais 
Agentes 
complexo o processo. 
Risco As certezas ou as incertezas influenciam nossas decisóes. 
. São aspectos sociais e culturais que interferem no processo 
Ambiente Pr 
decisório. 
i Surgem em função de choques de interesses entre setores de uma organiza- 
Conflitos E : 
ção ou entre decisores. 


Fonte: Elaborado com base em Lachtermacher, 2009, p. 4. 


É nesse contexto que os modelos utilizados em PO representam uma ferramenta altamente 
qualificável para o nosso trabalho de gestão, em todos os âmbitos — gerencial, estratégico ou 
operacional —, pois, além de fornecerem subsídios para uma decisão individual, facilitam a 
comunicação entre os decisores, entre os vários setores de uma organização (quando da de- 
cisão em grupo) e permitem quantificar as variáveis envolvidas na decisão. Portanto, tornam 
o processo de tomada de decisão objetivo”. 


A modelagem apresenta processos que são usados em várias situações. Entretanto, vamos 
nos focar no estudo dos problemas que requerem uma solução, uma decisão. Nessa cir- 
cunstância, você pode fazer a modelagem da situação para realizar várias simulações dos 
possíveis panoramas que irão configurar-se diante de determinadas alternativas e, assim, visu- 
alizar as consequências de uma tomada de decisão. Portanto, podemos dizer que a modelagem 
oferece diretrizes para as decisões nos vários âmbitos da organização, ou seja, a modelagem 
oferece informações para chegarmos à melhor decisão ou à melhor solução para o problema. 


2 O termo objetivo é usado aqui no sentido que lhe dá o Dicionário Houaiss da língua portuguesa em sua terceira notação, isto é: 
“diz-se do que está no campo da experiência sensível independentemente do pensamento individual e perceptível por todos os obser- 
vadores” (Houaiss; Villar, 2009). 
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Você deve estar se perguntando se nào foi sempre assim. Não, a elaboração de modelos tor- 
nou-se viável a partir das inovações tecnológicas, ou seja, da disponibilidade de bancos de 
dados, microcomputadores, entre outros instrumentos atualmente disponíveis. Surgiu, assim, 
uma nova condição para o gerenciamento de problemas ou para a gestão das organizações. 
Conforme Lachtermacher (2009), deixamos de lado a intuição (instrumento anterior dos res- 
ponsáveis pelas decisões) e passamos a utilizar somente os modelos de situações, criados a 
partir da seleção dos dados relevantes e da simulação das possíveis soluções. 


Todavia, na opinião do referido autor, com a qual concordamos, a postura de abandono da 
intuição não é a melhor escolha. O recomendável é unir as duas fontes de informações e uti- 
lizar a intuição como ferramenta auxiliar nos processos de: 


e classificação da relevância da informação, a fim de utilizá-la ou não; 
e escolha dos cenários a serem estudados; 

e validação do modelo sugerido para a situação; 

e análise dos resultados dos estudos. 


É possível percebermos nessa descrição operacional (reveladora da interferência do aspecto 
subjetivo) que, se, de fato, atribuirmos apenas aos dados da modelagem a responsabilidade 
pela decisão, estaremos negando a capacidade humana de análise e de ponderação sobre os 
fatos. Veja essa dinâmica de interação entre os modelos utilizados em PO e a intuição do de- 
cisor na Figura 1.1. 


Figura 1.1 — Processo de tomada de decisão 


Mundo real 


Mundo simbólico 


Fonte: Lachtermacher, 2009, p. 3. 


Você pode acompanhar na figura, pela direção das setas, a concepção que reserva para a in- 
tuição a função de selecionar aspectos do modelo que forem considerados pertinentes, bem 
como a de validar perspectivas de resultados. Observe que, enquanto todo o processo ocor- 
re no espaço do chamado “mundo real” (inclusive o “mundo simbólico”), a intuição situa-se 
fora de tal espaço. 
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Uso de modelos em PO 


No que se refere aos modelos, você pode encontrar na literatura e na prática de gestão trés ti- 
pos: os físicos?, os análogos e os matemáticos ou simbólicos. No entanto, neste estudo, quan- 
do nos referimos a modelos, estamos tratando dos modelos matemáticos, que sáo os mais 


utilizados nas atividades gerenciais. 


Considerando a necessidade de quantificar informações, Lachtermacher (2009) descreve as 
características fundamentais dos modelos matemáticos: 


e as grandezas são representadas por variáveis de decisão; 
e as relações entre as variáveis são apresentadas por expressões matemáticas. 


O que ocorre é o seguinte: no processo de enunciação do problema e resolução para encon- 
trar a solução ótima, aplicamos a modelagem, ou seja, procuramos definir os aspectos que 
nos permitem delinear (no caso de programação linear, inteira, mista ou não linear) as se- 
guintes condições: 


e Variáveis do problema - São os fatores que você quer controlar e aqueles que você de- 
seja saber quanto valem. São esses os elementos que irão representar as variáveis de 


decisão (incógnitas encontradas na solução do modelo). 


e Parâmetros do problema - Correspondem aos valores fixos do problema. Normalmente, 
são valores financeiros dos dados ou das variáveis que se apresentam no problema, co- 


mo o custo fixo para a produção de um determinado produto. 


e Restrições — Você as encontra em razão das próprias variáveis do problema. São aque- 
les aspectos que limitam as combinações e consequentemente as variáveis de soluções 


possíveis. 


e Função” objetivo — E uma função matemática que você usa para maximizar ou mini- 
mizar, dependendo do seu objetivo em relação ao problema focado. Ela define a qua- 


lidade da solução, sempre em razão das incógnitas encontradas. 
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Vamos observar esse processo em um exemplo prático. Considere o exercício resolvido a seguir. 


Nenhuma parte desta 


Digamos que uma empresa — vamos chamá-la de Aves & Flora — produza ração para aves. 
Ela fabrica dois tipos: a ração Super (para a fase inicial) e a ração Extra (para a fase de cres- 
cimento). Em ambas, mistura 10% de “núcleo de frango de corte”, conforme recomendações 


3 Modelos físicos são, como exemplifica Lachtermacher (2009, p. 5), “os de aeronaves e casas utilizados por engenheiros”, enquanto 
os análogos são aqueles que representam “as relações utilizando diferentes meios”. No caso deste, podemos citar como exemplo “os 
mapas rodoviários, que representam as rodovias de uma região por meio de traços sobre um papel”. 


4 Função é uma relação entre dois conjuntos na qual existe uma ligação unívoca entre os elementos. No ensino de matemática, 
define-se função como “relação entre dois conjuntos que abrange todos os elementos do primeiro e associa a cada elemento deste 
primeiro conjunto somente um elemento do segundo” (Houaiss; Villar, 2009). 
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dos setores especializados. A diferença de custos está nos 90% de matéria-prima restante pa- 
ra completar as embalagens de 10 kg. Isso significa que: 


e na Super são utilizados 60% de milho moído (6 kg) e 30% (3 kg) de farelo de soja; 


e na Extra são utilizados 70% (7 kg) de milho moído e 20% (2 kg) de farelo de soja; 
e aembalagem de 10 kg da ração Super custa R$ 12,00 e a da ração Extra custa R$ 18,00; 
e o farelo de soja custa R$ 3,00 por quilograma, enquanto o milho moído custa R$ 2,00; 


* aempresa conta com um suprimento mensal (estoque) de 100.000 kg de milho moído 
e 20.000 kg de farelo de soja. 


A área administrativa dessa empresa quer saber qual a quantidade de ração que a empresa 


deve produzir para maximizar o lucro. O que você faria? 


Você pode aplicar aqui os passos que vimos anteriormente, ou seja, analisar esse problema 
nos seguintes aspectos: 


e Verificar as variáveis de decisão — No caso, trata-se das quantidades de ração de cada 
tipo (Super e Extra) a serem fabricadas. 


e Identificar os parâmetros — Nesse problema, são os preços de compra da matéria-pri- 
ma e de venda das rações (produto final), bem como as quantidades específicas que 
são utilizadas dessas matérias-primas (milho moído e farelo de soja) em cada um dos 
tipos de ração. 


e Levantar as restrições — São representadas nessa situação pelas quantidades de esto- 
que, ou seja, há um limite para a produção estabelecido pelo volume de milho moído 
e de farelo de soja disponíveis. 


e Definir a função objetivo — Corresponde a uma função matemática que determina 
qual volume das rações Super e Extra a empresa deve produzir para maximizar o lucro. 


Esses são aspectos de estudo para uma tomada de decisão que podemos aplicar em diferen- 
tes contextos operacionais ou sistemas, o que corresponde a uma diversidade muito grande 
de práticas, que incluem atividades econômicas, sociais e ambientais. Nessas circunstâncias, 


podemos trabalhar realidades em que ocorrem: 


e função objetivo e variáveis de restrições lineares e não lineares; 
e variáveis de decisão contínuas ou discretas; 
e parâmetros que representam uma determinação ou uma probabilidade. 


Mas, ao falarmos em modelos e modelagem, devemos destacar que as funções linear, mista 
e não linear, entre outras, são técnicas matemáticas de uso na PO. Nesse contexto, a estrutu- 
ração do modelo tem relação direta com o sistema representado, que é real. 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal 


Perguntas & respostas 


O que é um sistema? 


Consideramos que as definições que melhor se adéquam ao contexto matemático da PO, 
entre as encontradas atualmente, relacionadas com diversas áreas, sejam estas duas: “1. con- 
junto de elementos, concretos ou abstratos, intelectualmente organizados [...] 1.3 distribui- 
ção e classificação de um conjunto de elementos segundo uma ordem estabelecida Exs.: s. 


taxonómico, s. decimal" (Houaiss; Villar, 2009, grifo do original). 


Prevalece no estudo de PO, como vocé irá perceber, ao analisar os exemplos, a ideia de uma 
estrutura em que é observado um conjunto de elementos que se interligam, estabelecendo 
uma situação específica, um contexto, ou seja, a realidade concreta de fatos e fatores. Nesse 
sentido, a modelagem é um processo que utiliza ferramentas matemáticas para organizar os 
dados de uma circunstância de tal forma que sejam criados modelos facilitadores para uma 
tomada de decisão. 


Figura 1.2 — Exemplo de modelos em PO 


Cadeia, projeto e rede de valor - Comparativo entre os modelos 
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CADEIA DE Infraestrutura 
VALOR PROCESSOS Gestão dos recursos humanos 
DE SUPORTE Desenvolvimento da tecnologia 
Compras 
PROCESSOS Logística Gigas Logística Marketing e Soniye 
DE NEGÓCIO de entrada de saída vendas | pós-venda 
PROJETO REDE DE 
DE VALOR | VALOR 
Levantamento Soluções Gestão do contrato e da rede 
4 eanálisedo — parao —a PROCESSOS z d 
í 2 rovisionamento do servico 
problema E polane Escolha da DE NEGÓCIO En NIS 
i solução do Operação da infraestrutura 
É : = == : problema 
^ Controlee | a 
avaliação Implantação SEA 
Infraestrutura Infraestrutura da empresa 
Gestão dos recursos humanos Gestão dos recursos humanos 
PROCESSOS 
Desenvolvimento da tecnologia DE SUPORTE Desenvolvimento da tecnologia 
Compras Compras 


Fonte: Adaptado de Santos, 2006. 


Na Figura 1.2, podemos ver, por meio de um diagrama, os processos de uma situa- 
ção real, evidenciada como um modelo, na qual você pode observar que, dentro da 
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cadeia de valor, instala-se o projeto de valor, possibilitando a organização da rede de va- 
lor. Interligando as três etapas do modelo, os processos de suporte se mantêm inaltera- 


dos; no entanto, os processos de negócio apresentam possibilidades diferenciadas para 


a busca da solução. 


Para facilitar os procedimentos, ou seja, a aplicação dessa ferramenta, costumamos dizer 


que a PO é dividida em seis fases, na seguinte ordem: 


6. 


formulação do problema; 

construção ou alteração do modelo; 
cálculo do modelo; 

teste do modelo e da solução; 
estabelecimento e controle das soluções; 


implantação e acompanhamento. 


Se você observar cada uma dessas fases, vai perceber que elas, de fato, são as etapas a serem 


transpostas, enfrentadas e executadas na solução do problema apresentado para encontrar 


a solução ótima. Vamos, portanto, fazer um resumo descritivo dos seis passos listados. Você 


irá perceber que o quinto e o sexto passo apresentam alguns aspectos interdependentes. 


1. 


Formulação do problema — Nessa fase, determinamos o que pretendemos fazer. É o 
momento de definições, no qual devemos desenvolver o problema cuidadosamente, 
com clareza e coerência, delimitando os objetivos a alcançar e identificando as limi- 
tações ou restrições do produto em estudo, inclusive esboçando e acolhendo críticas 
aos possíveis caminhos que desejamos alcançar. É o período adequado para fazermos 
as proposições possíveis, verificando os registros e coletando novas informações, com 
máxima precisão e veracidade das informações. 


Perguntas & respostas 
E se a formulação for feita de modo equivocado? 


Você deve lembrar que, se a formulação for feita de modo equivocado, fugindo do objeto 


de estudo, todas as outras fases ficarão comprometidas e poderão levar a erro. 


Construção ou alteração do modelo — Partimos para o enfrentamento dessa fase com 
a concepção de que o modelo é uma representação da realidade estudada; portanto, 
ele depende do estudo do problema já levantado na fase inicial. É importante lembrar 
que os elementos constitutivos do modelo são oriundos de dados da empresa ou do 
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mercado, ou ainda, do cenário que está sendo analisado. Aqui prevalece a modelagem 
matemática, definida por equações e inequações matemáticas, seja na função objetivo, 
seja em suas restrições, como veremos mais adiante. Cabe lembrar que nessa fase é 
preciso separar as variáveis decisivas, também chamadas de variáveis controláveis, das 
variáveis não controláveis. Para exemplificar, podemos dizer que uma variável contro- 
lável é aquela que, em uma situação de produção, representa a quantidade produzida. 
Já a variável não controlável poderia ser a demanda dessa quantidade ou o preço de 
mercado praticado por certo produto. 


3. Cálculo do modelo - Também podemos chamar essa etapa de resolução do modelo. É 
a fase na qual encontramos a solução ou as soluções do modelo por meio de diver- 
sas técnicas de resolução, desde as mais simples, nas quais os problemas também são 
simples, até as técnicas de programação linear mais sofisticadas. O método simplex 
é uma das técnicas que podemos utilizar quando se trata de problemas com mais de 
duas variáveis controladas. Existem também muitos recursos computacionais que 
permitem fazer o cálculo com extremo rigor, confiabilidade e rapidez. Podemos des- 
tacar alguns, como Solver (Microsoftº Office Excel), Lindo, Lingo, OMP, Maple, PLM 
e WinOSB. 


4. Teste do modelo e da solução — Nessa fase, verificamos se os resultados obtidos pelas 
soluções encontradas do modelo matemático servem para o modelo real do problema, 
após sua implantação, ou ainda, se ele permite verificar se serão necessárias novas so- 
luções para melhorar ainda mais. 


5. Estabelecimento e controle das soluções — Por que controle? Na construção do mode- 
lo, bem como na sua execução, devemos identificar alguns parâmetros ou valores fixos 
para a solução do problema. Esses parâmetros devem ser controlados, para que, caso 
sofram desvios durante o processo, o modelo possa ser corrigido. 


6. Implantação e acompanhamento — E a fase em que verificamos com precisão o que con- 
seguimos, ou seja, o que encontramos nas fases anteriores, desde que tenhamos feito 
todo o acompanhamento do processo. Podemos também, nessa etapa, fazer correções 
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ou ajustes no modelo, caso sejam necessários. 


E interessante que você saiba que a pesquisa operacional atua em vários níveis do conheci- 
mento, que podem ser formulados matematicamente, adotando-se uma técnica específica 


para cada situação, como nos seguintes casos: 


e em setores cujo objetivo é otimizar a quantidade produzida para alcançar o menor cus- 
to ou a maior receita (aumentar ou diminuir), entre os quais podemos citar a agricul- 


tura, a indústria química ou a produção industrial; 


e nas indústrias moveleira e metalúrgica, que procuram minimizar o desperdício com 
problemas de cortes de chapas ou madeiras, tendo, assim, um melhor aproveitamento; 
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e em carteiras de investimentos que possam trazer uma opção melhor para um determi- 


nado investimento, em função de se obter uma melhor rentabilidade; 


e em situações de transporte para otimizar o tempo e o custo no que se refere tanto a flu- 
xo de transporte como a fluxo para obtenção do caminho mínimo. 


Para cada situação, como já dissemos, utilizamos uma ferramenta (procedimento) ou técnica, 
conforme as particularidades identificadas. Entre elas, podemos destacar as técnicas descri- 
tas a seguir. 


e Programação linear (PL) — Para Passos (2008), consiste na programação matemática que 
procura otimizar situações-problema sujeitas a certas restrições — maximizando ou mi- 
nimizando. Cabe ressaltar que existem outros tipos de PL, com características bem defi- 
nidas, como as que se classificam como inteiras, mistas ou dinâmicas. 


e Teoria das filas - Tem como finalidade trabalhar com situações de congestionamento de 


sistemas, tempo de espera em filas, entre outras. 


* Teoria dos grafos — E direcionada para os fluxos máximos, o caminho mais curto, a ro- 
teirização de caminhos para veículos, o planejamento e a programação de projetos — 
Pert e CPM. 


e Simulação — Atua com modelos representativos, analisando o comportamento de vari- 
áveis e abrangendo situações que favorecem a lei do acaso. Como exemplo, podemos 
citar a simulação de Monte Carlo. 


* Teoria dos jogos — Baseia-se em estratégias para persuasão e tomada de decisão. 


Você pode usar todas essas técnicas para o estudo do problema e a quantificação de dados 
como subsídios para tomar decisões na empresa. No entanto, em nosso estudo, iremos privi- 


legiar a PL, a qual veremos no próximo capítulo. 


O uso de ferramentas quantitativas que ofereçam subsídios objetivos para a atividade de 
gerenciamento foi o enfoque que destacamos neste capítulo. Nesse contexto, a PO apresen- 
ta-se como uma ciência prática. Ela estabelece parâmetros decisórios confiáveis, nos quais 
considera fatores e cenários dos problemas e propicia, por meio de modelos matemáticos, 
a possibilidade de visualizarmos as possíveis soluções para situações cujas variáveis, restri- 
ções e funções objetivo são lidas a partir de cálculos, os quais são modelados a partir das fases 
de estruturação do problema. Por isso, também é necessário que você tenha clareza quanto 
aos conceitos de método, sistema e solução ótima, examinados neste estudo, de forma a po- 
der aproveitar esse instrumental nos processos de tomada de decisão. 
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Questoes para revisao 


1. Em relação às fases da PO, assinale (V) para verdadeiro e (F) para falso: 


Construção e alteração do modelo. 
Estabelecimento e controle das soluções. 


Cálculo do modelo. 


() 

() 

() 

( ) Consulta aos operários. 
( ) Teste do modelo e da solução. 
( ) Compra de software. 

( ) Formulação do problema. 

() 


Implantação e acompanhamento. 
2. Em relação aos passos para se analisar um problema de PO, é preciso: 


I) analisar os investimentos financeiros. 
II) verificar as variáveis de decisão. 
III) levantar as restrições e fazer a definição da função objetivo. 


IV) identificar os parâmetros. 
Assinale a alternativa verdadeira: 


a) Apenas as afirmativas I e III estão corretas. 
b) Apenas as afirmativas I e IV estão corretas. 
c) Todas as afirmativas estão corretas. 


d) Apenas a afirmativa IV está incorreta. 

3. Assinale a alternativa correspondente aos fatores que interferem na tomada de decisão: 
a) Importância; agentes; risco; ambiente; conflitos. 

b) Importância; diagnóstico; risco; ambiente; conflitos. 

c) Importância; agentes; risco; legislação; conflitos. 

d) Importância; agentes; fator humano; ambiente; conflitos. 


4. Como já comentamos, na opinião de Lachtermacher (2009), a postura de abandono da in- 
tuição no processo de tomada de decisão não é a melhor escolha; o recomendável é fazer a 


união das duas fontes de informações. Qual a sua opinião sobre isso? 


5. O que são os problemas enfrentados pelas organizações? 


23 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal 


‘leuad o6ipo5 op p8, "ue ojad opiund e 8664/0196 U 127 eu opi2ejeqeise awya e sieJoyne sojlaJip sop oeSejol^ y "seJeqesgJelu| &Jojip3 ep oe5ezuojne ejagud e was euo) no oraw Jenbjenb Jod eprznpoudej Jes ejepod og5eotjqnd ejsep aued eunuueN 


capítulo 


A programação linear no contexto 
da pesquisa operacional 
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Conteüdos do capítulo 


* Como utilizar a programação linear na resolução de problemas de pesquisa operacional. 
e Formulação de um problema matemático padrão em progamação linear com objetivos e 


restrições. 


Após o estudo deste capítulo, você será capaz de: 


1. fazer a modelagem matemática de um problema na busca pela solução ótima; 
2. aplicar a programação linear como ferramenta de modelagem matemática na pesquisa 


operacional. 


A programação linear (PL), segundo Passos (2008), é uma técnica de otimização aplicada em 
sistemas de equações e inequações lineares representativos de modelos previamente elaborados. 


Podemos dizer, então, que é uma aplicação matemática utilizada por profissionais para o 
melhor aproveitamento possível da produção, de modo a evitar desperdícios de produtos, 
matérias-primas ou mão de obra, consistindo na modelagem e solução de problemas de oti- 


mização de uma função linear, diante de restrições também lineares. 


2.1 Aplicação da função objetivo e das restrições em PL 


Como vimos no capítulo anterior, quando falamos em problemas de otimização, referimo- 
-nos àqueles em que queremos maximizar (aumentar) ou minimizar (diminuir) uma função 
(podendo ser uma função de receita, custo, demanda, produção, entre outras). Nessas situa- 
ções, precisamos levantar a função objetivo e também as restrições que o sistema analisado 


nos apresenta. 


A função objetivo em PL 


Você pode elaborar a definição da função objetivo, por exemplo, no sentido de maximizar 


o lucro ou minimizar o custo. 


E importante ressaltar que, nos problemas de PL, a função objetivo e as restrições são sempre 
equações ou inequações lineares. Nesse contexto matemático, os problemas ficam com sua 


linguagem modificada, ou seja, passam para a linguagem matemática. 
Desse modo, a função objetivo pode ser escrita nas duas formas que apresentamos a seguir: 


1. Se o problema for maximizar z: 


max Zz=Cx, E ex, Cd tata CX 
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2. Se o problema for minimizar z: 


mın z = cx, ar CX, pieno CX 


Observação: Esse é um aspecto matemático utilizado em PL que pode causar estranheza aos 
não matemáticos, pois a minimização de uma função z é equivalente à maximização dessa 
função em sua versão negativa —z. Como vimos nos itens 1 e 2, ambas são somatórias. 


Em ambos os casos, c, C, ..., c, são números reais, e x,, X, ..., x, são as variáveis do problema. 
1^2 n 1. 2 n 
Nos problemas de PL, além da função objetivo, que consiste em uma expressão matemática 
que representa a meta a ser alcançada, temos também as restrições. 
As restrições em PL 


Para você entender melhor esse conceito, imagine uma indústria de laticínios que deseja oti- 
mizar a sua produção, maximizando o lucro. Nesse caso, as limitações de ordem prática en- 
contradas para fazer a otimização da produção (função objetivo) constituem as restrições do 
problema em PL, que são: 


e quantidade disponível de matéria-prima; 
e capacidade do setor produtivo; 

e mão de obra; 

e limitações no preço. 


Em outros problemas de PL podem existir diferentes limitações, tais como limitações de lo- 
calidade ou de espaço físico e distância entre localizações. Nesse contexto, podemos encon- 
trar, por exemplo: 


e um agricultor que deseja plantar diversas culturas, mas que tem um limite de espaço 


a ser cultivado, ou seja, uma restrição espacial; 


e um investidor que deseja diversificar suas aplicações, mas que possui apenas certa 
quantia a ser aplicada, isto é, tem uma restrição de capital; 


e uma transportadora que tem como objetivo otimizar as entregas, mas que pode contar 
apenas com determinado número de veículos, o que significa uma restrição de quan- 
tidade disponível (de veículos). 


No cenário das restrições, encontramos as restrições de igualdade ou as de desigualdade. Veja 


como são representadas nas equações e inequações a seguir. 


1. Restrição de igualdade (equação): 


ane ra ca a Di 
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2. Restrições de desigualdade (inequações): 


AX; 


ou 


AX; 


Para Sb: 


PES oos PASS Elo 


2? 1 


Forma geral ou padrão 


A forma geral de um problema de PL é caracterizada pela padronização com o objetivo de 


facilitar o entendimento. Veja como se estrutura a forma geral ou padrão: 


max z= CX, +C,X, Tuc CX 


sa. axtaXt+..ta x € b, 
da A V ra X SD. 
à X, tà Xt. LX Eb. 
x,20,x,20,.., x, 20 

Em que: 


os termos a,, b, c são os coeficientes das equações e inequações que descrevem no proble- 
ma os números de quantidade, valor e custos (considerei=1,2,3,..,mej=1,2,3,...,n). 
Você se lembra das operações com matrizes em que esses termos indicam a posição dos 

elementos de uma matriz? Se a resposta for “não”, veja o Capítulo 8, no qual rememo- 


ramos alguns aspectos básicos das operações matriciais; 


as variáveis X, X, ..., X, são escolhidas de forma que as restrições sejam satisfeitas e 
a função objetivo otimizada; 


o termo “s.a.” significa “sujeito a”; 


já os termos b, b,, ..., b são chamados parâmetros da função, nos quais se representam 
as limitações das restrições; 


as restrições x, 2 0, x, 2 0, ..., x, 2 0 são chamadas de restrições de não negatividade. Elas 


ocorrem porque não podemos ter quantidades negativas de produtos ou de recursos. 


Com o intuito de ilustrar melhor as ideias apresentadas, vamos detalhar um exemplo de pro- 
blema de PL. 
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Com a utilização da técnica de PL, aplicaremos os conteúdos de PO em uma indústria (ima- 


ginária) de brinquedos. 


Vamos supor que essa indüstria fabrica dois modelos de veículos em miniatura: caminhone- 
tes e esportivos. Para a fabricação dos brinquedos, as principais matérias-primas empregadas 
são o plástico e o alumínio. Nesse processo de fabricação ocorre o seguinte: 


e O modelo esportivo consome 400 g de plástico e 300 g de alumínio. 
e A caminhonete requer 700 g de plástico e 150 g de alumínio. 


e A disponibilidade mensal de plástico é de 1 t e a disponibilidade mensal de alumínio 
é de 600 kg. 


e O lucro unitário referente à caminhonete é de R$ 12,00 e o lucro unitário referente ao 
modelo esportivo é de R$ 15,00. 


Suponha que toda a produção de miniaturas esportivas seja vendida e que a empresa consiga 
vender, por mês, no máximo 700 caminhonetes. Diante dessa situação, determine a quantida- 
de de cada modelo que deve ser produzida de maneira que o lucro seja máximo. 


Lembre-se de que, para formular o problema, é preciso identificar quais são as variáveis, a 
função objetivo e as restrições. 


Formulando a função objetivo 


Nesse problema, as variáveis são as quantidades de cada modelo a serem produzidas. Dizemos, 


então, que: 
e x, éa quantidade de caminhonetes; 
e x, éa quantidade de modelos esportivos. 


Com as informações transpostas para a linguagem matemática, podemos facilmente obter 


a função objetivo. Vejamos: como o lucro unitário referente à caminhonete é de R$ 12,00 e o 
lucro unitário referente ao modelo esportivo é de R$ 15,00, para obtermos o lucro total, de- 


notado por z, basta multiplicarmos os lucros unitários pelas respectivas quantidades que se- 
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rão produzidas, ou seja: 

z = 12x, + 15x, 

Como o objetivo é obter o lucro máximo, a função objetivo é dada por: 
max z = 12x, + 15x, 


É importante ressaltar que o lucro unitário é a diferença entre o preço de venda praticado pe- 
la empresa e o custo de produção de cada item. 
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Formulando as restrições 


Para determinar quais são as restrições do problema, primeiro você precisa verificar quais são 

os fatores que limitam a produção. No caso dessa fábrica de brinquedos, as restrições se re- 
ferem às quantidades disponíveis de plástico e de alumínio. Há também uma restrição em 

relação à quantidade máxima de caminhonetes que poderá ser comercializada. O número 

de restrições para esse problema é, portanto, igual a três. 


Na formulação matemática, podemos escrever a primeira restrição, a referente à quantidade 


de plástico que será consumida, como: 
0,7x, + 0,4x, < 1.000 


Nesse cômputo, temos a quantidade consumida de plástico na fabricação de caminhonetes, 
pois sabemos que a produção de uma caminhonete de brinquedo requer 700 g de plástico 
(0,7 kg), enquanto cada modelo esportivo requer 400 g de plástico (0,4 kg). Para sabermos o 
total de plástico que será utilizado na produção, basta multiplicar 0,7 por x, e 0,4 por x, e, em 
seguida, somar essas quantias, obtendo a expressão: 


0,7x, + 0,4x, 


Como a quantidade máxima de plástico que a indústria tem disponível é 1 tonelada (1.000 kg), 
a soma 0,7x, + 0,4x, não pode ultrapassar essa quantidade. Por esse motivo, escrevemos que 
0,7x, + 0,4x, tem de ser menor ou igual a 1.000, ou seja: 


0,7x, + 0,4x, € 1.000 


Da mesma maneira, podemos obter a segunda restrição, a referente ao consumo de alumí- 
nio. Como a caminhonete requer 150 g de alumínio (0,15 kg) e o modelo esportivo requer 
300 g (0,3 kg) de alumínio, temos que o total de alumínio que será consumido na produção 
dos modelos é: 


0,15x, + 0,3x, 

Sabendo que a disponibilidade mensal de alumínio é 600 kg, a segunda restrição fica assim: 
0,15x, + 0,2x, < 600 

Finalmente, a terceira restrição, relacionada à produção máxima de caminhonetes, é dada por: 
x, < 700 

Portanto, a formulação do problema de PL proposto a que chegamos é dada como: 


max z = 12x, + 15x, 


s.a. 0,7x, + 0,4x, < 1.000 
0,15x, + 0,3x, < 600 
x <700 
x,20,x,20 
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Observação: Para resolvermos um problema de PL, caso existam, as restrições > devem ser 
trocadas por x. Para isso, basta multiplicarmos cada restrição de > por (-1) e, em seguida, in- 
vertermos a desigualdade > por x. Por exemplo: 


A restricáo: 
iX $a x tum ox b 
é equivalente à restricáo: 


x €-b 


“AX 7 AX, e AX 1 


Como já vimos, as restrições x, >0,x,>0,... x, 2 0 são chamadas de restrições de não negativida- 
de. Elas ocorrem porque não podemos ter quantidades negativas de produtos ou de recursos. 


Vamos formular um exemplo de fácil percepção para a formulação do modelo matemático 
chamado de forma-geral ou forma-padrão. 


Vamos considerar na formulação matemática do modelo-padrão uma empresa de artefatos 


de madeira que fabrica porta-trecos e potes de madeira. 
Ocorre que nessa produção a empresa gasta: 
e duas unidades de madeira para fabricar uma unidade de porta-trecos; 
e duas unidades de madeira para fabricar uma unidade de pote de madeira. 


Sabemos que o total disponível de madeira é de seis unidades e que o lucro unitário por por- 
ta-trecos de madeira é de R$ 2,80, enquanto o lucro unitário por pote é de R$ 5,00. 


Diante dessa situação, quanto devemos produzir de porta-trecos ou de potes de madeira 
de modo a maximizar o lucro por hora da empresa? 


Resolução 


Primeiro passo - Identificamos as variáveis que queremos encontrar, ou seja, as variáveis que 


nos permitem alcançar o objetivo do problema, nomeando-as: 


e x quantidade produzida de porta-trecos; 


e x,- quantidade produzida de potes de madeira. 
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(| | Exercício resol 


Segundo passo — Procuramos identificar a função objetivo do problema. 


Se formos maximizar, usamos a simbologia z, Ou max z, lembrando que essa simbologia 
deve ser seguida do sinal de igual por ser uma equação. Assim, na montagem da função objeti- 
vo, somamos a multiplicação das variáveis identificadas pelos respectivos valores das receitas. 


Temos a seguinte informação no problema: “O lucro unitário por porta-trecos 
de madeira é de R$ 2,80 e o lucro por pote é de R$ 5,00”. Logo, podemos escrever: 


max z = 2,80x, + 5x, 


Terceiro passo — Formulamos agora as restrições do problema, identificando-as como as limi- 
tações desse problema. Podemos observar que temos a restrição de matéria-prima, pois a in- 
formação é a de que a empresa gasta duas unidades de madeira para fabricar uma unidade de 
porta-trecos e duas unidades de madeira para fabricar uma unidade de pote de madeira, sen- 


do que o total disponível de madeira é de seis unidades. Logo, nossa restrição se resume em: 
2x,+2x,<6 
Portanto, a formulação do problema de PL proposto é a seguinte: 
max z = 2,80x, + 5x, 
s.a 2x,+2x,< 6 
x,x,20 


Chegamos, assim, ao formato que chamamos de forma-padráo de um problema na PL. 


Um dos métodos mais famosos utilizados para se obter a solução de um problema de PL é 
o método simplex, o qual veremos no próximo capítulo. Inicialmente, podemos empregar 
o método gráfico para resolver problemas mais simples de PL (com duas variáveis e algu- 
mas restrições). 


Agora vamos praticar um pouco, colocando os problemas propostos na sequéncia no 


formato-padrão. 


Uma empresa, depois de uma estruturação no processo de produção, ficou com disponibili- 
dade de recursos. Após uma apuração, constatamos que tinha disponíveis três tipos de recur- 
sos, que chamaremos de R, R, e R, e que eles possibilitariam a fabricação de dois produtos, 
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que chamaremos de P e P,. Sabemos que o lucro obtido com o produto P, seria de 150 u.m.' 
e que P, daria 174 u.m de lucro. 


A tabela a seguir representa o uso dos recursos utilizados pelo setor de produção. 


Produto R,porunidade | R,porunidade | R, por unidade 
P, 2 4 6 
P, 6 2 4 
Disponibilidade 100 80 150 
de recursos 


Nessa situação, a questão que enfrentamos é saber que quantidade de produtos P, e P, deve 
ser fabricada para que a empresa tenha maior lucro. 


Resolução 
max z = 150 P, + 174 P, 


sa. 2P, + 6P, < 100 


ireitos autorais é crime estabelecido na Lei n° 9.610/1998 e punido pelo art. 184 do Código Penal. 


4P * 2P, « 80 


6P, + 4P, < 150 


2.3 Interpretação geométrica e solução gráfica 


Como todo problema de PL apresenta restrições lineares, podemos representar as restrições 
de um problema de duas variáveis em um sistema de eixos coordenados denominado pla- 


no cartesiano. 


Vamos relembrar: um par ordenado (x, y) representa um ponto no plano, como você já apren- 
deu em geometria analítica. 
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Figura 2.1 — Pontos coordenados no plano cartesiano 
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Assim, utilizando seus conhecimentos matemáticos, você pode dispor os dados em um grá- 
fico para construir um modelo. 


Para representarmos nesse sistema de eixos coordenados cada restrição de desigualdade, es- 
crevemos uma equação associada a cada uma delas, contendo os mesmos coeficientes e o mes- 
mo termo independente. No caso do exemplo anterior, aquele da indústria de brinquedos 
que fabrica dois modelos de veículos em miniatura (caminhonetes e esportivos), a restrição 
0,7x, + 0,4x, < 1.000 está associada à equação 0/7x, + 0,4x, = 1.000. Da mesma forma, a restri- 
ção 0,15x, + 0,3x, < 600 está associada à equação 0,15x, + 0,3x, = 600, e a restrição x, < 700 está 


associada à equação x, = 700. 
Mas como obter a representação da primeira restrição? 


Para obter a representação da primeira restrição, basta você construir uma pequena tabela 
atribuindo dois valores quaisquer para a variável x,. Em seguida, você deve calcular os res- 


pectivos valores da variável x,. 
E de que maneira fazemos esse cálculo dos valores das variáveis? 


Esses pontos da forma (x, ex,) são chamados de pares ordenados. Uma maneira simplificada 
de obtermos os pares ordenados é atribuirmos o valor zero à variável x, para, dessa maneira, 
encontrarmos o respectivo valor de x,. Em seguida, atribuímos à variável x, o valor zero pa- 
ra encontrarmos o valor de x,. 


A vantagem de utilizarmos esse artifício é que os pontos obtidos estarão sobre os eixos coor- 
denados, facilitando o processo de representação gráfica de cada uma das restrições. Como 
as restrições são de desigualdade, é importante ainda considerar em qual semiplano a região 
factível estará. Um método bastante prático é considerar a origem do plano cartesiano, o par 
ordenado (0, 0). Se, ao substituirmos os valores de x, e x, por zero, a desigualdade for ver- 


dadeira, a região factível estará no semiplano que contém a origem. Caso contrário, a região 
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factível estará no semiplano oposto, ou seja, o que nào contém a origem. Esse processo é usa- 


do também para as demais restrições. 
Como desenvolver o cálculo desse processo? 
Fazendo: 

e x,=0, temos que x, = 2.500; 

e x,=0, temos que x, = 1.428,57. 


Assim, temos a tabela com os valores de x, e x, e os respectivos pares ordenados: 


X, x, ex x) 
0 2.500 (0; 2.500) 
1428,57 0 (1.428,57; 0) 


Logo: 
e o ponto (0; 2.500) está localizado no eixo da variável x, pois x, é igual a zero; 
e oponto (1.428,57; 0) fica sobre o eixo x, pois x, é igual a zero. 


Assim, como o ponto (0, 0), origem do sistema, satisfaz a inequação 0,7x, + 0,4x, < 1.000, te- 
mos que a região factível está no semiplano que contém a origem. 


Figura 2.2 - Desenho da primeira restrição 


2.500 0,7x,+ 0,4x, < 1.000 


4 


1.428,57 


No mesmo plano cartesiano em que representamos a primeira restrição, representaremos 
também a segunda e a terceira restrição. Para representarmos a segunda restrição, faremos 
x, = 0. Nesse caso, x, = 2.000. 


Observe que, se fizermos x, = 0, teremos x, = 4.000, o que resultará em um ponto muito dis- 


tante da região factível. Nesse caso, apenas para fins de representação geométrica, faremos 
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x, = 1.000. Substituindo esse valor na equação 0,15x, + 0,3x, = 600, temos que x, = 1.500, resul- 
tando no par ordenado (1.000; 1.500). 


Portanto, a tabela contendo os valores de x, e x, e os respectivos pares ordenados fica assim: 


X, X, (xy Xj) 
0 2.000 (0; 2.000) 
1.000 1.500 (1.000; 1.500) 


Considerando esses valores, você pode observar que a região factível está localizada abaixo 
da segunda restrição, pois a origem (0, 0) satisfaz a restrição 0,15x, + 0,3x, = 600. 


Figura 2.3 - Desenho da primeira e da segunda restrição 
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E, finalmente, a terceira restrição é uma reta vertical paralela ao eixo x, que passa pelo ponto 
x, = 700. É fácil perceber que isso ocorre, pois, considerando-se a equação x, = 700, para qual- 
quer valor de x,, x, sempre é igual a 700. 


Nesse cálculo, como a origem (0, 0) satisfaz a desigualdade x, < 700, temos que a região fac- 
tível está localizada à esquerda da terceira restrição. 


Figura 2.4 — Desenho com todas as restrições (total de três) 


x, ^ 
2.500 07x, + 4x, « 1.000 ! 
T. 
a al 0,15x, + 0,3x, < 600 
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v 


v 


< 1.000 1.428,57 x 
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Observe que em cada figura com suas respectivas restrições há uma seta indicando em que 
sentido a desigualdade é satisfeita. Considerando as limitações impostas pelas restrições do 
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problema, obtemos uma região denominada região factível (ou região viável). Qualquer ponto 


dessa região satisfaz as restrições e pode ser uma solução do problema. 
Figura 2.5 - Região factível 


X, 4 


Região 


factível 


E qual é a solução ótima? 


Nesses cálculos, para encontrar a solução ótima, você precisa considerar também a função 
objetivo. 


Nesse estágio do processo, é importante lembrar que a função objetivo do exemplo é dada 
por z = 12x, + 15x, bem como o fato de que é possível representar graficamente a função ob- 
jetivo atribuindo-se valores aleatórios para z. 


Considerando essas condições, escolhemos, então, ao acaso, alguns valores para z: 


z-0 
z = 15.000 
z = 30.000 


Dessa maneira, é possível obtermos, respectivamente, as seguintes equações: 
12x, + 15x,=0 

12x, + 15x,= 15.000 

12x, + 15x,= 30.000 

Graficamente representamos essas equações por retas pontilhadas. 


O processo para representar graficamente as equações segue o mesmo princípio utilizado na 
representação das restrições. Vamos retratar a equação 12x, + 15x,= 0 construindo uma tabe- 
la com dois valores atribuídos aleatoriamente para x, e calcular x,: 


37 


5 
E 
: 
5 
F 
E 
8 
: 
E 
3 
5 
E 
& 
E 
E 
: 
5 
P 
3 
E 
E 
i 
E! 
E 
E 
[4 
E 
e 
$ 
8 
8 
E 
e 
- 
£ 
E 
L4 
E 
e 
£ 
E 
E 
8 
& 
H 
E 
2 
E 
e 
8 
: 
e 
E 
E 
5 
B 
B 
8 
* 
E 
H 
E 


Nenhuma parte desta 


X, X, (x x) 
0 0 (0; 0) 
500 -400 (500; -400) 


A equação 12x, + 15x,= 15.000 resulta na seguinte tabela: 


X, X, (x x) 
0 1.000 (0; 1.000) 
1.250 0 (1.250; 0) 


E a tabela construída a partir da equação 12x, + 15x, - 30.000 é: 


X, X, (Xy xj) 
0 2.000 (0; 2.000) 
2.500 0 (2.500; 0) 


Representando graficamente as três equações, temos a figura a seguir. 


Figura 2.6 — Representação gráfica de um problema de PL 


X, A 
~ 
z=30.000 
2.000" . 
z 15.000 im mm | 
De Região EA l Y 

factível ER - 

z-0 E. T é 
TER 500 < 1250 ^ B 2.500 X 

v K 


-400 T 


Dessa forma, você pode perceber que, à medida que o valor de z aumenta, as retas que repre- 
sentam a função objetivo ficam cada vez mais próximas do vértice V. Vale a pena lembrar que 
toda figura em um plano cartesiano sempre apresentará vértices que indicam seus cantos. Em 
PO, o vértice que aumenta o valor da função objetivo é chamado de gradiente. 


Assim, podemos concluir que a solução do problema pode ser obtida fazendo-se a interseção 


das retas que se interceptam nesse vértice. São elas: 


0,7x, +0,4x,=1.000 e 0,15x,+0,3x, = 600 


Resolvendo o sistema | 
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Obtemos a solução x, = 400 e x, = 1.800, a qual também pode ser encontrada por meio de mé- 
todos algébricos de sistemas de equações. 


Substituindo esses valores na função objetivo, temos que o lucro máximo é z = 31.800. 


Talvez você esteja pensando: o que faço com tanta informação? Por isso vamos retomar essa 
situação (processo) mais detalhadamente com a resolução de um novo exemplo. 


Nesse contexto, vamos partir do pressuposto de que você conseguiu deixar seu problema na 
forma-padrão, como temos a seguir, a qual se refere à maximização do lucro de certo produ- 
to, obedecendo a certas restrições: 


max Z = 2x *5x, 
s.a x <6 
x, $5 


x *2x, < 14 


X; X, 20 
Para você resolver graficamente essa situação, o primeiro passo é igualar as restrições, trans- 


formando-as em equações. Para isso, basta substituir o valor das desigualdades (<, 2) por 
uma igualdade, ou seja: 


x,76 
xX,-5 
x, + 2x,- 14 


Devemos agora traçar as retas das referidas equações no plano cartesiano. Para isso, basta 
igualar uma variável de cada restrição a zero, achando, assim, os pontos coordenados. 


A abordagem que fizemos da PL no contexto da PO teve como intento levantar as fases de ela- 
boração de uma PO, isto é, as etapas correspondentes à definição do problema e à construção, 
reconstrução, solução e validação do modelo. Nesse contexto, apresentamos e trabalhamos 
com a PL por ser uma ferramenta, uma técnica de planejamento que possibilita, por meio da 
linguagem matemática (equações, inequações, interpretação geométrica e solução gráfica), 
encontrar o modelo para a situação real, considerando-se o lucro máximo ou o custo mínimo. 
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Nenhuma parte desta 
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"t Questoes para revisao 


1. Resolva graficamente o seguinte problema de PL: 


max  z-3x,*2x, 


s.a. —x, + 2x, <4 
-3x, + X, <6 
X,* X, € 10 


2. Resolva graficamente o problema de PL a seguir: 


max z= 4x, t 3x, 


s.a. 2x, <9 
x, S7 

E sx SE 

X, X 240 


3. Suponha que uma fábrica de implementos agrícolas produza os modelos A e B, os quais 
proporcionam lucros unitários da ordem de R$ 16,00 e R$ 30,00, respectivamente. A exigên- 
cia de produção mínima mensal é de 20 unidades para o modelo A e de 120 para o modelo B. 
Cada tipo de implemento requer certa quantidade de tempo para a fabricação das partes que 
os compõem, para a montagem e para os testes de qualidade. Especificamente, uma dúzia de 
unidades do modelo A requer 3 horas para fabricar, 4 horas para montar e 1 hora para testar. 
Você ainda constata que uma dúzia de unidades do modelo B requer 3,5 horas para fabricar, 
5 horas para montar e 1,5 hora para testar. Agora, considere que, durante o próximo mês, a 
fábrica tem disponíveis 120 horas de tempo de fabricação, 160 horas de montagem e 48 horas 
de testes de qualidade. Formule o problema de programação de produção como um modelo 
de PL (forma-padrão) e resolva-o graficamente. 


4. Com base no que foi apresentado, você poderia formular uma finalidade para a função ob- 
jetivo no processo de PO? Ainda em relação à função objetivo, você concorda que, ao atribuir- 


mos valores aleatórios para z, seja possível representá-la graficamente? 


5. Para que serve a forma-padrão? 
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O método simplex e variantes 
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Conteüdos do capítulo 


e Conceituação e características do método simplex. 
e Passos para resolver um problema de programação linear utilizando o método simplex. 
e Utilização de variáveis artificiais na busca pela solução em programação linear. 


e Conceituação e uso do método Big-M. 


Após o estudo deste capítulo, vocé será capaz de: 


1. identificar e definir as variáveis de um problema de programação linear; 
2. resolver um problema de programação linear com o método simplex; 
3. introduzir variáveis artificiais na resolução de problemas; 


4. resolver um problema de programação linear pelo método Big-M. 


No capítulo anterior, vimos a forma geral de um problema de programação linear (PL) e a sua 
interpretação gráfica. Neste capítulo, abordaremos um importante método amplamente utilizado 
na resolução de problemas lineares de otimização: o método simplex, desenvolvido em 1947 por 
George B. Dantzig (1914-2005), professor emérito de pesquisa operacional e ciência da computação. 


3.1 Como e por que utilizar o método simplex? 


O método simplex é uma importante ferramenta destinada a resolver problemas de PL. Esse 
método consiste em buscar, caso existam, uma ou mais soluções partindo-se de uma solução 
básica factível, gerando uma sequência de soluções factíveis. Quando essa sequência é com- 
pletada, a solução ótima é obtida. 


Nesse sentido, o professor Lisboa (2002, p. 15) afirma que esse método “caminha pelos vér- 
tices da região viável até encontrar uma solução que não possua soluções vizinhas melhores 
que ela. Esta é a solução ótima”. 


Você deve observar, e sobre isso o autor ao qual nos referimos também alerta, que há duas 
situações em que não é possível chegarmos à solução ótima. Esses casos ocorrem quando: 


1. não há uma solução que devamos ou que possamos executar por apresentar restrições 


de incompatibilidade; 


2. não encontramos um máximo ou um mínimo. Nesse caso, ocorre que uma das variá- 
veis pode se estender ao infinito (apresentar essa tendência), embora as restrições se- 


jam satisfeitas; consequentemente, teremos um valor sem limites para a função objetivo. 


E importante conhecer esse método justamente porque ele apresenta, excetuando-se as duas 
situações descritas anteriormente, a possibilidade de você resolver, por meio de um esquema 
de equações lineares, o modelo de PL. 
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Passos para resolver um problema de PL utilizando o método simplex 


E necessário dispor de uma sistemática que determine diretrizes para um procedimento ob- 


jetivo, como: 
e definir qual sistema de equações devemos usar na resolução; 


e indicar o próximo sistema resolutivo que devemos utilizar para encontrarmos uma 
solução melhor; 


e identificar uma solução ótima, quando a encontrarmos. 


Como já dissemos, você pode empregar o método simplex para resolver um problema de PL 
após deixá-lo na forma-padrão. Ele é a sistemática que atende a essa necessidade. No entan- 
to, nesse processo, há uma sequência de passos que devemos observar para chegarmos à so- 
lução, a saber: 


1. acrescentar as variáveis de folga ou excesso no modelo de PL (lembrando que o sinal < 


acrescenta a variável de folga; se o sinal for >, diminui a variável de excesso); 


2. montar a tabela com os coeficientes da função objetivo, os coeficientes das restrições, 
os termos independentes; 


3. identificar quais são as variáveis básicas e quais são as não básicas, utilizando, para is- 


so, os recursos matemáticos apresentados nos exemplos e nos modelos; 


4. verificar qual variável entrará na base (basta escolher o coeficiente mais negativo da 
função objetivo); 


5. verificar qual variável irá sair da base (divida cada termo independente pelos respecti- 
vos coeficientes positivos da coluna pivô, referente à variável que irá entrar na base, e 


considere o menor valor obtido); 


6. identificar o elemento pivô (interseção da coluna da variável que irá entrar na base com 


a linha da variável que irá sair da base); 
7. seo pivô for diferente de 1, multiplicar a linha pelo inverso multiplicativo do atual pivô; 


8. utilizando operações elementares (multiplicação de linha por um escalar e soma de li- 
nhas), exceto o pivô, transformar em zero os elementos da coluna pivô da variável que 
está entrando na base, colocando os resultados obtidos em um nova tabela; 


9. enquanto houver coeficientes negativos na linha referente à função objetivo, repetir os 


passos do terceiro ao oitavo. 
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Perguntas & respostas 


O que são variáveis básicas e não básicas? 


As variáveis básicas são aquelas cujos valores são calculados pelo sistema de equações. Já 


as variáveis não básicas são aquelas zeradas. 


E as variáveis de folga ou excesso? 


São aquelas que acrescentamos no problema, para poder resolvê-lo pelo método proposto. 


Você pode representar essa operacionalização expressa nos nove passos por meio de uma ta- 


bela como a apresentada a seguir. 


Tabela 3.1 — Representação genérica da tabela utilizada no método simplex 


Z C E; C. 
Xy un f fin 
x ü ü ü 


b2 21 


ü 


bm ml 


Ago 


mn 


E importante atentarmos para o fato de que, quando dizemos representação genérica, quere- 


mos constituir uma representação para qualquer problema, de forma generalizada, como se 


fosse um modelo a ser seguido. 


Para ilustrarmos esses procedimentos de resolução de um problema de PL por meio do mé- 


todo simplex, utilizaremos como exemplo a mesma indústria de brinquedos apresentada no 


Capítulo 2 para ilustrar um cálculo de PL. Lembra-se? A situação era esta: 


e A empresa fabrica dois modelos de veículos em miniatura: caminhonetes e esportivos. 


e As principais matérias-primas empregadas para a fabricação dos brinquedos são plás- 


tico e alumínio. 


e O modelo esportivo consome 400 g de plástico e 300 g de alumínio. 


e A caminhonete requer 700 g de plástico e 150 g de alumínio. 


e A disponibilidade mensal de plástico é de 1 t e a disponibilidade mensal de alumínio 


é de 600 kg. 


e Olucro unitário referente à caminhonete é de R$ 12,00. 


e O lucro unitário referente ao modelo esportivo é de R$ 15,00. 
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Agora que recordamos as condições em que a fábrica opera, vamos supor que, por més, toda 
a produção de miniaturas esportivas seja vendida, enquanto a venda de caminhonetes seja 
no máximo de 700 caminhonetes. Com esses dados, determine a quantidade de cada modelo 


que deve ser produzida de maneira que o lucro seja máximo. 


No capítulo anterior, você viu que a formulação do modelo-padrão de PL para o problema 
era a seguinte: 


max z= 12x, + 15x, 


s.a. 0,7x, + 0,4x, « 1.000 
0,15x, + 0,3x, < 600 
x « 700 


1 


x, 20,x, 20 


Como é visivelmente perceptível que as restrições são todas com sinal de desigualdade (me- 
nor igual), ou seja, inequações, devemos, nesse caso, acrescentar novas variáveis de folga a 


cada restrição. Assim, temos: 


max z= 12x, + 15x, 


s.a. 0,7x, + 0,4x, + x, = 1.000 
0,15x, + 0,3x, pA = 600 
a +x,=700 


x, 20,x, 20 


Quais sáo os procedimentos que adotamos para realizar esses cálculos? 


Vamos apresentar passo a passo o processo pelo qual realizamos esses cálculos. Para construir 


a primeira tabela, utilizamos o modelo representado a seguir: 


Figura 3.1 - Modelo de tabela para o uso do método simplex 


Função - z m. 
s Coeficientes da função objetivo 
objetivo 
Termos e - 
. Coeficientes das restrições 
independentes 


Sabemos que toda tabela apresenta linhas e colunas (mesmo que não fiquem evidenciadas 
pelo traço da linha) e denominamos as linhas, para identificá-las, como linha um, linha dois etc. 
e suas colunas como primeira coluna, segunda coluna etc. 


Colocamos na linha zero (primeira linha da tabela) o valor de z (que sabemos ser a função 
objetivo e cuja solução inicial tem valor de z = 0) e os coeficientes da função objetivo, todos 
com os sinais trocados. Alguns autores preferem colocar os valores da função objetivo na úl- 
tima linha da tabela, que não altera a solução ou os procedimentos de resolução, é apenas 
uma formalização diferenciada feita por pessoas diferentes. Nós utilizaremos a primeira linha. 


No exemplo, os coeficientes a serem colocados são: 12, -15, 0, 0 e 0. 
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Vocé deve ficar atento para o fato de que completamos os campos vazios sempre com zero para 
identificar os campos nulos. Assim, temos o seguinte: 


e Na linha um colocamos a quantidade de recursos disponíveis referente à primeira res- 
trição. Lembre-se de que são somente os coeficientes: 0,7; 0,4; 1. O valor de 1.000 entrará 
na primeira coluna, pois é o coeficiente do termo independente ou do recurso disponível. 


e Na linha dois escrevemos os seguintes valores: 15; 0,3; 0; 1 e 600, que são, respectiva- 
mente, a disponibilidade de recursos para a segunda restrição e os seus coeficientes. 


e Finalmente, na linha três registramos os valores da última restrição: 1; 0; 0; 0; 1e 700, 


como a disponibilidade de recursos. Com esses procedimentos, o nosso quadro fica assim: 


1 2 3 4 5 

z=0 -12 -15 0 0 0 
X 1.000 0,7 04 1 0 0 
i 600 15 0,3 0 1 0 
xS 700 1 0 0 0 1 


Nesse processo de resolução, inicialmente, as variáveis básicas são as variáveis de folga: x, x, 


e x.. As variáveis não básicas são x, e x, ou seja, as variáveis originais do problema. 


Lembre-se de que, para decidirmos qual variável entrará na base, precisamos verificar qual 
variável não básica irá fornecer o maior lucro possível. Para isso, basta verificar qual variável 
tem o coeficiente mais negativo na linha zero. Nesse caso, a variável que entra na base é x, cujo 
coeficiente na linha zero é igual a -15. O próximo passo é determinar qual variável básica irá 
sair da base. Para isso, iremos dividir cada termo independente pelos respectivos coeficientes 
positivos (valores negativos ou nulos não são considerados) da coluna referente à variável x,: 


1.000 : 0,4 = 2.500 


600 : 0,3 = 2.000 
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Note que a divisão de 700 por 0 não é possível. Como o menor resultado obtido é 2.000 (di- 
visão de 600 por 0,3), a variável que irá sair (—) da base é x, pois é a variável básica que tem 
o elemento unitário na linha referente a essa divisão. 


Escolhemos aqui o menor valor da divisão em virtude das limitações impostas inicialmente 


Primeira iteração — resolução na Tabela 1 


pelo problema. Logo, a primeira tabela fica assim: 


x, X; X. x, X; 
z=0 -12 -15 0 0 0 
1.000 0,7 0,4 1 0 0 
600 0,15 0 1 0 
700 1 0 0 0 1 


Variáveis básicas: x, x, e x, 


Variáveis nào básicas: x, e x, 


Entra: x, 


Observe que o pivô (elemento da interseção da coluna da variável que entra na base com a li- 
nha da variável que sai), nesse caso, é igual a 0,3. Precisamos transformar o pivó em 1. Como 


3 
03-2 10' basta multiplicar todos os termos da linha dois por 2 (inverso multiplicativo de E ) 


Pmhados | (00.10. | o5. 102 | 93. 10- | 0.10. 1.10- .10- 
multiplicada 3 3 3 3 3 3 
por 10 . 2.000 0,5 1 0 3,33 0 
3 
Que resulta em: 
Nova linha dois: | 2.000 0,5 1 0 3,33 0 


A nova linha dois obtida irá substituir a linha anterior. Iremos utilizar a nova linha dois para 


zerar os demais elementos não nulos da coluna referente à variável x,. 
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Iniciaremos pela linha zero. Com o intuito de zerar o elemento a, = —15, multiplicaremos a 
nova linha dois por 15. Em seguida, somaremos os resultados obtidos com os respectivos ele- 
mentos da linha zero: 


Linha zero: 0 -12 -15 0 0 0 
Nova linha dois 
Rn 2.000 - 15 = 0,5-15= B H 3,33-15= CS 
multiplicada 30.000 75 1.15=15 0-15=0 50 0:15=0 
por 15: 
Nova linha ze- 
ro (soma da li- 
nha zero com a | 0 + 30.000 = 12+7,5= 15 +15 = E " o 
nova linha dois | 30.000 45 0 PED votados: | DD 
multiplicada 
por 15): 
Portanto, a nova linha zero é: 
30.000 —4,5 0 0 50 0 


Para zerar o elemento a,, = 0,4, iremos multiplicar a nova linha dois por -0,4 e, em seguida, 


somaremos os resultados obtidos com os respectivos elementos da linha um: 


Linha um: 1.000 0,7 0,4 T 0 0 
1 doi PA " ; E 
Da con ostoa- Lews | oca | CO | oco 
ultiplicada por (-0,4) = P -0A o (-0,4) = 0 
-0,4: —800 —1,33 


Nova linha um (so- 


madalinhaum | , 9. 800 04-04= 0 Mies = 
com a nova linha -200 0,7-0,2=0,5 n 140-21 133 0-020 
dois multiplicada ' 


por -0,4): 


Portanto, a nova linha um é: 
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200 0,5 0 1 -1,33 0 


Como o elemento a,, é igual a 0, nào é necessário alterar essa linha. 


Assim, a primeira iteração se resume na tabela a seguir: 


1 2 3 4 5 

z = 30.000 4,5 0 0 50 0 
X; 200 0,5 0 1 -1,33 0 
x, 2.000 0,5 1 0 3,33 0 
X: 700 1 0 0 0 1 


Aqui podemos ver a primeira solução ótima com z = 30.000. 


Como temos ainda valor negativo na variável não básica (x,), devemos proceder a uma no- 


va iteração. 


Segunda iteração — em busca de nova solução 


Atualizando os elementos com os resultados recém-obtidos, podemos construir a seguinte tabela: 


X X X, X X 


7 > 3 4 5 

z=30000 | -45 0 0 50 0 
à, 200 0,5 0 1 -1,33 0 
x, | 2.000 0,5 1 0 3,33 0 
x, 700 1 0 0 0 1 


Ainda temos um coeficiente negativo na linha zero, como já vimos. Assim, a variável que irá 


entrar na base é x,. Para decidirmos qual variável sairá da base, faremos as seguintes divisões: 
200 : 0,5 e 700:1 


O menor resultado obtido é 400, que é o valor da divisão de 200 por 0,5, correspondente à li- 
nha um. Logo, a variável que irá sair da base é x,. 


z = 30.000 —4,5 


x 200 


[5i 
O 
O1 
[e] © 
I 
a 
Q 
[69] 
© © [9] 


2.000 0,5 1 


© 
w 
Q 
[69] 


x 700 1 0 


© 
[e] 
Hm 


Variáveis básicas: x, x, e x, 


Variáveis não básicas: x, e x, 
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Entra: x, 
Sai: x, 


Pivô: a, 2 0,5 


z = 30.000 
zu 1 ea : ; 
O elemento a, é igual 0,5 = PA Nesse caso, multiplicaremos a linha um por 2 (inverso mul- 
TUNE 1 
tiplicativo de z ): 
Linha um -1,33 -2 = 
multiplica- | 200 - 2 = 400 05.221 0-220 1.2=2 66 0-220 
da por 2: ' 
Assim: 
Nox 400 1 0 2 -2,66 0 
um: 
Zerando o elemento da linha zero: 
Linha zero: 30.000 -45 0 0 50 0 
Nova linha 
um multipli- DO menus SO | o.4509 | 29655 | 94500 
1.800 0 =-12 
cada por 4,5: 
Nova linha ze- 
ro (soma da li- 
30.000 
nha zero dim + 1.800= —4,5+4,5= 0+0=0 049-9 50-12- 040-0 
a nova linha 0 38 
E 31.800 
um multipli- 
cada por 4,5): 


Faremos o mesmo para zerar os coeficientes da coluna referente à variável x, nas demais linhas: 


Nova linh 
Min » E 400-(-05)-| 1-(-05)- | 0-C05- | 2-(05)= |-2,66-(-05) | 0- (-0,5) = 
E -200 -0,5 0 -1 =1,33 0 
cada por -0,5: 
Linha dois 2.000 0,5 1 0 3,33 0 
Nova linha 
dois (soma da 
nova linha um | -200 + 2.000 | -0,5 +0,5 = E | 1,33 + 3,33 = E 
multiplicada | =1.800 0 CESSE || res 4,66 Midi 
por -0,5 com a 
linha dois): 
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Nova linha um 
multiplicada 
por -1: 

Linha trés 700 1 0 0 0 jT 


Nova linha trés 
(soma da nova li- B B o. 
nhaummulüpl- 40*70- 1-0 | 04000 uiis ano 0+1=1 

300 -2 2,66 
cada por -1 com 
a linha trés): 


400-(-D)=| 1:(-D= | 0-CD- | 2-(-1)= |-266-(1D=| 0-(-1)= 
-400 -1 0 -2 2,66 0 


Todo esse processo, como vocé pode ver, resulta na tabela a seguir: 


1 2 3 à : 

z - 31.800 0 0 9 38 0 
x 400 1 0 2 -2,66 0 
x, 1.800 0 1 -1 4,66 0 
x, 300 0 0 -2 2,66 1 


Variáveis básicas: x, x, e x, 
Variáveis nào básicas: x, e x, 
Solução ótima: 

x, = 400 

x, = 1,800 

x, = 300 

z = 31.800 


Como todos os coeficientes da linha zero são positivos ou nulos, chegamos à seguinte solu- 
ção ótima para o problema de PL: 


x, =400 
x, = 1.800 
z = 31.800 


Portanto, podemos dizer que a solução ótima para atender à função objetivo será produzir 
400 caminhonetes e 1.800 esportivos, totalizando R$ 31.800,00 de lucros mensais. 
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Como foi a montagem desse modelo de PL? 


Você acompanhou o processo de resolução do problema anterior por meio das iterações (ta- 
belas), nas quais resolvemos um modelo de PL por meio de um método de solução de siste- 
mas de equações lineares. Ficou evidenciado o seguinte: 


e Formulamos o modelo apresentando o problema: os recursos, a disponibilidade, a si- 
tuação atual do processo, o lucro atual. 


e  Montamos um modelo com variáveis de decisão do problema: 
x, - quantidade a produzir de caminhonetes; 
x, - quantidade a produzir de carros esportivos. 

e Apresentamos matematicamente a função objetivo: lucro. 
z = 12x, + 15x, 


e No que se refere às restrições, apresentamos a relação lógica que há no problema e in- 
troduzimos as variáveis de folga, transformando as inequações em equações. Dessa 
maneira, a estrutura lógica que era: 


emprego dos recursos < disponibilidade 

transformou-se em: 

emprego dos recursos + folga = disponibilidade 

relação que traduz como condição o seguinte raciocínio: 

se emprego do recurso < disponibilidade, então a folga > 0 
se emprego do recurso = disponibilidade, então a folga = 0 


Nessa montagem ou modelagem, a variável de folga pode ser expressa por uma variável cuja 
forma seja igual à produção de cada produto. 


Portanto, para resolver o problema introduzindo as variáveis de folga, chegamos ao seguin- 
te modelo: 


z=12x + 15x, 

Em que: substituindo x, e x, por 400 e 1.800 respectivamente (x, = 400, x, = 1.800), teremos: 
z = 12 . 400 + 15 . 1.800 

z = 4.800 + 27.000 


z = 31.800 
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|^ | Exercício resolvido 


Suponha que você tenha uma marcenaria que produz dois tipos de produtos: armário e mesa. 


Para produzir um armário, você usa 3 m? de madeira, 2 L de cola e 50 parafusos. O valor de 
venda lhe rende R$ 300,00 por unidade. 


Para produzir uma mesa, você usa 2 m? de madeira, 1 L de cola e 10 parafusos. O valor de 
venda lhe rende R$ 200,00 por unidade. 


Você tem no estoque da empresa 15 m? de madeira, 10 L de cola e 70 parafusos para traba- 
lhar. Maximize seu lucro. 


Resolução 


Primeiro devemos montar o modelo de PL na forma-padrão, lembrando que o objetivo é ma- 


ximizar o lucro. 


Portanto, identificamos: 


“4. |X, = quantidade de armários a ser produzida 
e as variáveis: 
x, = quantidade de mesas a ser produzida 
e a função objetivo: max z = 300x, + 200x, 
3x, + 2x, < 15 (quantidade de madeira em m?) 
e as restrições: (3x +x, € 10 (quantidade de cola em L) 


50x, + 10x, € 70 (quantidade de parafusos) 


Tendo esses dados, montamos o problema na forma-padrão: 
max z = 300x, + 200x, 
s.a. 3x, +x, <15 


3x, EX <10 
50x, + 10x, <70 


x, X, 2 0 
Tabela 1 
y 
xi X, X4 Xx, Xs 
z=0 —300 —200 0 0 0 
x, | 15 3 2 0 0 
X; 10 3 1 0 1 0 
x, | 70 10 0 0 1 > 
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Variáveis básicas: x, x, e x, 
Variáveis nào básicas: x, e x, 
Entra: x, 

Sai: X, 


Pivô: a, - 50 


z-0 
Tabela 2 
y 
*; x: X, xX, x, 
z-420 | 0 | -140 0 0 6 
x | 10 0 14 1 0 | -0,06 
x, | 58 0 0,4 0 1 | -006 
x | 14 1 0 0 0,02 


Variáveis básicas: x, x, e x, 
Variáveis nào básicas: x, e x, 
Entra: x, 

Sai: x, 

Pivô: a4, = 0,2 


z = 420 


Tabela 3 (final) 


X; X; Xs x: XE 

z-1400 | 700 0 0 0 20 
x, 1 -7 0 1 0 -0,2 
x, 3 -2 0 0 1 -0,1 
x, 7 5 1 0 0 0,1 


Variáveis básicas: x, x, e x, 
Variáveis nào básicas: x, e x, 
Solução ótima: 


x, =0 (VB) 
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x, = 1 
x,=3 
x,=0 (VB) 
z = 1.400 


Resultado: x, = 0; x, = 7; z = 1.400. 
Essa solução foi encontrada pelo método simplex. Portanto, 
e a produção de mesas deverá ser de 7 unidades; e 


e a produção de armários deverá ser de 0 unidade. 


ilização de variáveis artificiais na busca pela solução em PL 


Você irá perceber, diante de determinadas situações, que a utilização do método simplex não 
oferece a solução factível. Quando isso ocorre? Sempre que o problema: 


e apresentar restrições do tipo “=” ou “>”; 


e não apresentar solução inicial factível. 


Nessas condições, você pode aplicar a técnica das variáveis artificiais. Mas o que é a técnica 
das variáveis artificiais? 


Podemos dizer que essa técnica representa a elaboração de um problema auxiliar. Neste, in- 
serimos uma nova variável — a variável artificial —, a qual irá ocupar todos os espaços de res- 
trições nos quais não houve a possibilidade de adicionarmos uma variável de folga (como fa- 
zemos no método simplex). Nessas circunstâncias, a variável artificial passa a ser a variável 


básica para a nova equação e torna possível identificarmos uma solução inicial. 


Lembre-se de que fazemos uso de variáveis artificiais, as quais são um simples artifício ma- 
temático, justamente para suprir a ausência de uma solução básica inicial em um problema 


de PL, condição esta que inviabiliza a aplicação do método simplex na sua forma prototípica. 


Falando em linguagem matemática, essa é a solução para a resolução de uma equação quando 
estamos diante de um problema no qual a matriz A das restrições não apresenta uma subma- 
triz identidade. 


Qual é a forma de operação que essas variáveis realizam? Qual é a finalidade de tal técnica? 
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A finalidade, como já foi visto, é obtermos uma solução básica inicial. A forma como isso ocorre 
no processo de resolução da equação auxiliar, na qual está inserida a variável artificial, é por 
meio da interação com o método simplex, forçando a anulação dessas variáveis uma a uma. 


Existem dois métodos para implementar essa técnica: o Big-M e o método das duas fases. O 
método que você pode utilizar na busca pela maximização ou minimização expressa na fun- 
ção objetivo é o Big-M. 


Método Big-M 


Esse método se baseia no acréscimo de uma variável artificial, em um processo semelhante à 

“folga”, visto no método simplex. A diferença é que nesse cálculo iremos forçar essa variável 
a se anular no fim. O “M” representa um número arbitrariamente grande, que poderemos re- 
presentar por 10.000, 100.000 ou mais, por isso Big-M, ou seja, Grande M, que é o coeficiente 
de penalização atribuído às variáveis. 


Conforme Castillo ([2001?]), nesse método, você pode observar que “as variáveis artificiais são 
“fortemente” penalizadas na função objetivo do problema de PL de modo a provocar 'rapi- 
damente' o seu anulamento. Assim como coeficientes das variáveis artificias na f.o. é introdu- 
zido um parâmetro M (uma constante positiva arbitrariamente grande)” [grifos do original]. 


Utilizamos: 
e -Mnos problemas cuja função objetivo é a maximização de uma condição; 
e M nos problemas cuja função objetivo é a minimização de uma condição. 


Nessa operação, quando você for selecionar os vetores para entrar na base da equação, é im- 
portante escolher apenas entre aqueles que não são artificiais, pois a penalidade de “M” apli- 
cada às variáveis artificiais impossibilita a reentrada destas. 


Condições de resolução do método Big-M 
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e Você pode considerar a solução básica admissível do problema auxiliar equivalente à 


do problema original quando conseguir anular as variáveis artificiais. 


e Quando a solução ótima que você encontrar para o problema auxiliar corresponder a 


uma solução básica admissível do problema original, pode considerá-la ótima para o 


problema original. 


e Sempre que você for determinar a variável que entra, deve escolher a correspondente 
ao maior custo reduzido, entre as de custo reduzido positivo, da mesma forma que 


vocé faz no método simplex. 


e No entanto, quando você for determinar a variável que sai, deve escolher a que alcan- 
çar o menor quociente. 


Agora, vamos transformar isso em linguagem matemática expressa na forma de equações. 


Vejamos o exemplo. 


Suponha que, no exemplo do Capítulo 2 (fábrica de miniaturas de automóveis), exista ainda 
a restrição de que pelo menos 500 caminhonetes devem ser produzidas para atender a uma 
demanda mínima. Nesse caso, temos de acrescentar a restrição x, > 500 para indicar que a va- 
riável x, deve ser maior ou igual a 500. Nesse caso, a formulação do problema é: 


max z = 12x + 15x, 


sa. — 07x, + 04x, < 1.000 
0,15x, 03x, < 600 


x « 700 
x 2500 
x,20,x,20 


Note que a função objetivo e as restrições são as mesmas. Apenas acrescentamos uma nova restrição. 


Para resolver o problema utilizando o método simplex, precisamos adicionar as variáveis de folga: 


maxz = 12x, + 15x, 


s.a. 0,7x, + 0,4x, + x, 
0,15x,+0,3x, +x, = 1.000 
X, tX. = 600 
X, =X -700 
x, > 0, x 6 =500 


Mas note que, como a restrição x, > 500 é de valor maior ou igual, para obtermos a igualdade, a va- 
riável de folga x, tem sinal negativo. Por esse motivo, não conseguimos formar uma solução inicial. 


Para que possamos ter uma solução básica inicial, é necessário, portanto, acrescentar mais 
uma variável à última restrição, mas com coeficiente positivo. Essa variável é denominada 
variável artificial. 


Como essa nova variável não existe de fato no problema original, precisamos fazer com 
que ela saia da base. Para isso, é necessário acrescentar uma penalidade muito alta na função 
objetivo para forçar a saída da variável artificial da base. Em outras palavras, teremos um 
custo extremamente alto para manter a variável artificial na base. Esse custo alto é represen- 
tado pela letra M. Logo, o problema pode ser escrito como: 


maxz = 12x, + 15x,+ 0x, + Ox, + 0x,+ 0x, — Mx, 


s.a. 0,7x, + 0,4x, + x, = 1.000 
0,15x, + 0,3x, FX = 600 
X, +X, = 700 
X =X, tX, = 500 
x,20,x,20 
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Escrevendo o problema sob a forma de tabela, temos a tabela inicial: 


1 2 3 n 5 6 7 

z-0 -12 -15 0 0 0 0 M 
X. 1.000 0,7 0,4 1 0 0 0 0 
X, 600 0,15 0,3 0 1 0 0 0 
x: 700 1 0 0 0 1 0 0 
X. 500 i 0 0 0 0 -1 


Assim como no método simplex, trocamos os sinais dos coeficientes da função objetivo para co- 


locá-los na tabela. As variáveis básicas são: x, x, X„ e X,. As variáveis não básicas são: x, X, € x.. 


Como x, é uma variável básica, precisamos que seu coeficiente na função objetivo seja igual 
a zero. Para que isso aconteça, basta multiplicar a linha 4 por -M e somar com a linha zero: 


Linha zero: 0 -12 -15 0 0 0 0 M 
Linha quatro multiplica- | 500- (-M) = | 1- (-M) = 0- (MO 0- (-M) =| 0- (M) =| 0- (M) = -1- (M) =| 1- (-M)= 
da por -M: -500M -M o 0 0 0 M -M 
Nova linha zero (soma 
da linha zero com a li- | 0+(-500M) =| -12-M= | —15+0= 0+M= |-M*M- 
NA 0-020, 0+0=0 | 0+0=0 
nha quatro multiplicada -500M -12-M -15 M 0 
por -M): 


Substituiremos, então, a linha zero pela nova linha zero. 


Observe que, como a variável x, está na base, o valor de z é -500M, um prejuízo gigantesco. 
Para decidirmos quem entrará na base, basta considerar, na função objetivo, o coeficiente 
mais negativo. Nesse caso, x, entra na base. Dividindo cada termo independente pelos res- 
pectivos coeficientes positivos (valores negativos ou nulos não são considerados) da coluna 
referente à variável X; temos: 


1.000 : 0,7 = 1.428,57 
600 : 0,15 = 4.000 
700 : 1-2 700 
500:1- 500 


O menor resultado obtido foi o da divisão de 500 por 1. Logo, a variável que sairá da base é x,. 
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xi x, Xz X, Xs Xe x. 

15 | 0 0 0 M 0 

X, 04 | 1 0 0 0 0 

x, 03 | 0 1 0 0 0 

x, 0 0 0 1 0 0 
x, 0 0 0 o | a 1ı |> 


O pivô, nesse caso, é igual a 1. Como vimos anteriormente, precisamos agora zerar os demais 


elementos da coluna pivô. Seguindo os passos do método simplex, construímos a seguinte tabela: 


z = 6.000 


x 650 = 


x 525 


o 
[es 
[e] 
Hm 


[5] 
[se] ce ce c 
o 
NS 
ra 
[=] © 
© [=] © 
© 
MI 
| 
O 
MI 


x 200 


500 1 0 


© 
© 
© 


-1 1 


A variável que entra na base é x, e a variável que sai é x,. Realizando as operações para trans- 


formar o pivô em 1 e para zerar os demais elementos da coluna pivô, chegamos à tabela fi- 
nal, ou seja: 


x x, X4 x, X. X6 x, 
z-3037 | 0 0 |375| 0 0 | 142 |-14254M 
x, 1.625 0 " v NE 0 | 17 | A75 
x, 37,5 0 o k 1 0 [-0375 | 0,375 
x, 200 0 0 0 0 1 1 2 
x, 500 1 0 0 0 0 E 1 


Como todos os coeficientes da função objetivo são números não negativos, obtivemos a so- 


lução ótima, que é: 


x, = 500 
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x, = 1.625 
x, 7 37,5 
x, = 200 


z = 30.375 
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Ou seja, para que o lucro da empresa seja o maior possível, é preciso fabricar 500 caminho- 
netes e 1.625 carros esportivos. Para essas quantidades, o lucro máximo é de R$ 30.375,00. 


Note que x, indica a sobra de alumínio e x, a quantidade de caminhonetes que deixou de 
ser produzida em relação ao máximo que poderia ser vendido. 


Para saber mais 


GOMIDE, F. Algoritmos numéricos de busca em otimização. Conteúdo para a aula de 
Planejamento e Análise de Sistemas de Produção da Faculdade de Engenharia Elétrica e 
de Computação da Unicamp. Disponível em: <http://w ww .dca.fee.unicamp.br/-gomide/ 


courses/EA044/transp/EA 044 BuscaNumericaOtimizacao.pdf>. Acesso em: 20 nov. 2009. 


Nesse material, você encontra um excelente resumo sobre algoritmos de busca, aplicações, 


problemas de maximização e minimização, bem como sobre o método BIG-M. 


Coeficiente de penalização atribuído às variáveis M 


Na aplicação do método Big-M para o algoritmo simplex, você pode chegar às seguintes 
situações: 


e eliminação de todos os vetores artificiais da base; logo, você obteve uma solução 


admissível; 


e impossibilidade de eliminar todos os vetores artificiais da base, sendo que os custos 
reduzidos são não positivos, o que significa que você encontrou um problema impos- 
sível de solucionar ou, então, encontrou uma solução admissível por meio da substi- 
tuição de uma solução inicial degenerada ou da eliminação de restrições redundantes. 


“1 


= 


— 


/ | 


l | 


L | Exercicio resolvido 


Suponhamos - e aqui vamos utilizar um exemplo apresentado pela professora Gladys Castillo 
([2001?]) — que temos a seguinte situação como um problema de PL na forma-padrão: 


max Z = 3x, *5x, 


sa. X <4 
2x, <12 
3x, + 2x, =18 


x,x,20 
Como vocé pode perceber: 


e na terceira restrição há um sinal de igualdade; 


e para reduzir o problema na forma-padrão, devemos adicionar duas variáveis de folga, 


isto é, x, e x, pois a restrição do problema é uma restrição de igualdade. 
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maxz- 3x, *5x, 


S.a. X EX ES 


> 
X, X,, X, X, 20 


Como resultado dessa interferéncia, criamos um problema: a matriz A náo contém uma ma- 
triz identidade. 


Vamos, portanto, verificar o que detectamos nessa fase de resolução do problema: 


* avariável de folga para utilizarmos como variável básica inicial para a equação 3 nào 
foi encontrada; 


* amatriz A, representativa do sistema de equações, não apresenta uma submatriz iden- 
tidade B. 


Ao detectarmos tais condições no processo de resolução do problema, torna-se difícil a iden- 
tificação de uma SBA (solução básica apresentada) inicial. 


Como já dissemos, o procedimento que usualmente você pode adotar nesse caso é aplicar a técnica 
das variáveis artificiais, pois esta se constitui em “um procedimento integrado no método sim- 
plex que permite ultrapassar o desconhecimento de qualquer SBA inicial num problema de PL na 
forma-padrão” (Castillo, [2001?]). 


Portanto, em todas as restrições nas quais não conseguimos adicionar uma variável de folga 
(como ocorre no simplex), a variável artificial passa a ser a variável básica! para a equação 
em foco. 


O que fazer? Podemos dizer que essa técnica se traduz pela elaboração de um problema au- 
xiliar, por meio do qual é inserida uma nova variável — a variável artificial. 


Lembre-se de que não devemos confundir as variáveis artificiais com as variáveis de folga, pois 
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aquelas não apresentam nenhum valor econômico, correspondem a mero artifício matemático. 


Como esse problema não contém uma submatriz identidade, adicionamos uma variável ar- 


i 

tificial Xs E 

E 

maxz  -3x,*5x,- Mx, 2 
s.a. x, +X, =4 
2x, *x, =12 
11 2x, +x, -18 


X, X,, X, X, X, Z0 


1 Uma variável pode ser caracterizada como básica uma vez que apresente coeficiente 1 na equação apresentada e coeficientes nulos 
nas demais. 
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Agora é possível identificar uma SBA inicial: X, - (0, 0, 4, 12, 18). 


Lembramos que SBA inicial significa a primeira solução básica apresentada pelo método simplex. 


Na busca por uma solução ótima nos processos de planejamento das atividades de uma or- 
ganização com a utilização da PO, vimos as características específicas de utilização dos mé- 
todos simplex e do Big-M. Por meio da aplicação das operações de maximização, destacamos 
os passos de resolução de situações-problema com a aplicação de ambos os métodos para a 
elaboração de modelos de resolução. Assim, construímos exemplos práticos para subsidiar 
as atividades relativas à identificação, estruturação e resolução matemática de problemas que 
influenciam nas decisões do cotidiano de uma organização. 


1. Uma pequena empresa de artefatos de couro apresenta dois segmentos de produtos: car- 
teiras e bolsas em couro. Cada carteira representa um lucro líquido de R$ 28,00 e cada bol- 
sa, de R$ 92,00. Para produzir uma carteira, a empresa utiliza 0,25 m? de couro e, para pro- 
duzir uma bolsa, são necessários 2,4 m? de couro. A empresa tem um estoque disponível de 
16.500 m? de couro cru. Sabemos ainda que a empresa precisa produzir pelo menos 100 bol- 
sas e 150 carteiras para atender a alguns pedidos. Qual quantia a empresa deve produzir de 
cada produto para que o lucro seja o maior possível? 


2. Uma marcenaria possui 30 peças de madeira e dispõe de 44 horas de trabalho por semana 
para confeccionar biombos ornamentais. Dois modelos venderam muito bem no ano passado, 
de maneira que ela se limitou a fabricar esses dois tipos. O modelo I requer 2 peças de ma- 
deira e 1 hora de trabalho, enquanto o modelo II necessita de 1 peça de madeira e 2 horas de 
trabalho. Os preços dos modelos são, respectivamente, R$ 120,00 e R$ 80,00. Para atender à 
demanda, a marcenaria deverá produzir a quantidade mínima de 5 biombos do modelo I e 7 
biombos do modelo II. Quantos biombos de cada modelo a empresa deve montar se desejar 
maximizar o rendimento obtido com as vendas? 


3. Um fabricante produz dois tipos de aço: normal e especial. Uma tonelada de aço normal re- 
quer 2 horas no forno de soleira aberta e 5 horas de molho; 1 tonelada de aço especial requer 
2 horas no forno de soleira aberta e 3 horas de molho. O forno de soleira aberta está disponí- 
vel 8 horas por dia, e o molho está disponível 15 horas por dia. O lucro em 1 tonelada de aço 
normal é de $ 120,00 e de $ 100,00 para 1 tonelada de aço especial. A empresa precisa produzir 
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diariamente no mínimo 2 toneladas de aco normal e 1 tonelada de aco especial. Determine 


quantas toneladas diárias de cada tipo de aço devem ser produzidas para maximizar o lucro. 


4. Já apresentamos os passos de resolução de situações-problema com a aplicação do méto- 
do simplex e do método Big-M para a elaboração de modelos de resolução. Nesse processo, 
quais foram as mudanças e operações aplicadas? 


5. Quais são as principais características do método Big-M? 
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capítulo 


Utilização do WinQSB na resolução 
de problemas de programação linear 
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Conteüdos do capítulo 


e Como utilizar o WinOSB na resolução de problemas de programação linear. 


Após o estudo deste capítulo, vocé será capaz de: 


1. resolver problemas de programação linear utilizando o WinQSB. 


O WinOSB, um software gratuito desenvolvido por Yih-Long Chang, é um pacote de ferramentas que 
objetiva servir de suporte para a tomada de decisões baseada em problemas de pesquisa operacio- 
nal (PO). Sua primeira versão, intitulada Quantitative System for Business (QSB), era bastante abran- 
gente, mas com uma interface simples, baseada no ambiente DOS. A versão atual conta com uma 


interface mais amigável, incluindo recursos gráficos mais avançados do que os da primeira versão. 


4.1 Utilizando o WinQSB na resolução de um problema de PL 


Usaremos o exemplo da fábrica de miniaturas (Capítulo 2) para ilustrar as principais funções 
do WinQSB utilizadas na resolução de problemas de programação linear (PL). Relembrando, 


as informações são as seguintes: 
e A empresa fabrica dois modelos de veículos em miniatura: caminhonetes e esportivos. 


e As principais matérias-primas empregadas para a fabricação dos brinquedos são plás- 
tico e alumínio. 


e O modelo esportivo consome 400 g de plástico e 300 g de alumínio. 

* Acaminhonete requer 700 g de plástico e 150 g de alumínio. 

* A disponibilidade mensal de plástico é de 1 tonelada e a de alumínio, de 600 kg. 
e Olucro unitário referente à caminhonete é de R$ 12,00. 

e O lucro unitário referente ao modelo esportivo é de R$ 15,00. 

e A produção de caminhonetes não pode ultrapassar 700 unidades. 


Você viu também que a formulação do problema era a seguinte: 
max z = 12x, + 15x, 


sa. 0,7x,+0,4x, «1.000 


0,15x, + 0,3x, < 600 
x * «x 700 
x, 20,x, 20 


Você vai perceber que é muito simples começar a utilizar o WinOSB. Após a instalação, basta 
clicar em Iniciar, Programas, WinQSB. Você verá um menu com várias opções. Para resolver 
os problemas de PL, utilize a opção Linear and Integer Programming. 
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Figura 4.1 - Menu do WinOSB 


Acceptance Sampling Analysis 
Aggregate Planning 

Decision Analysis 

Dynamic Programming 

Facility Location and Layout 


Forecasting and Linear Regression 
Goal Programming 

Inventory Theory and System 

Job Scheduling 

Linear and Integer Programming 
MarKov Process 

Material Requirements Planning 
Network Modeling 
Nonlinear Programming 
PERT CPM 

Quadratic Programming 
Quality Control Chart 
Queuing Analysis 

Queuing System Simulation 


A EI PERELE Ein 


Ao clicar nessa opção, aparecerá a tela representada na figura a seguir. 


Figura 4.2 - Tela de abertura da opção Linear and Integer Progamming 


a Linear and Integer Programming 


Linear and Integer 
Programming 


Version 2.00 
Copyright O Yih-Long Chang 


Licensed to. User 
Company User 


Pltoaso wait while loading program 


Em seguida, clique em File e escolha a opção New Problem. 
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Figura 4.3 - Menu File, opcào New Problem 


a 


$8 Linear and Integer Programming 


Agora, você conhecerá um dos passos mais importantes para a resolução de um problema de 
PL por meio do WinOSB. Escolha um nome para o problema e escreva-o no campo denomi- 
nado Problem Title. O número de variáveis deve ser colocado no campo Number of Variables, 
e o número de restrições, em Number of Constraints. O critério da função objetivo (maximi- 


zar ou minimizar) deve ser indicado no item Objective Criterion, selecionando-se a opção 


Maximization ou Minimization, de acordo com o problema. A forma de entrada dos dados de- 
ve ser indicada em Data Entry Format. O ideal é deixar selecionada a forma matricial, indica- 
da por Spreadsheet Matrix Form. Em seguida, é necessário indicar o tipo de variável do pro- 


blema. A primeira opção, chamada Nonnegative continuous, indica variáveis contínuas e não 
negativas. A segunda opção é Nonnegative integer. Utilize essa opção quando o problema se 


referir a variáveis inteiras. A terceira opção, que indica variáveis binárias — Binary (0,1) —, de- 
ve ser marcada em problemas cujas variáveis podem assumir exclusivamente um dos se- 
guintes valores: 0 ou 1. Finalmente, a última opção — Unsigned/unrestricted — indica variáveis 
irrestritas, utilizadas em problemas nos quais as variáveis podem assumir valores positivos, 
negativos ou zero. 


Figura 4.4 — Tela das especificações do problema 


LP-ILP Problem Specification 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.810/1998 e punido pelo art. 184 do Código Penal. 


No caso do problema de miniaturas, temos duas variáveis inteiras, pois nào faz sentido a pro- 
dução de miniaturas fracionadas - o número de restrições do problema é igual a trés. Como 
o problema é de maximização, marque essa opção. Logo, o preenchimento ficará como re- 
presentado na figura a seguir. 


Figura 4.5 — Tela das especificações do problema com os campos preenchidos 


LP-ILP Problem Specification 


Feito isso, é só clicar em OK. Você terá, então, uma tabela na qual poderá inserir os dados do 
problema: coeficientes da função objetivo e das restrições e termos independentes!. 


Figura 4.6 — Tabela para a inserção dos dados do problema 


Preenchendo-se os respectivos campos, a tabela ficará conforme a figura a seguir. 


Figura 4.7 — Tabela preenchida com os dados do problema 


1 Right Hand Side (R.H.S.), que representa o lado direito das inequações. 
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Note que o separador de casas decimais utilizado no WinQSB é o ponto, e não a vírgula, que 
é o padrão ao qual estamos acostumados. 


Após preencher os dados do problema, é necessário que você obtenha a solução. Para isso, 
clique na opção Solve and Analyze e, em seguida, em Solve the Problem. 


Figura 4.8 - Menu Solve and Analyze, opção Solve the Problem 


Utilities Window 


4 Solve the Problem O 


F Solve and Display Steps 
| Graphic Method 


Change Integer Tolerance 
Specify Solution Quality 
Specify variable Branching Priorities 


Ao fazer essa opção, aparecerá uma mensagem informando que o problema foi resolvido. 
Clique, entào, em OK. 


Figura 4.9 - Mensagem de alerta 


Linear and Integer Programming ES 


The problem has been solved. 
Optimal solution is achieved. 


A solução do problema é fornecida em uma tabela similar à apresentada na Figura 4.10. 


Figura 4.10 - Solução do problema 


Reduced Basis Allowable 
Cost Status Min. cti) 


0 basic 7.5000 
o basic 6.8571 


Constraint 


c 
cz 


Os valores de x, e x, são, respectivamente, 400 e 1.800, e o valor de z é 31.800 (Objective Function 
(Max.) =). O WinQSB também informa os custos reduzidos e as folgas das variáveis, além do 
preço sombra, elementos abordados no Capítulo 7. 
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E possível não só escolher nomes para as variáveis e as restrições, como também adicionar 
ou excluir variáveis e restrições, além de alterar o critério da função objetivo. Essas opções 
podem ser encontradas no menu Edit. 


Figura 4.11 — Menu Edit 


Format Solve and Analyze Results 


Cut Ctrl+x 

Copy Ctrl+C 
| Paste Ctri-V 
é Clear 

Problem Name 


Variable Names 
Constraint Names 
Objective Function Criterion 


Insert a Variable 
Delete a Variable 


Insert a Constraint 
Delete a Constraint 


O tamanho da fonte, o alinhamento dos dados e outras opções podem ser alterados no me- 
nu Format. 


Figura 4.12 — Menu Format 


Solve and Analyze Results 


Number 
Font 
Alignment 


Row Height 
Column Width 


Switch to Normal Model Form 
Switch to Dual Form 


Também é possível obter a solução gráfica para problemas de duas variáveis. Nesse caso, 
escolha a opção Graphic Method no menu Solve and Analyze ou clique no ícone : 
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Figura 4.13 - Menu Solve and Analyze, opção Graphic Method 


Results Utilities — wWinc 
Solve the Problem 
Solve and Display Steps 
Graphic Method 


Perform Parametric Analysis 
Alternative Solution 


Change Integer Tolerance 
Specify Solution Quality 
Specify variable Branching Priorities 


Escolhendo essa opção, abrirá uma janela com as opções para a variável a ser plotada no ei- 
xo horizontal e para a variável referente ao eixo vertical. Clicando-se em OK, o WinQSB fará 
o gráfico. 


Figura 4.14 - Seleção de variáveis do método gráfico 
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Select Variables for Graphic Method 


X (horizontal) axis Y (vertical) axis 


x1 


ES ee 
x2 


No gráfico, você verá cada restrição do problema e também a função objetivo, além da solu- 
ção ótima desejada. 


Figura 4.15 — Solução gráfica 


x2 Constraint 
2.500,00 


Objetive Function -------- Feasible Area 


Optimal Solution 


2.250,00 OBJ = 31.800,00 


X1= 400,00 


2.000,00 X2 = 1.800,00 


1.750,00 


1.500,00 


1.250,00 


1.000,00 


750,00 


500,00 


250,00 


0,00 


0,00 900,00 1.600,00 2.400,00 3.200,00 4.000,00 


x1 


A utilização de ferramentas computacionais facilita, e muito, a obtenção de respostas para 
problemas de otimização. Neste capítulo, vimos como utilizar o software WinQSB na resolu- 
ção de problemas de PL. Muitas empresas se valem de softwares comerciais com esse intui- 
to; em outros casos, em decorrência de diversas particularidades, é necessário desenvolver 
softwares específicos para determinadas empresas. 


1. Uma marcenaria fabrica estantes e balcões e deseja otimizar a sua produção, com o objetivo 
de maximizar seu lucro mensal. Atualmente, a empresa possui 1.000 m? de madeira e 500 m? 
de vidro, as principais matérias-primas utilizadas na produção desses artigos. Para fabricar 
uma estante, são necessários 10 m? de madeira e 4 m? de vidro. Cada balcão requer 7 m? de 
madeira e 1 m? de vidro. A cada mês, a capacidade máxima de produção está limitada a 80 es- 
tantes e 100 balcões. Sabe-se também que o lucro para cada estante é de R$ 750,00, enquanto o 
lucro para cada balcão vendido é de R$ 640,00. Formule a questão como um problema de PL. 


2. Uma indústria alimentícia produz massa para macarrão e para pastéis, usando como prin- 
cipais ingredientes farinha e ovos. Para produzir um pacote de massa para macarrão, são uti- 
lizados 1 kg de farinha e 2 ovos; para um pacote de massa para pastéis, são necessários 1/2 


73 


£ 
S, 
E) 
k 
3 
$ 
E 
8 
E: 
É 
5 
e 
8 
2 
E: 


Lei n? 9.610/1994 


uer meio ou forma sem a pi 


io poderá ser reproduzida p: 


74 


kg de farinha e 3 ovos. As quantidades de farinha e de ovos disponíveis por mês são, respec- 
tivamente, 120 mil kg e 250 mil unidades. Sabe-se também que, em razão das limitações das 
máquinas, a indústria pode produzir no máximo 110 mil pacotes de massa para macarrão e 
170 mil pacotes de massa para pastéis. O lucro líquido referente a cada pacote de massa para 
macarrão é de R$ 0,70, e o lucro líquido referente a cada pacote de massa para pastéis é de R$ 
0,45. Com o objetivo de maximizar o lucro dessa indústria, determine qual deve ser a produ- 
ção mensal de cada um dos dois tipos de massa. 


3. Uma empresa fabrica dois artigos de camping: sacos de dormir e barracas. Cada saco de dor- 
mir requer 2 horas para que os tecidos sejam cortados, 5 horas para costurá-los e 1 hora para 
impermeabilizá-los. Cada barraca requer 1 hora para que sejam cortados os tecidos, 5 horas 
para costurá-los e 3 horas para impermeabilizá-los. Dados os recursos limitados da empresa, 
ela dispõe de 14 horas para o corte, 40 horas para a costura e 18 horas para a impermeabiliza- 
ção. A margem de lucro é de R$ 50,00 por saco de dormir e de R$ 30,00 por barraca. Formule 
e resolva o problema. 


4. Uma emissora de rádio está fazendo algumas alterações na programação com o objetivo 
de aumentar a audiência. Para isso, pretende reorganizar o número de apresentações sema- 
nais de cada programa. O programa A, dedicado à música pop, tem 50 minutos de música e 
10 minutos de comerciais, atingindo uma audiência de aproximadamente 22 mil ouvintes. O 
programa B, dedicado à música clássica, tem 45 minutos de música e 15 minutos de comer- 
ciais, alcançando aproximadamente 15 mil ouvintes. A emissora, com esses programas, pre- 
tende ter pelo menos 105 minutos semanais de comerciais e, no máximo, 435 minutos de mú- 
sica. Determine quantas vezes cada programa será apresentado semanalmente de modo que 
a respectiva audiência semanal seja maximizada. 


5. O hipoclorito de sódio, produto obtido a partir de uma reação de cloro com hidróxido de só- 
dio, é uma substância muito utilizada como alvejante ou desinfetante. Uma indústria química 

fabrica três produtos para diferentes finalidades: alvejante de uso doméstico, com 2,5% de hipo- 
clorito de sódio; clorador de uso geral, com 12% de hipoclorito de sódio; e clorador para pisci- 
nas, com 30% de hipoclorito de sódio. Cada litro de hipoclorito de sódio custa R$ 3,00. O preço 

de venda do litro de alvejante é R$ 2,30. Os custos unitários totalizam R$ 1,00, referentes à em- 
balagem, ao rótulo e a demais componentes, exceto o hipoclorito de sódio. De maneira similar, 
cada litro de clorador de uso geral custa R$ 0,90, mais o custo referente ao hipoclorito de sódio, 
e é vendido por R$ 3,90. Cada litro de clorador para piscina custa R$ 0,85, mais o custo referente 

ao hipoclorito de sódio, e é vendido por R$ 4,70. A disponibilidade semanal de hipoclorito de 

sódio corresponde a 700 mil litros. Por questões contratuais, a produção mínima de alvejante 

deve ser de 2 milhões de litros, e a produção de clorador para piscinas nào pode ultrapassar 40% 

da produção de clorador de uso geral. A indústria tem à disposição 7 milhões de embalagens 

por semana. O objetivo é determinar qual deve ser a produção semanal que maximize o lucro. 
Formule e resolva o problema tendo como base o modelo de PL. 
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capítulo 


A pesquisa operacional (PO) 
e os problemas de otimização em redes 
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Conteüdos do capítulo 


e Utilização da pesquisa operacional para resolver problemas de transporte. 
e Aplicação dos métodos do custo mínimo, de Vogel e do canto noroeste na solução 
de problemas de transporte. 


e Resolução de problemas que envolvem árvores mínimas. 


Após o estudo deste capítulo, você será capaz de: 


1. identificar as variáveis básicas de problemas de transporte; 

2. elaborar um modelo matemático para a resolução de um problema 
de transporte; 

3. aplicar os métodos do custo mínimo de Vogel e do canto noroeste na otimização 
de um problema de transporte; 


4. resolver problemas que envolvem árvores mínimas. 


Problemas que envolvem transportes são interessantes temas da pesquisa operacional (PO), 
pois apresentam diversas variáveis de origem (fábricas, depósitos, lojas etc.) e de destino (fi- 
liais, postos de distribuição, obras etc.). 


Esse emaranhado de rotas (ir e vir) é justamente o que nos permite formular problemas de 
transportes, os quais podem ser resolvidos com problemas de programação linear (PL). Sobre 
essa situação, Passos (2008, p. 182) afirma que “esses problemas buscam encontrar rotas pa- 
ra a transferência de cargas de diferentes origens para diferentes destinos minimizando os 


seus custos”. 


Nesse contexto, sabemos que o custo de transporte impacta diretamente os preços dos pro- 
dutos para o cliente final. São os chamados custos logísticos, como tarifas das distâncias e vo- 
lume de carga que se quer transportar. Essas variáveis interferem diretamente na logística, 
o que leva o departamento de logística de uma empresa à busca por tirar o máximo proveito 


das vantagens de localização. 


Existem métodos específicos para a resolução desses problemas, bem como softwares exclusi- 
vos que calculam e apresentam melhores cenários. Mostraremos algumas situações que usu- 


almente permitem compreender melhor esses cenários por meio da PL. 


Para a resolução desses problemas, existem vários métodos, como o Vogel e o canto noroes- 
te. No entanto, o método do custo mínimo, na maioria das vezes, possibilita que se obtenha 
a solução mais rapidamente. Mas, para que você amplie seus conhecimentos, apresentaremos 


esses três métodos com exemplos de resolução de problemas. 
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É importante que vocé observe, para se inserir no contexto, que, na maioria das vezes, o que 
pretendemos é minimizar o custo — processo designado como Min C, o qual é compreendi- 
do por origens e destinos. Ainda na elaboração desses cálculos, as disponibilidades ou ofer- 
tas (origens) e as demandas (destinos) podem aparecer sempre com as relações de igualdade, 
menor igual ou maior igual. 


Ou seja, a quantidade existente na oferta (origem) pode ser maior ou igual ou menor que 
a quantidade demandada no destino. 


Quando a quantidade de origem é igual à quantidade de destino, dizemos que os sistemas 
se encontram equilibrados. 


quantidade de origem = quantidade de destino — sistema equilibrado 


Temos sistemas desequilibrados quando as quantidades de origem são diferentes das quan- 
tidades de destino. 


Vejamos o esquema que representa essa desigualdade logo a seguir. 


Figura 5.1 — Sistemas desequilibrados 


> | Destino 1 
Origem 1 
8 x Destino 2 
Origem 2 Pw q 
EB s 
E assi E Destino 3 
Origem 3 
Destino 4 


Quando elaboramos esse tipo de representação geométrica em um problema com dados de 
custos, utilizamos as setas para indicar o custo unitário pela carga a ser transportada para os 
respectivos destinos (como você verá em exemplos, logo a seguir). Você pode observar nessa 
figura que o produto de uma empresa sai (origem) de trés lugares diferentes e vai (destino) 
para quatro lugares distintos; portanto, há uma desigualdade numérica entre as fontes de 
origem (três) e os pontos de chegada ou destino (quatro). Essas são circunstâncias que carac- 


terizam um sistema em desequilíbrio. 
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Problemas sobre quantias a serem transportadas (sistema equilibrado) 


Vamos tornar essa operação prática, aplicando os métodos de resolução em exemplos. Para 
isso, vamos imaginar algumas situações-problema relacionadas a transportes nas quais são 


considerados as quantidades, a origem e o destino. 


Solução pelo método do custo mínimo 


Vamos supor uma situação na qual estejam envolvidos uma indústria de motocicletas e o seu 
processo de distribuição da produção. Nesse contexto, vamos considerar que a nossa neces- 
sidade é saber qual é a quantidade de motocicletas que deverá ser transportada de cada 
fábrica para cada depósito de modo que o custo total de transporte seja o menor possível. 


Vamos aos dados! Digamos que essa indústria de motocicletas possui quatro fábricas distri- 
buídas em diferentes cidades do mundo: Amsterdã, Cairo, Nova Iorque e Cidade do México. 
Cada fábrica apresenta uma capacidade semanal de produção conforme a seguir: 


Fábricas Amsterdã Cairo | Nova Iorque Cidade do México 


Capacidade 70 105 130 65 


As motocicletas são transportadas para quatro depósitos para que, em seguida, sejam distri- 
buídas às lojas. As demandas de cada depósito são as seguintes: 


Depósitos Califórnia Teerã Rio de Janeiro Madri 


Demandas 50 80 150 90 


Outro dado necessário para realizarmos essa operação diz respeito aos custos referentes ao trans- 
porte de uma motocicleta, de cada fábrica para cada depósito. Assim, apresentamos na próxima 
tabela os custos em dólar: 


Depósitos 
Califórnia | Teerã | Rio de Janeiro | Madri 
Amsterdã U$ 8 U$ 8 U$5 U$2 
3 Cairo U$ 10 U$5 U$ 10 U$ 10 
S Nova Iorque U$7 U$ 10 U$6 U$ 8 
Cidade do México U$ 4 U$ 10 U$ 10 U$5 


Note que a formulação do problema de transporte pode ser dada como a formulação de um 
problema de PL. Nesse caso, cada variável (x) apresenta dois Índices (i e j): 


e į- primeiro número (índice) indica a origem (fábrica); 


e j- segundo número (índice) representa o destino (depósito). 
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Por exemplo: 
e x, indica a quantidade a ser transportada da fábrica 1 (Amsterdã) para o depósito 1 (Califórnia), 


e x, indica a quantidade a ser transportada da fábrica 1 (Amsterdã) para o depósito 2 
(Teerã) e assim por diante. 


Ao formularmos o problema por PL, devemos definir como objetivo a minimização do custo to- 


tal do transporte das motocicletas para os depósitos. 


Figura 5.2 — Transporte das motocicletas 


Amsterdã Califórnia 2 
70 50 E 
Cairo Teerã E 
105 80 z 
E 
8 
Rio de 3 
Nova lorque Janeiro E 
130 150 i 
Cidade do N a 2 ] | Madri 
México u M NN 
- 90 
a 84 
65 s]- 


Min z = 8x, + 8X, + 5X; + 2x, + 10x, + 5X, 10x,, + 10x,, + 7X4 + 10x, + 6x. + 8X, + 4x, + 


10x,, + 10x,, + 9X, 
Sujeito a: 
e Restrições de capacidade de produção: 
Xa POLL É 70 
Xat Xat Xat X x 105 


X4 t Xy t X4 * X, < 130 
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X, *X54* X4 * Xy € 65 


e Restrições de demanda de depósito: 


Xu + Xz T Xs i Xa 2 50 
Xp * X» i X» Es Xp 2 80 

Xi Tj Xs T Xs u Xag 2 150 
Xu E Xy E Xu T Xu 2 90 
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e Restrições de nào negatividade: x, 2 0, inteiros 


Em que: x, é a quantidade de suprimentos a serem transportados da fábrica (i) para o depó- 
sito (j). Ou seja, 


e i= fábrica; 

e j= depósito; 

e x= quantidade. 
Lembre-se também de que: 


e minz significa “minimizar” a função objetivo, no caso, os custos do transporte. 


Cálculo pelo método do custo mínimo 


Para facilitar o fluxo do problema apresentado pelo cálculo do método do custo mínimo, uti- 
lizaremos a tabela a seguir (essa tabela tem o mesmo significado matricial, que será visto no 
Capítulo 8). 


Figura 5.3 - Modelo genérico para o cálculo pelo método do custo mínimo 


Destinos Capacidade 
Destino 1 | Destino2 | Destino3 | Destino 4 de oferta 
Origens (disponibilidade) 


Origem 1 


LH LI E] 
LS AS LAT JL 
E LI L 

ll 


Origem 4 


Demanda 
(capacidade) 


Nos quadrados menores, colocamos os valores dos custos de transporte das motocicletas das 
origens (linhas) para os destinos (colunas). Aliás, podemos nomear as linhas e as colunas do 
quadro para facilitar a visualização, como você pode ver a seguir. 
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Destinos Rio de Capacidade 
Califórnia Teerã ; Madri (disponibili- 
Origens Janeiro dade) 


Amsterdã 


Cairo 


Nova Iorque 


Cidade do 
México 


© 
EEB 
ERES 
EE 


UE 
[e] [e] 
JE 
JE 


Demanda 
(necessidade) 


Agora, preenchemos os quadros maiores com as capacidades de cada fábrica e com as deman- 
das dos destinos. Na última coluna, denominada capacidade (disponibilidade), colocamos os 
valores correspondentes à disponibilidade que cada fábrica tem. Na última linha, registramos 


os valores correspondentes à demanda (necessidade) de cada revenda. 
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Destinos Rio de Capacidade 
Califórnia Teerã : Madri (disponibili- 
Origens cii dade) 
"" e 
Amsterdà 70 
Cairo 105 
E E 
Nova Iorque 130 
q 65 i 
México i 
370 * 
p 50 80 150 90 
(necessidade) 370 


Observação importante: Veja que a última célula fechou em 370/370 (que são os somató- 
rios de demandas e ofertas). Quando os somatórios ficam iguais, dizemos que o sistema está 
equilibrado. 


Agora, alocamos os valores de modo a obter a solução inicial. Basta alocar o máximo possível 


de unidades (comecando pela célula de menor custo e seguindo para as próximas células de 
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menor custo), sempre respeitando as capacidades de cada fábrica e a demanda de cada de- 
pósito. Nesse exemplo, faremos o seguinte: 


e Primeiro, colocamos 70 unidades na célula (1,4), cujo custo de transporte é igual a 2. 
Observe que essa é a célula de menor custo. 


e No próximo quadro, para nào nos perdermos, riscamos o valor 70. Observe que 70 é 
a capacidade máxima referente à primeira linha. Como a demanda para essa coluna é 
de 90 unidades, precisamos ainda transportar 20 unidades. O menor custo correspon- 


dente da coluna está na célula (4,4), cujo valor é igual a 5. 


e Portanto, basta alocar nela as 20 unidades restantes, o que faz a capacidade da quar- 
ta linha baixar para 45 (era 65 — 20 = 45). Para nào nos perdermos, estamos riscando e 
sempre colocando o novo valor, como vocé pode observar a seguir. 


Destinos : i 
Lees z Rio de : Capacidade 
Califórnia Teerã Janeiro Madri (disponibilidade) 


AH 


105 


Origens 


Amsterdã 


Cairo 
Nova Iorque 


E 
y 


© 


= A 
m 
Q 
© 


Cidade do nm E 0 Ps 
México x. 45 
Demanda 50 80 150 90 370 


(necessidade) Pd 370 


Completando novamente o quadro, a próxima célula de menor custo é a (4,1), com um custo 


unitário de transporte igual a 4. Vocé já deve ter observado que sempre fazemos o novo cál- 
culo partindo da linha ou coluna da última iteração realizada. A demanda referente a essa 
coluna é de 50 unidades, mas a capacidade da fábrica, nesse momento, é de apenas 45 unida- 
des. Note que, das 65 unidades disponíveis, 20 já estão alocadas na célula (4,4): 65 — 20 = 45. 


As 5 unidades que faltam para atender à demanda devem ser alocadas na célula (3,1), cujo 


custo unitário de transporte é 7, deixando nosso quadro com a seguinte composição: 
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Destinos 


Califórnia Teerà dns de Madri Capac iga "E 
Origens Janeiro (disponibilidade) 
Amsterdà ES P 
0 
Cairo 105 


Nova Iorque 


7 
E 130 


Qo 
- ES B 
E 


e] B B 
A 
[2i x 5 [2i 

= 


125 
Cidade do 10 p 
México 
20 
Demanda 0 150 sad 370 
(necessidade) 


370 


Como a capacidade é de 130 e conseguimos alocar apenas 5 motocicletas em decorrência de a 
demanda do depósito ser de 50 (50 — 45), podemos adicionar as 125 unidades restantes na célula 
(3,3) com um custo de 6 por unidade (lembre-se: menor custo). Seguindo esse raciocínio, basta 
preencher as demais células que ainda apresentam demanda e capacidade, até o momento em 
que todas as demandas estejam completas. Vejamos, então, o último valor que falta. 


Destinos ; i 
Califórnia Teerã Rig de Madri Capacidade 
Origens Janeiro (disponibilidade) 


Amsterdã 


Cairo 


© 


N 
pa 
Ha 
o 
q 


qq 


H 


Nova Iorque 


oo 


Ha 
XN 5 a 
m 


125 
Cidade do A 
México A 
Demanda n EM EE 
(necessidade) "m d 370 


Na tabela anterior podemos dizer que cada H |é uma celula matricial; desse modo, podemos rela- 


cionar ao elemento x,, ou seja o primeiro + [significa x, e assim sucessivamente. Mas voltemos 
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à resolução. Observe que as respectivas somas das linhas e das colunas, as quais correspondem 


ao total das capacidades e das demandas, é que determinam o gasto total de transporte. 


Dessa maneira, obtivemos a solução inicial para o problema, cujo objetivo era customizar o 
gasto de transporte, ficando assim o custo total (CT): 


CT = (70 - 2) + (80 - 5) + (25 - 10) + (5 - 7) + (125 - 6) + (45 - 4) + (20 - 5) = 
140 + 400 + 250 + 35 + 750 + 180 + 100 = 1.855 


O que precisamos, agora, é verificar se essa solução inicial é a melhor possível, ou seja, a solução 
que atende às necessidades dos depósitos respeitando-se as capacidades de cada fábrica, com um 
custo total de transporte que seja o menor possível, ou se ainda é possível realocar determinadas 
quantidades para poder, então, minimizar o custo total referente ao transporte das motocicletas. 


O processo adotado para verificar se a solução inicial é a melhor possível consiste em calcular os 
valores das variáveis u, e Vy conhecidas como variáveis duais. As variáveis u, estão relacionadas às 
origens, e as variáveis v, estão relacionadas aos destinos. 


Temos, então, um problema dual associado ao problema inicial (denominado, também, pro- 
blema primal): 


max z = 70u, + 105u, + 130u, + 65u, + 50v, + 80v, + 150v, + 90v, 


s.a. u, +v <8 
u, +v, <8 
u, +v, <5 
u, +v, <2 

u, +v <10 
u, +v, <5 
u, +V, <10 
u, *v, €10 
u, +v < 
u, +V, <10 
u, HV <6 
u, +v, <8 
u,+V, <4 
u, +V, <10 
u,+ +V, <10 
u, + +v,<5 


uy V, irrestritos de sinal 
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A base para um problema de transporte é formada por m + n - 1, vetores linearmente 
independentes. 


Note que, para cada variável básica do problema inicial, temos sempre uma restrição de igual- 
dade para o problema dual. Dessa maneira, para qualquer solução básica de um problema de 


transporte, utilizamos um valor u, ou M conhecido. Temos: 
u, t V, =C. 
1 J 1j 


Para as variáveis não básicas, os custos reduzidos d; podem ser calculados como segue: 


Para verificar se a solução é ótima ou nào, analisamos a seguinte condição: 


e se todos os custos reduzidos para as variáveis nào básicas forem maiores ou iguais a 


Zero, ou seja, GUY, 2 0, a solução é ótima; 


e no entanto, se existir algum custo reduzido negativo, ou seja, € -u - V, « 0, então, a 


respectiva variável não básica deverá entrar na base, resultando na troca dessa base. 
Importante: c, são os valores do custo encontrados na solução inicial. 
Cálculo dos valores de u, e ME 


Para facilitar o cálculo, esquematizamos a tabela a seguir, com a identificação das variáveis 
desconhecidas. 


x,,é a célula do custo de valor 2. Substituindo o i e o j das variáveis u, e v pori=1 ej=4 da 
variável x, temos u, e v,. Como precisamos de um valor conhecido, sempre sugerimos u,- 0. 
Vejamos o cálculo na tabela. 
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Para variáveis básicas, cuu, Ni 0. 
p/x,, 2-u,-v,-0 > p/u,=0 v,=2 
p/ X, 5-u,-v,-0 > 5-10-v,=0 | v,--5 
p/ Xy 10-u,-v,-0 > [10-u,-0-0 | u,-10 
p/ X, 7-u,-v,-0 — 7-u,-1-0 476 
p/ X; 6-u,-v,-0 > 6-6-v,=0 ,70 
p/ X, 4-u,-v,=0 = 4-3-v,-0 v,-1 
x 5-u,-2=0 E 
P/ Xw 5-u,-v,=0 > ae cat) u,=3 


Para o cálculo dos custos reduzidos, procedemos do mesmo modo, calculando as variáveis 
não básicas por d; =c, — u, - v,. As variáveis nào básicas são aquelas células da tabela que fi- 


caram em branco, ou Seja: X, X, Xiy Xa Xay Xay Xay Xy € X. 
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p/X, d,- 8-0-1 > d,- 7 
P/ Xiy d= 8-0- (25) > d,- 13 
p/ Xy d,- 5-0-0 > d,- 5 
P/X,y d, = 10-10-1 > d,- -1 
P/ X, d, - 10-10-2 E d. -2 
PI X», d, - 10-6 - (-5) > do= 9 
p/ Xay d,- 8-6-2 > d, - 0 
P/ X, d,- 10 - 3 - (-5) E d, - 12 
p/ Xy d,- 10-3-0 > d,- 7 


Como existem custos reduzidos negativos (-1 e —2), a solução inicial ainda não é a solução 
ótima. 


O que fazer, então, para chegar à solução ótima? 


Precisamos realizar alguns passos para melhorar a solução. Para determinar qual variável não 
básica entrará na base, consideraremos a que tem o custo reduzido mais negativo. Mas por 
que isso? Porque é essa variável que representará uma maior economia no custo de transporte 
para a atual iteração. Como x, é a variável que irá entrar na base, devemos construir, na tabela 
dos custos, um percurso de avaliação formado inicialmente pela variável que está entrando na 
base e pelas variáveis básicas correntes, no qual cada vértice do circuito deverá coincidir com 
uma variável básica. Partindo da última iteração (ou última tabela do custo mínimo), temos: 


D ti » . 
od Califórnia Teerã Rio de Madrid Capacidade 
Origens Janeiro (disponibilidade) 
8 » 
Amsterdã 70 
70 

| 

Cairo | 105 
0 5 | | 
ri 


8 
Nova Iorque [5] 130 
aues [5] [5 
Cidade do 65 


México 45 20 
Demanda 
(necessidade) 50 80 150 90 
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Definido o caminho, alocaremos na célula (2,4) o quanto for possível, somando e subtraindo, 
respectivamente, as possíveis quantidades de tal maneira que, no percurso de avaliação, ne- 
nhuma restrição de oferta e demanda seja violada. Após esses cálculos, uma variável básica 


será reduzida a zero. Isso significa que ela deixará de ser básica. Nesse caso, é a célula x,,, em 


31” 
que temos 5 — 5 = 0, com custo igual a 7. 


Destinos i i 
TE Rio de . Capacidade 
ep janeiro | Mada | (as sponibilidade) 


Origens 


Amsterdã 70 
70 
10 
Cairo 105 
80 25-5=20 | 0+5=5 
ES 
Nova lorque 198 
5-5=0 125 +5 =130 
; 
Cidade do 65 
México 454-5250 pu 
15 
Demanda 
(necessidade) 50 80 150 90 


Portanto, temos agora as novas quantidades a serem transportadas. Precisamos verificar se 


essa é a solução ótima ou se ainda podemos minimizar os custos de transporte. 


Destinos É : 
Califórnia | Teerã à» Madrid Wh "- 


Origens Janeiro (disponibilidade) 


PID JIR 
Amsterdã 


70 


105 
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E 

70 

5 

z| [o] [s] 
130 


130 
Cidade do 4 65 
México 50 15 
Demanda 
(necessidade) 50 80 150 90 
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Fazendo novamente a análise das variáveis básicas por c,- u,- v = 0, temos: 


p/x,, 2-u,-v,-0 > p/u- 0, =2 
p/ X, 5-u,-v,=0 > v,=-3 
P/ X, 10-u,-v,-0 > v,=2 
p/x,, 10-u,-v,-0 > u,-8 
p/ Xs, 6-u,-v,-0 — u,=4 
p/ X, 4-u,-v,-0 > v,71 
p/x,, 5-u,-v,-0 — u,-3 


Como obtivemos novas variáveis básicas, precisamos calcular os custos reduzidos das vari- 


áveis nào básicas. 


Cálculo dos custos reduzidos: 


Para variáveis nào básicas, d; 7Cj-U,- V 


p/ X, d, = 8-0-1 > d= 7 
p/ X, d,- 8-0-(-3) > d,- 11 
P/ Xy d= 5-0-2 > d= 3 
p/ X, dy = 10-8-1 > d,- 1 
p/X,, d. um 7-4-1 > d,- 2 
p/ X, d, - 10-4- (-3) > d,- 9 
p/ X, du 8-3-2 > d= 3 
p/ X d, = 10-3- (-3) > d= 10 
p/ X, d, = 10-3-2 > dE 5 


Como todos os custos reduzidos são positivos, a solução obtida é a melhor solução possível. 


Logo, as quantidades a serem transportadas são aquelas que apresentaram o melhor custo 


reduzido de transporte, como podemos observar na tabela a seguir. 


Califórnia Teerã sien Madri 
Amsterdã — — — 70 
Cairo — 80 20 5 
Nova lorque = = 130 = 
| ANC PB l e|: 


O custo total de transporte será dado pela multiplicação entre quantidades a serem transpor- 


tadas e os respectivos custos unitários, isto é: 
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CT=2.70+5.80+10-.20+10-5+6-130+4.50+5-15 

CT = 140 + 400 + 200 + 50 + 780 + 200 +75 

CT=1.845 

Portanto, o custo total é menor do que o valor de 1.855 encontrado na solução inicial. 


Contudo, se você encontrar em um problema de transporte uma capacidade diferente da quan- 
tidade demandada, ou seja, um sistema não equilibrado, basta criar um depósito fictício (quan- 
do a capacidade for maior do que a demanda) ou uma origem fictícia (quando a demanda for 
maior do que a capacidade). Lembre-se, ainda, de que, para o depósito ou para a origem fictícia, 


os custos de transporte são iguais a zero. 


Cálculo pelo método de Vogel ou método das penalidades 


Para facilitar a comparação entre os métodos, usaremos o mesmo exemplo com o mesmo 


propósito, ou seja, minimizar o custo de transporte. 


Esse método se baseia em forçar o transporte que se encontra na linha ou na coluna que tem 
a maior penalidade e o menor custo. A penalidade é a diferença entre os menores custos que 


se encontram na linha ou na coluna da tabela. Observe a tabela a seguir. 


Destinos : i 
Califórnia Teerã Ear Madri Capacidade 
Origens Janeiro (disponibilidade) 
Amsterdã 70 
Cairo 105 
m 
Nova Iorque 130 
M 65 
México 
370 
Demanda 50 80 150 90 
(necessidade) 370 


A essa tabela acrescentamos uma nova linha e uma nova coluna que representarão as pena- 
lidades, determinadas pela diferença entre os menores custos. 
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5 
£ 
E 
E 
$ 
E 
3 
E 
E 
E: 
3 
2 
E 
E 
o 
E 
s 
E 
3 
E 
3 
g 
8 
g 
$ 
e 
E 
E 
g 
8 
j- 
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Destinos Rio de Capacidade | Penalidades 
Califórnia Teerà Madri (disponibi- | (menor custo 


Origens Teneno lidade) da linha) 


Amsterdã 70 5-223 


Cairo 105 8-523 


a 
e 


pá 
© 
N 
l 
jon 
[i 
ma 


130 
Nova lorque 
90 
Cidade do 9» 0 i -— 
México 45 
370 
o 50 80 150 90 
(necessidade) 370 
Penalidades 
(menor custo 7-423 8-523 6-521 5-223 
da coluna) 


Agora escolhemos a maior diferença na linha ou na coluna para iniciar o processo de pena- 
lizacáo, descarregando as necessidades ou disponibilidades no quadro. 


Como temos várias penalidades (linhas e colunas) iguais, ou seja, a maior diferença foi 3, re- 
solvemos escolher a última coluna (Madri), fazendo, assim, o descarregamento no menor 
custo (célula x,,). 


Destinos Rio de Capacidade 
Califórnia Teerã Tani Madri (disponibili- | Penalidades 
Origens ^a dade) 


Amsterdã 70 5-2=3 


[e] 


Cairo 105 


Nova Iorque 


= 
o 
oo 
l 
q 
LU 
[9] 


zi 
y 
N 
l 
O 
ll 
pa 


130 
Cidade do 1 0 65 -[ 
México r 
370 
e 9 80 150 — |90-70-20 


370 


Penalidades 7-423 8-523 6-521 5=2=3 
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$ 
8 
E 
e 
z 
E 
g 
E 
& 
E 
5 
B 
8 
E 
e 
E 
É 
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penalidades excluindo a quarta coluna. 


Dessa forma, anulamos a coluna que sofreu o descarregamento. Agora, calculamos as novas 


Destinos 


Origens 


Califórnia 


Teerà 


Rio de 
Janeiro 


Madri 


Capacidade 
(disponibili- 
dade) 


Penalidades 


B 


Amsterdà 


© 
© 


m 
e 
N 
oo 
l 
q 
[i 
[99] 


105 


i : i i ] 

Nova Iorque 130 : 
5-421 E 

Cidade do i i 65 
México £ 
Demanda AM 370 E 

: 50 80 150 È 
(necessidade) 20 370 s 
7 gênio -5-s |Nembed | 5-2-3 : 

Penalidades É 


Percebemos que as novas penalidades (em cinza escuro) também se repetem. Para facilitar, 


adotaremos entre as penalidades repetidas a que apresenta o menor custo. Assim, utilizare- 


ização da Editora InterSaberes. A violac: 


mos a primeira coluna, pois apresenta o custo com valor igual a 4. 


Destinos Rio de Capacidade 
Califórnia Teerã : Madri (disponibi- | Penalidades 
Origens xil lidade) 
Amsterdà i 
70 0 
Cairo 105 
L5] , | ? ] i 
Nova Iorque 130 
cidäde do E 65-50=15 
México 
50 
Demanda yd 370 
idad 80 150 20 
(necessidade) 0 370 
7-4=3 8-5=3 6-5=1 5-2=3 


Penalidades 


licação poderá ser reproduzida por qualquer meio ou forma sem a prévia 


E 
E 
E 
2 


Fazemos o devido carregamento na célula, ou seja, na localidade (4, 1)(quarta linha e primei- 
ra coluna), de 50 unidades. 


Indo para nova iteração (lembra-se desse termo no método simplex?), achamos as novas pe- 
nalidades. Como temos ainda a igualdade nas penalidades e nos custos de menor valor igual 


a 5, podemos escolher qualquer um ao acaso. 


Vamos escolher a disponibilidade que ficou com 15 unidades a serem entregues, com custo 


igual a 5, na célula x,, conforme a tabela a seguir. 


44^ 


Destinos "ue N Rio de Capacidade Ê 

. Califórnia Teerá . Madri (disponibi- | Penalidades 3 
Origens Janeiro lidade) : 
70 -2- i 

Amsterdà : i 

70 E 

8-523 E 

5 10 10 i 

Cairo 105 : 


7-6=1 
7 10 
L5] 130 
Nova Iorque 
Cidade do México 


ização da Editora InterSaberes. A violac: 


a 


50 ts) 


E A 370 E 

Demanda EN 8 
(necessidade) E 150 A " : 
5 E 

7-4=3 | 8-5=3 | 6-5=1 | Bored 1 

8-5=3 | 6-5-1 | 5-2-3 É 

Penalidades E 
8-5=3 | 6-5-1 | 5-2-3 i 

= 8-5=3 | 6-5=1 — f 


g 
E 
E 
2 


Calculamos agora a nova penalidade: 


Destinos DN i Rio de Capacidade 
. Califórnia Teerã . Madri (disponibi- | Penalidades 
Origens Janęiro lidade) 
P» 
70 E 
Amsterdà : 
70 o 
8-5=3 
5 10 10 
; 105 8-5-3 
Cairo 
25 8-5-3 
80 : 
8-5-3 i 
7-621 J 
n i 
7-621 é 
Nova Iorque 130 E 
7-621 Ê 
7-6=1 i 
ê 
cR : 
Cidade do 15 5-4=1 E 
México = E 
i 50 Js Ü 5-4=1 ] 
D d p d 150 a 370 E 
Dn EE 0 20 : 
necessidade 
4 4 5 370 
7-4-3 | 8-5-8 | 6-5=1/| 5-2-3 É 
; 7-4=3 8-5=3 | 6-5-1 | 5-2-3 E 
Penalidades F 
7-4=3 8-5-3 | 6-5=1 | 5-2-3 E 
EM 8-5-3 | 6-5-1 - d 
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Calculando a penúltima iteração, temos: 


Destinos DEM - Rio de i Capacidade . 
. Califórnia Teerã : Madri (disponibi- | Penalidades 
Origens Janeiro lidade) 
Amsterdã 0 i 
70 " 
SEDIS 
1 1l 
pé 8 E 5 5 3 
Cairo 25 
80 25 0 a 
SEDIS $ 
mn E i i : 
7-6=1 E 
Nova Iorque 130 E 
7-621 E 
7-62] E 
5-4-1 E 
4 10 10 5 Br E 
Cidade do 5-4=1 = 
"E 15 8 
México 5-4-1 E 
50 15 0 i 
9 E 
p sr 150 zy E E 
Demanda 0 0 125 2” E 
(necessidade) 5 370 i 
7-4=3 | 8-5-3 | 6-5-1 | 5-2=3 i 
] 7-4- 8-523 | 6-5=1 | 5-2=3 E 
Penalidades 8 
— 8-5=3 | 6-5=1 | 5-223 8 
— 8-523 | 6-5=1 — 


Partindo para a última iteração, alocamos as penalidades com valor igual a 1, no menor custo. 


O valor cuja penalidade é igual a 1 se refere aos valores que se encontram na célula x,,, cujo 


custo de transporte é igual a 8. Poderíamos, se quiséssemos, fazer a iteração pela célula x,,, 


pois também tem iteração igual a 1. 
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Destinos ra F Rio de Capacidade 
. Califórnia Teerã : Madri (disponibi- | Penalidades 
Origens Janeiro lidade) 
70 ES 
Amsterdà : 
70 o 
Cairo 25 


8-5=3 

8-5=3 

8-5=3 

80 25 0 Bem 
7 10 7-6=1 g 
vo 7-6=1 F 
Nova Iorque 16 Abe E 
5 : 
7-621 E 
5-421 E 
D) i 
Cidade do a 5-421 3 
México 5-4-1 - 
50 15 0 E 
50 80 150 » 1 
Demanda 20 370 E 
0 0 125 s $ 
(necessidade) 370 
-4= 8-5=3 | 6-5=1 | 5-2-8 E 
; —4= 8-5=3 6-5=1 | 5-2=3 E 
Penalidades $ 
— 8-523 6-5=1 | 5-2=3 E 
— 8-5=3 | 6-5-1 — E 


Logo, para melhor visualização, deixamos a célula x,, para a última iteração da tabela a seguir. 


B 
8 
E 
e 
z 
E 
g 
E 
& 
E 
5 
B 
8 
E 
e 
E 
É 
3 
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Capacidade 


Desti i 
l estinos Califórnia Teerã Rio de Madri (disponibili- | Penalidades 
Origens Janeiro dade) 
Es] E 
70 == 
Amsterdã 0 m 
70 — 
| 8-523 
Cairo 8-5-3 
80 25 0 
8-523 
EI JE: 
Nova Iorque a 
125 5 0 NA 
7-6=1 
Cidade do M 
México Polo 
xi 50 15 0 2 
Demanda 0 0 0 0 370 
(necessidade) 
5 
7-4=3 | 8-5-3 | 6-5=1 | 5-2=3 7-4=3 
7=4=3m|=8=5=80 |96-5=) | 5-2=3 7-4-3 
Penali 
enalidades 7-423 | 825-23/!| 6-521] | 5-2-3 — 
= 8-5=3 | 6-5=1 = E 


Assim, escrevendo de cima para baixo e da esquerda para direita, encontramos os valores 2, 
5, 10, 6, 8, 4, 5, que são os custos de transporte e os valores das quantidades, apresentados na 
mesma ordem: 70, 80, 25, 125, 5, 50, 15. Finalizando, precisamos multiplicar o valor de cada 
custo pelo seu respectivo valor de quantidade: 


CT=2.70+5.80+10-.25+6.125+8.5+4.50+5.15= 
CT = 140 + 400 + 250 + 750 + 40 + 200 + 75 = 
CT=1.855 


Essa é a solução inicial obtida pelo método de Vogel. Para determinar a solução ótima, basta cal- 
cular os valores de u, e v; do mesmo modo quecalculamos anteriormente, seguindo os mesmos 
passos adotados no método do custo mínimo. 


Cálculo pelo método do canto noroeste 
Por fim, veremos o último método proposto neste estudo, embora existam vários outros. 


Já que o nome é método do canto noroeste, vamos seguir à risca essa orientação desde a primeira 
iteração. 


ireitos autorais é crime estabelecido na Lei nº 9610/1998 e punido pelo art. 184 do Código Penal. 


8 
E 
e 
8 
: 
g 
E 
& 
E 
5 
E 
8 
8 
* 
E 
É 
8 
$ 
g 
E 
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Vocé deve estar se perguntando: mas como funciona isso se usamos a mesma tabela? 


Vamos à resolução e você verá que esse método é tão simples quanto os outros. Utilizaremos 
o mesmo exemplo que foi visto nos cálculos dos métodos anteriores e partiremos da tabela 
inicial (a seguir). 


Destinos : í 
Califórnia Teerã Rio de Madri Capac ica de 
Origens Janeiro (disponibilidade) 
[5] 
Amsterdã 70 
E 
Cairo 105 
E 
Nova Iorque 130 
e 65 
México 
370 
Demanda 50 80 150 90 
(necessidade) 370 


Vamos para a primeira célula da tabela, ou seja, para a célula que indica a origem de Amsterdã 
e Califórnia, pois é a primeira, está no canto e a noroeste, para descarregar a disponibilidade 
de 50, como pode ser visto a seguir. 
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Destinos A 5 
Califórnia Teerã fio de Madri Rana Ade 
Origens Janeiro (disponibilidade) 
i 
Amsterdã 
20 
m 
Cairo 105 
ES T i 
Nova Iorque 
pco 65 
México 
Demanda aid zh PU 
(necessidade) " 870 
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Como zeramos a coluna da demanda (necessidade), nào iremos mais trabalhar com ela, pois 
cumprimos o que deveríamos fazer. Veja que sobraram 20 unidades de capacidade (disponibi- 
lidade) para serem atendidas na primeira linha. Seguindo da esquerda para a direita, descarre- 
gamos a capacidade que sobrou na mesma linha. 


Destinos Rio de Capacidade 
Califórnia Teerà A Madri (disponibilida- 
Origens Janeiro de) 
Amsterdã 20 
20 0 
Cairo 


Nova Iorque 130 


= 
e 
e 


Cidade do 
México 


na 
Ha 
o 
q 


65 


o 


e 
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1 370 


o 


Demanda o 


(necessidade) 


q gi 


60 370 


Vocé percebe que conseguimos anular a primeira linha e a primeira coluna por meio do canto 
noroeste? E desse modo que iremos continuar fazendo, com todas as outras linhas e colunas, 


ização da Editora InterSaberes. A violac: 


até terminar todas. 


O novo canto noroeste se encontra na posição (2,2) (lembra-se dos índices de matrizes para 
todas as tabelas?). 


Destinos , ] i 
Califórnia Teerà od d Madri Capa ida Be E 

Ori Janeiro (disponibilidade) E 
rigens E 
B : 

Amsterdà é 
0 E 

Cairo f 


130 
Nova Iorque 
25 
Cidade do 
Mexico 65 


9f 370 


370 


Demanda 
(necessidade) 


d ge 150 
0 


60 105 


A última iteração da tabela apresenta como canto noroeste o quadrado que tem o valor do 
custo unitário igual a U$ 6, que corresponde à posição (3,3), ou seja, a célula de Nova Iorque 
e Rio de Janeiro, conforme a tabela a seguir. 


Destinos 


Califórnia | Teerã id Madri Capacigade 
Origens Janeiro (disponibilidade) 
T 
Amsterdã 20 
[U 20 0 
Cairo 


q 


A 
B E m 
N 
BG 
c 
à F 


q 
© 


Nova Iorque 


gı 
= 
a qo en) 
e 


pee e À 65 
Mexico 
65 
0 150 
50 
Demanda 0 105 E 370 
(necessidade) Q 


0 0 370 


Assim, ao final, o cenário apresenta os custos 8, 8, 5, 10, 6, 8, 5 e as quantidades 50, 20, 60, 45, 105, 
25, 65, proporcionando o seguinte custo de transporte: 


CT = (8 - 50) + (8 - 20) + (5 - 60) + (10 - 45) + (6 - 105) + (8-25) + (5 - 65) 
CT - 400 + 160 + 300 + 450 + 630 + 200 + 325 
CT=2.265 


Note que essa é apenas a solução inicial. Sabemos que a solução pode ser melhorada. Para isso, 
basta calcular os valores de u, e v, como fizemos primeiramente no método do custo mínimo e 


no método de Vogel, seguindo, assim, os critérios de deixar as variáveis básicas CG-U,-V;> 0. 


Vamos analisar os dados do problema a seguir e aplicar os cálculos do canto noroeste para 


chegarmos a uma solução ótima. 


A indústria de notebooks Mellow possui três fábricas localizadas nas cidades de Santa 
Maria, São João e Santa Fé. Cada uma dessas fábricas produz, por dia, 12.000, 6.000 e 
4.000 unidades, respectivamente. Os notebooks produzidos são enviados para quatro 
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distribuidores localizados em Monte Claro, Floresta Grande, Campina Verde e Flores. As 
demandas dos distribuidores são de 4.000, 8.000, 3.000 e 5.000 notebooks, respectivamente. 
Os custos de transporte são os seguintes: 


e Monte Claro Pimpa campina Flores 
Fábricas Grande Verde 
Santa Maria 70 74 62 62 
Sào Joào 64 62 68 72 
Santa Fé 68 65 64 66 


Determine quanto deve ser transportado de cada fábrica para cada depósito de modo que 
o custo total seja o menor possível. 


Resolução 


Para que os dados fiquem mais claros, vamos sintetizar a tabela. 


Solução inicial obtida pela regra do canto noroeste: 


E 
Maria 4.000 8.000 0 


E E 

São João 6.000 
E 

Santa Fé 


Produz 


12.000 


4.000 


22.000 


Necessidade 
4.000 8.000 3.000 5.000 2.000 22.000 


Número de variáveis básicas: m +n -1=7 
a) Cálculo dos valores de u, e Vi 


Para variáveis básicas, Q -u-vj- 0. 


i 
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P/ Xy 70-u,-v,-0 > p/u, - 0, v, 770 
p/X, 74-u,-v,-0 > v,=74 
p/ X, 62-u,-v,-0 — v,= 62 
p/ X, 68-u,-v,-0 — u,-6 
p/ x, 72-u,-v,-0 — v,7 66 
p/ Xy, 66-u,-v,-0 — u,=0 
p/ Xa, 0-u,-v,-0 = 5-0 
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b) Cálculo dos custos reduzidos 


Para variáveis nào básicas, d; =c 74- V 


P/ X, d= 62-0-66 > d= -4 
Pl Xy, de= 0-0-0 > d= 0 
p/ X,, d,- 64 —6— 70 > d,, -12 
p/ X, d,- 62-6-74 > d,- -18 
p/ Xy, d,- 0-6-0 > de= -6 
p/ X4, d,- 68 —0- 70 > d,- -2 
p/ X», do = 65-0-74 > d „= -9 
p/ Xs, d= 64-0-62 > d= 2 


Primeira iteração 
Variável que entra na base: x,,. 


Pelo percurso de avaliação, determinamos qual variável sai da base, ou seja: 


is ae E A T i 
E —. n — L i : 
po | 4.000 | 
2.000 2.000 a 

4.000 8.000 3.000 5.000 2.000 

que resulta em: 
é 

—. EN [nd ii j 
— —. — EN iod : 
po | 4.000 
2.000 2.000 E 

4.000 8.000 3.000 5.000 2.000 : 


Fim da primeira iteração. 

Seguindo esses passos, a solução ótima é dada por: 

x, = 2.000; x,, = 3.000; x,, = 5.000; x,. = 2.000; x, = 2.000; x,, = 4.000; x,, = 4.000 

z = 70 - 2.000 + 62 - 3.000 + 62 - 5.000 + 0 - 2.000 + 64 - 2.000 + 62 - 3.000 + 65 - 4.000 


z = 1.272.000 
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Problemas de árvore mínima consistem em situações nas quais o objetivo é interligar todos 
os nós de um grafo! com o menor custo possível. Aparecem com muita frequência quando 
precisamos, por exemplo: 


e interligar computadores de uma rede com o menor gasto em cabeamento; 
* conectar caixas eletrónicos de um determinado banco; 

e instalar um sistema de segurança nas residências de um condomínio; 

e interligar máquinas, equipamentos e computadores de uma empresa. 


O algoritmo mais utilizado para que se obtenha a solução de um problema de árvore míni- 
ma é o de Kruskal. Ele se caracteriza por ligar todos os nós de um grafo de modo que a soma 
dos custos referentes às conexóes entre dois nós seja a menor possível. 


O princípio básico do algoritmo de Kruskal consiste em, a partir de um grafo com pesos, is- 
to é, com custos associados a cada arco, a cada iteração, marcar o arco de menor custo. Caso 
haja formação de ciclo, eliminam-se os respectivos arcos. Esse processo deve ser repetido até 
que todos os arcos estejam marcados ou eliminados. 


Para explicar melhor, considere o seguinte problema: uma companhia de TV a cabo está ins- 
talando decodificadores nas residéncias de um condomínio. A figura a seguir ilustra as loca- 
lizações de cada decodificador e os respectivos custos, em reais. Vamos determinar como o 
técnico da empresa deve proceder para que nào apenas todas as residéncias fiquem interli- 
gadas, mas também o custo total de instalação seja o menor possível. 


Figura 5.4 — Localização das residências do condomínio 
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Nenhuma parte desta 


O primeiro passo é determinar qual é a conexão de menor custo. Observando a Figura 5.4, 
você pode ver que o arco C-E tem um custo igual a 4. Assim, marcamos esse arco. 


1 Um grafo G (V, A) é uma estrutura G formada por um conjunto não vazio V, composto por objetos denominados vértices e um 
conjunto A de pares não ordenados de V, chamados de arestas. 
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Figura 5.5 — Primeira iteração do método de Kruskal 


O passo seguinte é determinar o próximo arco de menor custo. Observe que os arcos A-B e 
E-F têm um custo igual a 5. Nesse caso, escolheremos, ao acaso, o arco A-B. 


Figura 5.6 — Segunda iteração do método de Kruskal 


O próximo arco a ser escolhido é o E-F. 


Figura 5.7 — Terceira iteração do método de Kruskal 
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Como nào podemos fechar um ciclo, devemos eliminar o arco C-F. 


Figura 5.8 - Quarta iteracáo do método de Kruskal 


O próximo arco a ser escolhido é o B-D ou o D-E, ambos com custo igual a 6. Escolheremos, 
ao acaso, o B-D. 


Figura 5.9 — Quinta iteração do método de Kruskal 


O próximo passo é escolher o arco D-E, também de custo igual a 6. 


alquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal. 


Nenhuma parte desta publicação poderá ser reproduzida por qui 


Figura 5.10 — Sexta iteração do método de Kruskal 


Note que é necessário eliminar alguns arcos. Primeiramente, eliminaremos o arco B-E. 


Figura 5.11 — Sétima iteração do método de Kruskal 


Como ainda temos a possibilidade de fechar um ciclo, eliminaremos o arco A-E. 
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Figura 5.12 - Oitava iteração do método de Kruskal 


Devemos fazer o mesmo para o arco A-C, pois também forma um ciclo. 


Figura 5.13 - Nona iteração do método de Kruskal 


Dentre os arcos restantes, o de menor custo é o D-G. Marcamos, entáo, esse arco. 
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Figura 5.14 - Décima iteração do método de Kruskal 


Para que nào se forme um novo ciclo, devemos eliminar o arco E-G. 


Figura 5.15 — Décima primeira iteração do método de Kruskal 
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Devemos fazer o mesmo com o arco F-G. 


E 
E 
E 
e 
E 
E 
E 
g 
E 
8 
5 
B 
g 
E 
* 
S 
E 
S 
£ 
g 
E 


107 


Figura 5.16 — Décima segunda iteração do método de Kruskal 


Como todos os arcos foram marcados ou eliminados, chegamos a uma solução ótima do pro- 


blema de árvore mínima. 


Figura 5.17 — Solução final 


Os nós a serem interligados são A-B, B-D, C-E, D-E, D-G e E-F, com um custo mínimo de 5 + 
6+4+6+7+5, que totaliza R$ 33,00, referentes à instalação dos decodificadores. 
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Para saber mais 


Se você deseja aprofundar seus conhecimentos sobre os fluxos de rede e algoritmos de 
otimização de soluções, sugerimos o estudo matemático realizado pelo professor Paulo 
Feofiloff. 


FEOFILOFF, P. Fluxo em redes. 2003. Disponível em: «http://www .ime.usp.br/-pf/flows/ 
mynotes/FluxoEmRedes.pdf>. Acesso em: 20 nov. 2009. 
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At Sintese 


A PO localizou nos problemas relacionados a transportes (de matérias-primas, produtos 
industrializados, pessoas) uma área de aplicação privilegiada, os fluxos de rede. Nessas 
circunstâncias, na busca pelo custo mínimo, encontramos nas formulações dos mode- 
los matemáticos uma ferramenta fundamental para responder de forma organizada às 
demandas. Como vimos, os métodos utilizados podem variar (o do custo mínimo, o de 
Vogel, o do canto noroeste, entre outros), mas o propósito permanece, isto é, encontrar a solu- 
ção ótima (z), aquela que viabilize o lucro da organização em um ambiente de competitividade. 


A Questões para revisão 


1. Uma locadora de automóveis se defronta com um problema de alocação resultante dos con- 
tratos que permitem que os automóveis sejam devolvidos em localidades diferentes daquela 
onde foram alugados. Atualmente, há três agências de locação (origens), respectivamente com 
25, 25 e 50 carros excedentes, e quatro agências (destinos), que necessitam, respectivamen- 
te, de 15, 20, 30 e 35 veículos. Os custos unitários de transporte são dados na tabela a seguir. 


Destinos 
Origens 1 2 3 4 Disponibilidade 
1 10 5 6 7 25 
2 8 2 7 6 25 
3 9 3 4 8 50 
15 20 30 35 
Demanda 


Determine o custo mínimo de transporte desses veículos. 


2. Uma empresa que produz refrigeradores possui três fábricas, localizadas em São Paulo, 
Curitiba e Belém do Pará. Os refrigeradores são transportados para quatro depósitos locali- 
zados em Salvador, Rio de Janeiro, Campinas e Porto Alegre e, em seguida, distribuídos para 
as lojas. As capacidades de cada fábrica, as demandas de cada depósito e os custos unitários 
de transporte são dados a seguir. 
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Destinos 
Riod Port 
Origens | Salvador iba Campinas dn Disponibilidade 
Janeiro Alegre 
São 
Paulo 5 3 2 10 80 
Curitiba 6 5 4 4 40 
Belém 
do Pará 9 ? 5 ? di 
40 80 30 60 
Demanda 


Determine quanto deverá ser transportado de cada fábrica para cada depósito de modo que 


o custo seja o menor possível. 

3. Resolva o seguinte problema de transporte: 

min z=70x, + 74x,, + 62x, + 62x, + 64x, + 62x, + 68x, + 72x,, 

sa X,*X,*X4* X, < 16.000 
X, + XX AE x, < 8.000 


4.000 


x 
+ 
x 
IV 


x 
i 
x 
IV 


8.000 


3.000 


x 
2 
x 
IV 


x ix > 5.000 


X; > 0,inteiros 


Em que: x, é a quantidade de pneus a serem transportados da origem (fornecedor) i para o 
destino (localidade) j. 
4. Quais fatores identificam um sistema em equilíbrio ou em desequilíbrio? 


5. Por que devemos utilizar modelos matemáticos para solucionar problemas de logística? 


6. O grafo a seguir ilustra as localizações das residências de um condomínio fechado e as res- 
pectivas distâncias entre elas. Deseja-se interligar todas as casas para a instalação de um sis- 
tema de segurança. Determine quais conexões devem ser feitas para que o custo total refe- 


rente a elas seja o menor possível. 
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7. Um determinado banco está investindo na implantação de diversos caixas eletrônicos em 
uma pequena cidade, cuja estrutura é relativamente precária; portanto, será necessária a co- 
locação de dutos e cabos para interligar todos os caixas eletrônicos. O diagrama a seguir ilus- 
tra as localizações dos caixas eletrônicos, bem como os custos, em reais, relacionados aos ca- 
bos e aos dutos. 


va. D 476 


Determine quais conexões devem ser feitas de modo a minimizar o custo total referente à im- 


plantação desses caixas eletrônicos. 
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capitulo 


Utilização do WinQSB na resolução 
de problemas de otimização em redes 
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Conteúdos do capítulo 


e Como utilizar o WinQSB na resolução de problemas de otimização em redes. 


Após o estudo deste capítulo, você será capaz de: 


1. resolver problemas de otimização em redes utilizando o WinQSB. 


No Capítulo 4, você viu como resolver problemas de programação linear (PL) por meio do 
WinQSB. Agora, saberá como usar esse programa para resolver problemas de otimização em 
redes. Os procedimentos serão muito parecidos com aqueles abordados anteriormente. 


6.1 Utilizando o WinQSB na resolução de problemas de transporte 


Problemas de transporte podem ser resolvidos facilmente por meio do WinQSB. Para isso, cli- 
que em Iniciar, Programas e WinQSB. Depois, entre as opções apresentadas, escolha Network 


Modeling. 
Figura 6.1 - Menu do WinOSB 


Acceptance Sampling Analysis 
Aggregate Planning 

Decision Analysis 

Dynamic Programming 

Facility Location and Layout 
Forecasting and Linear Regression 
Goal Programming 

Inventory Theory and System 

Job Scheduling 

Linear and Integer Programming 


MarKov Process 


ZHBEBBXZEZ2ZZHBEEBNBEDZEH 


Material Requirements Planning 
Network Modeling 

Nonlinear Programming 

PERT CPM 

Quadratic Programming 

Quality Control Chart 

Queuing Analysis 

Queuing System Simulation 
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Após esse procedimento, aparecerá a tela representada na figura a seguir. 
Figura 6.2 — Tela de abertura da opção Network Modeling 


* Network Modeling BOSE 


Network 
Modeling 


Version 2.00 
Copyright O Yih-Long Chang 


Licensed to: User 
Company User 


Please wait while loading program 


Após a inicialização, basta acessar o menu File e selecionar a opção New Problem. 


Figura 6.3 - Menu File, opção New Problem 


Network Modeling 


à New Problem 
Load Problem 


Exit 


Observe que há vários problemas passíveis de serem resolvidos. Para a resolução de problemas 
de transporte, escolha a opção Transportation Problem. Você precisa, ainda, escolher a opção 


Minimization, pois nosso objetivo é minimizar o custo total de transporte. Escolha também a 


forma matricial para a entrada dos dados; para isso, selecione a opção Spreadsheet Matrix Form. 
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Figura 6.4 — Tela das especificações do problema 


NET Problem Specification 


E necessário que vocé escolha um título e informe o nümero de origens (Number of Sources) e 
o de destinos (Number of Destinations). Para isso, relembre o problema das motocicletas vis- 


to no Capítulo 5 a fim de usá-lo como exemplo para demonstrar a utilização do WinQSB na 
resolução de problemas de transporte. O objetivo era determinar o custo total de transporte 
das motocicletas produzidas nas quatro fábricas da empresa para os quatro depósitos. Assim, 
o nümero de origens é igual a 4 e o nümero de destinos também é igual a 4. Agora, é só cli- 
car em OK. 


Figura 6.5 — Tela das especificações do problema com os campos preenchidos 


NET Problem Specification 


[ 


9) 


Após o preenchimento dos dados iniciais, surgirá a tabela a seguir, com os espaços destina- 
dos aos custos unitários de transporte de cada origem para cada destino, bem como as capa- 
cidades (Supply) de cada origem e as demandas (Demand) de cada destino. 
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Figura 6.6 — Tabela para a inserção dos dados do problema 


Como inicialmente são apresentados nomes-padrão para as origens, você pode alterá-los cli- 
cando em Edit e em Node Names. 


Figura 6.7 — Alteração dos nomes dos nós 


Node Names for Motocicletas 


Coloque nas quatro primeiras posições as localidades das fábricas e, nas últimas, as locali- 
dades dos depósitos. 


Figura 6.8 - Nomes das localidades do problema 


Node Names for Motocicletas 


Criamos, então, uma nova tabela, com as localidades das fábricas e dos depósitos, como po- 
de ser visto na Figura 6.9. 
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Figura 6.9 — Tabela para a inserção dos dados do problema com os nomes das localidades 


From + To Califórnia Teerã Rio de Madri Supply 


O próximo passo é preencher a tabela com os custos unitários de transporte, as capacidades 
e as demandas. Use os mesmos dados apresentados no Capítulo 5. 


Figura 6.10 — Tabela preenchida com os dados do problema 


Após informar esses dados ao WinQSB, você obterá a resposta clicando em Solve and Analyze 


e, em seguida, em Solve the Problem. 


Figura 6.11 — Solução do problema 


A solução ótima consiste em transportar 70 unidades de Amsterdã para 
Madri, 80 unidades do Cairo para Teerã, 20 unidades do Cairo para o Rio de Janeiro, 5 unida- 
des do Cairo para Madri, 130 unidades de Nova Iorque para o Rio de Janeiro, 50 unidades da 
Cidade do México para a Califórnia e 15 unidades da Cidade do México para Madri, a um cus- 
to total de U$ 1.845,00. 


6.2 Utilizando o WinQSB na resolução de problemas de árvore mínima 


Veremos agora como resolver problemas de árvore mínima por intermédio do WinQSB. Os 
passos iniciais são bastante semelhantes aos utilizados para a resolução de problemas de 
transporte: clique em Iniciar, Programas, WinQSB e Network Modeling. A diferença é que ago- 


ra você deve selecionar a opção Minimal Spanning Tree. O problema é de minimização, e os 
dados são simétricos. 
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Figura 6.12 - Tela das especificações do problema 


NET Problem Specification 


ido pelo art. 184 do Código Penal. 


Como exemplo, utilizaremos o problema, já abordado no capítulo anterior, referente à inter- 
ligação de pontos de TV a cabo nas residências de um condomínio. As localizações das resi- 
déncias e os respectivos custos de instalacào em reais podem ser vistos a seguir. 


Figura 6.13 — Localização das residências do condomínio 
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Como o problema apresenta 7 residências, registre essa quantidade no campo Number of Nodes, 
destinado à colocação do número de nós do problema. 
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Figura 6.14 — Tela das especificações do problema com os campos preenchidos 


NET Problem Specification 


Clicando em OK, você verá a tabela a seguir para a colocação dos dados do problema. 


Figura 6.15 — Tabela para a inserção dos dados do problema 


E interessante alterar os nomes dos nós. Para isso, basta clicar em Edit e depois em Node Names 


a fim de fazer a devida alteração. 


Figura 6.16 — Alteração dos nomes dos nós 


Node Names for TV 


Após a alteração dos nomes, clique em OK. 
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Figura 6.17 - Nomes dos nós alterados 


Node Names for TV x 


Em seguida, é necessário registrar os custos apresentados na Figura 6.13. Para que todos os va- 
lores possam ser colocados corretamente, você deve escrever todas as conexões feitas a partir 

de cada um dos nós. Na linha referente ao nó A, insira os valores 5, 7 e 11 nas colunas referentes 

aos nós B, Ce E, respectivamente. Na linha referente ao nó B, insira os valores 5, 6 e 9 nas colu- 
nas referentes aos nós A, D e E, respectivamente. Como na tela das especificações do problema 

a opção Symmetric Arc Coefficients está marcada, quando você preencher o valor 5 na linha do 

nó A, coluna do nó B, automaticamente será preenchido o valor 5 na linha do nó B, coluna do 

nó A, e assim por diante. No entanto, é importante conferir sempre se os valores são colocados 

corretamente. Faça o mesmo para os demais nós até completar a tabela. 


Figura 6.18 — Tabela preenchida com os dados do problema 


Para obter a resposta, você precisa clicar em Solve and Analyze e, em seguida, em Solve the Problem. 


Figura 6.19 — Solução do problema 


Como solução, temos as seguintes conexões que geram o menor custo possível: A-B, E-C, B-D, 
D-E, E-F e D-G. O custo mínimo é de R$ 33,00. 
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q Sintese 


Neste capítulo, observamos como proceder para obtermos a solução de problemas de trans- 
porte e de árvore mínima por meio do WinQSB. O uso de softwares auxilia muito na resolu- 
ção de problemas de otimização em redes e, de modo mais abrangente, de problemas de PO. 


Questões para revisão 


1. Uma montadora possui duas unidades, uma em São Paulo e outra em Curitiba, com capa- 
cidade de produção de 25 mil e 20 mil unidades, respectivamente. Em virtude de uma gran- 
de demanda de exportação, é necessário enviar os automóveis produzidos para três portos: 
12 mil unidades para Santos, 16 mil para Paranaguá e 8 mil para Itajaí. Com base nos custos 
unitários de transporte apresentados a seguir e nas capacidades e demandas citadas, defina 
quantas unidades devem ser mandadas de cada unidade para cada porto de forma a manter 
o custo total referente ao transporte o menor possível. 


Santos Paranaguá Itajaí 
Curitiba R$ 400,00 R$ 150,00 R$ 280,00 
São Paulo R$ 110,00 R$ 380,00 R$ 410,00 


2. Uma grande rede de móveis e eletrodomésticos fará uma megapromoção de roupeiros. A 
empresa tem duas centrais de distribuição, uma em Colombo e outra em Campo Largo. As 
lojas da região estão situadas nas seguintes cidades: Curitiba, Almirante Tamandaré, Pinhais, 
Campo Largo e Araucária. Os custos unitários de transporte, capacidades e demandas são 
fornecidos a seguir. 


Custos unitários de transporte: 


Gombe o RS ppb CPO 


: Araucária 
Tamandaré Largo 


Colombo R$ 12,00 R$ 7,00 R$ 15,00 R$ 20,00 R$ 20,00 
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Campo 


R$ 19,00 R$ 24,00 R$ 22,00 R$ 5,00 R$ 11,00 
Largo 
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Capacidades: 


Origem Capacidade 
Colombo 1.300 
Campo Largo 860 
Demandas: 
Destino Demanda 
Curitiba 700 
Almirante Tamandaré 300 
Pinhais 400 
Campo Largo 540 
Araucária 310 


Determine quanto deve ser transportado de cada origem para cada destino de modo que o 
custo total de transporte seja o menor possível. 


3. Considere uma indústria de implementos agrícolas que tem, em uma determinada re- 


giáo, dois depósitos (A e B) e três revendas (P, Q e R). As capacidades de cada depósi- 
to, as demandas de cada revenda e o custo unitário do transporte dos implementos de ca- 
da origem para cada destino são demonstrados na figura a seguir. Determine, utilizando o 
WinOSB, quanto deve ser transportado de cada depósito para cada revenda de modo que o 
custo total de transporte seja o menor possível. 
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4. Encontre a árvore mínima associada à seguinte rede: 


E 22 
A 10 di H 
19 
E 
18 14 17 
G 12 
B m~~ 16 I 
16 D 
N 27 
F 


5. A figura a seguir mostra a localização dos computadores e equipamentos de uma empresa, 
bem como as distâncias entre eles. 


A E D 12 19 I = M 28 O 
15 
10 10 G H 11 11 
22 22 
B E 12 19 J N 
10 10 17 
22 
c SE K 
18 
È 


Determine quais conexões devem ser feitas de maneira a minimizar o custo total referente a 
cabos de rede necessários para conectar todos os computadores e equipamentos. 
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Conteüdos do capítulo 
e Efeitos na solução ótima quando ocorrem mudanças nas quantidades ou valores de um 
modelo de otimização. 


e Quando e por que utilizamos a simulação. 


Após o estudo deste capítulo, vocé será capaz de: 

1. verificar se alterações provocadas em restrições ou recursos de um problema produzem 
alterações e influências na solução ótima; 

2. identificar quando ocorre a sensibilização das variáveis analisadas; 


3. entender o processo de simulação. 


Usamos a análise de sensibilidade com a finalidade de verificar se algumas alterações intro- 
duzidas em determinados coeficientes de um problema (restrições ou recursos) de otimiza- 
ção — max (z) ou min (z) - alteram ou influenciam a solução ótima. O objetivo é saber em qual 
momento as variáveis analisadas se sensibilizam com algumas modificações que fazemos 
para levantar a hipótese de certeza de confiabilidade nos valores apresentados pelo modelo. 


7.1 Aplicações da análise de sensibilidade 


As modificações de que vamos tratar aqui acontecem acima ou abaixo dos valores já conheci- 
dos dos coeficientes da função objetivo ou nos valores dos coeficientes das restrições (recur- 
sos). Nesse contexto, Lachtermacher (2009, p. 97) avalia que a análise de sensibilidade deve 
responder a três perguntas: 


e Qualéo efeito de uma mudança em um coeficiente da função objetivo? 
e Qualéo efeito de uma mudança em uma constante de uma restrição? 
e Qualéo efeito de uma mudança em um coeficiente de uma restrição? 


O autor menciona que uma maneira simples de responder a essas perguntas é fazer as alte- 
rações na modelagem e encontrar nova solução por meio de novas otimizações. 


Pode ocorrer que, ao se fazer uma modificação nos valores ou nas quantidades da modela- 


gem, nada aconteça, não exercendo essa alteração influência no valor final. 


Para demonstrarmos isso de forma prática, veremos esse processo no problema a seguir. Nele, 
apresentamos a situação de uma empresa que fabrica dois produtos de limpeza e cujos gesto- 


res querem saber qual quantidade de produção irá maximizar o lucro da empresa. 


Temos, portanto, dois produtos de limpeza: o Limpa-Limpa e o Tira-Tudo, ambos fabricados 
pela empresa Vem Que Tem Ltda. Na sua produção, as quantidades e as características dos 


recursos utilizados correspondem ao expresso na tabela a seguir: 
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Produtos/recursos Limpa-Limpa Tira-Tudo 


Matéria-prima I 2un. 4un. 


Matéria-prima II 5un. 5un. 


Por sua vez, o estoque das matérias-primas I e II (recursos utilizados na produção) é de exa- 
tamente 400 e 600 unidades, respectivamente. Sabemos, ainda, que o retorno financeiro, ou 
o lucro, dos produtos Limpa-Limpa e Tira-Tudo é de R$ 10,00 e R$ 14,00, respectivamente. 


Esse é o contexto em que vamos observar a análise de sensibilidade, lembrando que os gesto- 


res querem saber qual quantidade de produção irá maximizar o lucro da empresa. 


Se resolvêssemos pelo método simplex, a última iteração, ou tabela final, ficaria assim: 


X, E E 
z=1.560 0 0 d = 
40 0 1 E 2 
100 1 0 E i 


Nessa iteração final, os coeficientes da função objetivo que se referem às variáveis de folga x, 
= E (ou 0,6, que é o resultado da divisão) e x, = J (ou 2,2) são chamados de valores duais ou 
shadow prices. Essa é a denominação utilizada por vários autores que se dedicam ao estudo 
de pesquisa operacional (PO). 


Você deve estar se perguntando: mas o que significa shadow prices? 


E a informação “preço sombra”. Esse preço é o valor de aumento no custo marginal de cada 
unidade. Lembre-se de que em um processo de produção o custo marginal corresponde ao 


acréscimo feito no n + 1, ou seja, no produto fabricado após uma produção inicial de n produtos. 


Podemos resumir esse processo do seguinte modo: 
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e a cada aumento de uma unidade de matéria-prima I utilizada, teremos um aumento 
de 0,60! (3) no lucro; 


e coma matéria-prima II, verificamos o mesmo processo, ou seja, quando há o consumo 
de uma unidade de produção a mais, o lucro gerado será de 2,20. 


E se, em vez de aumentarmos a unidade consumida de matéria-prima (I ou ID, nós a 
diminuirmos? 


1 Para facilitar os cálculos, utilizaremos os valores como unidade monetária, sem especificar a moeda. 
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Ora, se nós a diminuirmos, o efeito será o contrário, ou seja, ao diminuirmos a unidade de 
matéria-prima (I ou II), em virtude de uma customização no processo de produção, sem in- 
fluenciar totalmente a qualidade, deixaremos de ganhar os 0,60 ou 2,20, por unidade. 


Agora, se, em vez de 400, tivéssemos conseguido mais 100 unidades da matéria-prima I com 


nosso fornecedor, passando nossa restrição a ser < 700? 

Lembre-se de que você já está visualizando o problema na forma-padrão: 
Função objetivo: max z = 10x, + 14x, 

Restrições: 

2x, + 5x, < 400 

4x, + 5x, € 600 

x, X, 20 


Se mudarmos a primeira restrição, ficaremos com 2x, + 5x, € 700. Calculando tudo de novo, 


a nossa nova tabela final apresentará o seguinte resultado: 


x x X, x; 

z = 1.680 E 0 0 = 

FA 100 -2 0 1 -1 
" 120 B 1 0 i 


Veja que, nesse quadro, a produção do produto Limpa-Limpa (x,) deixa de ser fabricado. Se 
considerarmos que não podemos deixar de fabricar tal produto, o problema fica insustentável. 


Na análise de sensibilidade, para Passos (2008, p. 160, grifo nosso), devemos tomar cuidado, pois 


não podemos fazer crescer ou decrescer aleatoriamente os valores dos termos indepen- 
dentes (recursos) e manter a otimização do modelo. Existe um limite superior e inferior 
de acréscimo ou decréscimo nos valores, denominado intervalo de variação, onde o valor do 


preço sombra é mantido e, assim, o lucro permanece otimizado. 


Podemos, portanto, afirmar que é importante não perder o controle dos ajustes nos modelos. 
Quando ocorre a perda do controle, o processo fica livre para produzir alterações com valores 


quaisquer. E, quando isso acontece, dizemos que o modelo se tornou degenerado. 


Então, partimos para o recurso II, ou seja, aumentamos a matéria-prima II em 700 unidades 
e mantemos a primeira com o respectivo valor original. 
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A segunda restrição passou de: 
4x, + 5x, < 600 para 4x, + 5x, < 700 


Se resolvermos pelo método simplex, teremos como iteração final o seguinte quadro: 


x, x, x, X, 

3 11 

z= 1.780 0 0 5 5 

1 

x, 150 1 0 1 2 
x; 20 0 1 Z E 


Podemos observar que z = 1.780. Aumentou, portanto, o lucro. Em relação à mudança inicial, 
o preço sombra permaneceu com os mesmos valores de 0,6 e 2,2. 


Nesse processo, é interessante você atentar para as observações feitas por Santos (2009, p. 2), 
o qual resumidamente apresenta uma análise conclusiva sobre as dimensões das alterações 
nos quadros do modelo simplex: 


1. Uma alteração de um coeficiente de uma variável não básica da função objectivo apenas 


implica a modificação do coeficiente correspondente do quadro de simplex; 


2. Uma alteração de um coeficiente de uma variável básica da função objectivo apenas 


implica a modificação da última linha do quadro de simplex; 


3. Uma alteração dos termos independentes das restrições apenas implica a modificação da 


última coluna do quadro de simplex; 


4. Uma alteração dos coeficientes de uma variável não básica das restrições apenas implica 


a modificação da coluna correspondente do quadro de simplex; 


5. Uma alteração dos coeficientes de uma variável básica das restrições implica a modifica- 


ção de todo o quadro de simplex. 


A análise de sensibilidade ocorre por meio de dois tipos fundamentais. De acordo com 
Lachtermacher (2009), um se caracteriza pelo fato de avaliar a possibilidade de alterações e 
influências quando ocorre apenas uma alteração por vez na otimização do problema; o outro, 
por avaliar as alterações e influências na solução ótima de um problema quando ocorre mais 


uma alteração de forma simultânea. 


No primeiro caso, explica o autor, são estabelecidos “limites superiores e inferiores para to- 
dos os coeficientes da função-objetivo e para as constantes de restrição” (p. 97). Um aspecto 
interessante é que essa análise é realizada de forma automática pelo Microsoft? Office Excel 
e também pelo WinQSB . 
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No segundo caso, como é um procedimento mais complexo, nào é feita a análise de forma 
automática e precisamos interferir pelo processo de modelagem, procurando, mediante al- 
terações no problema, “nova solução por meio de uma nova otimização”, como foi visto nos 


quadros anteriormente demonstrados. 


Você deve lembrar que essas análises são realizadas após o processo de otimização inicial, vi- 
sando justamente avaliar a confiabilidade em relação aos coeficientes da função objetivo e às 
constantes de restrições. 


Para saber mais 


Essas alterações, assim como os problemas, podem ser expressas pela representação grá- 
fica da análise de sensibilidade. Se você tiver interesse em aprofundar seus conhecimen- 


tos sobre esse assunto, sugerimos a leitura do seguinte livro: 


LACHTERMACHER, G. Pesquisa operacional na tomada de decisões. 4. ed. São Paulo: 
Prentice Hall, 2009. 


Assim, se, por um lado, precisamos de um embasamento matemático (o Capítulo 8 irá auxi- 
liá-lo na rememoração dos cálculos matriciais), por outro, não podemos esquecer que o nos- 
so foco é a tomada de decisão, é a gestão. E a análise de sensibilidade é uma etapa de funda- 
mental importância para fazermos a análise da tomada de decisão. Nesse sentido, Eschenbach, 
citado por Silva e Belderrain (2009), elucida que a análise de sensibilidade tem uso efetivo 
nos processos de: “(1) tomar melhores decisões, (2) decidir quais dados estimados devem ser 
refinados antes de tomar uma decisão e, (3) concentrar-se nos elementos críticos durante a 
implementação”. 


Veremos, na sequência, como realizar a análise de sensibilidade em cada um dos elementos 
vistos há pouco: coeficientes da função objetivo, constantes das restrições e coeficientes das 
restrições. 


Alterações nos coeficientes da função objetivo 


Para analisarmos o que acontece com um problema de programação linear (PL), quando mu- 
danças são feitas na função objetivo, iremos utilizar um parâmetro « em cada coeficiente da 
função objetivo, separadamente, para determinar o intervalo dos possíveis valores de a, tais 
que as condições de otimalidade sejam satisfeitas. Duas abordagens simples, mas distintas, 
serão empregadas para os seguintes casos: 


e variáveis não básicas; 
e variáveis básicas. 


Com o intuito de ilustrar os procedimentos para que seja feita a análise de sensibilidade em 


um problema de PL, precisamos relembrar o exemplo da indústria de brinquedos, abordado 
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nos Capítulos 2 e 3. O problema foi resolvido pelo método simplex, resultando na tabela fi- 
nal reproduzida a seguir. 


Tabela final do problema da indüstria de brinquedos 


x, dis x: x X; 
z - 31.800 0 0 9 38 0 
x, 400 1 0 2 -2,66 0 
x, 1.800 0 1 -1 4,66 0 
x, 300 0 0 -2 2,66 1 


Variáveis básicas: x, x, e x, 
Variáveis nào básicas: x, e x, 
Todos os valores necessários para determinarmos as desigualdades envolvendo o parámetro 


a podem ser obtidos diretamente da tabela final: 


e Para a variável não básica x, basta observar que o custo reduzido é 9 - a, pois 9 é o 
coeficiente da variável x, na função objetivo. 


e Como9-a >0 faz com que a solução obtida seja ótima, é fácil, então, perceber que, nesse 


caso, a solução do problema é ótima para qualquer valor de o entre —% e 9, ou seja, a <9. 


e Seaformaior do que 9, significa que a solução ótima ainda não foi encontrada. Nesse 
caso, x, entra na base e mais iterações serão necessárias para encontrarmos a solução 


ótima. 


Em resumo, se os lucros relacionados à quantidade de plástico excedente (a variável x, está 
relacionada à folga da restrição referente à quantidade de plástico disponível) forem maio- 
res do que R$ 9,00, é vantajoso para a empresa vender o excedente e, consequentemente, ma- 
ximizar o lucro total. 


Considerando agora a variável não básica x, referente à quantidade de alumínio excedente, 
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temos que seu coeficiente na função objetivo é igual a 38. Logo, o custo reduzido é 38 — a. Para 
que a solução obtida pelo método simplex seja ótima, a deve estar entre —oo e 38, ou seja, a x 38. 


Se a for maior do que 38, x, entra na base e a nova solução ótima precisa ser encontrada. Isso 


significa que, se o lucro referente à venda do alumínio excedente for maior do que R$ 38,00, a 
empresa poderá aumentar seu lucro total. 


Os cálculos referentes à análise de sensibilidade das variáveis básicas seguem o mesmo princípio, 
mas envolvem também os custos reduzidos das variáveis não básicas”. De modo resumido, as desi- 
gualdades em relação ao parâmetro a podem ser obtidas diretamente da tabela final. Os coeficientes 


2 Mais detalhes podem ser vistos em Zionts (1974, p. 154). 
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de a são obtidos na linha referente à variável básica em questão e nas respectivas colunas re- 


ferentes às variáveis nào básicas. Como assim? E fácil! Vejamos na sequéncia. 


Relembrando a solução ótima do problema da indústria de brinquedos, temos que as variáveis 
básicas são x,, x, e x,. Podemos considerar a análise feita para cada variável básica como segue. 


Para x, variável que representa a quantidade de caminhonetes a serem produzidas, temos a 


seguinte expressão: 
(9 38)+a(2 —2,66)20 


Observe que 9 e 38 se referem aos coeficientes de x, e x, respectivamente, na função objetivo. 
Os termos 2 e -2,66 são os coeficientes da linha referente à variável básica x, nas colunas re- 
ferentes às variáveis não básicas x, e x,. 


Resolvendo a desigualdade que mostramos, temos: 
(9 38)+a(2 -2,66)>0 

(9 38)+ (2a 2,660) > 0 

(9+2a 38-2,660) >0 

Que resulta em: 

9+20>0 

38 - 2,660 > 0 


Da inequação 9 + 2a >0, associada à variável não básica x, (pois os termos que aparecem ne- 
la são da coluna de x,), temos que: 


2a > -9 
-9 

a2 

a2 —4,5 
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Resolvendo, então, a inequação 38 — 2,66a > 0, associada à variável nào básica x, temos: 


—2,66a 2 —38 


Nenhuma parte desta 


2,660. < 38 


a < -38 
— 2,66 


a < 14,28 


Obtivemos inicialmente a condição de que a deve ser maior ou igual a —4,5. Depois, vimos 
que a também precisa ser menor ou igual a 14,28 para que a solução original do problema 
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seja ótima. Reunindo as duas informações, temos que a deve estar entre —4,5 e 14,28, ou se- 
ja, —4,5 < a < 14,28. 


Como o lucro referente à venda das caminhonetes é R$ 12,00 e —4,5 < a < 14,28, a solução ori- 
ginal é ótima se o lucro referente à venda das caminhonetes estiver entre 12 — 45 e 12 + 14,28, 
ou seja, se o lucro c, estiver entre R$ 7,50 e R$ 26,28 (-7,50 € c, < 26,28). 


Isso significa que uma nova solução ótima poderá ser obtida se o lucro por caminhonete for 
menor do que R$ 7,50 ou maior do que R$ 26,28. Se a for menor do que 4,5, x, entra na base. 
Se a for maior do que 14,28, x, entra na base. 


Considerando agora a variável x,, que representa a quantidade de modelos esportivos que 
serão produzidos, temos a seguinte desigualdade: 


(9 38)+a(-1 466)>0 


Note que 9 e 38 são os coeficientes das variáveis x, e x, respectivamente, e —1 e 4,66 são os 
coeficientes da linha referente à variável básica x, nas colunas referentes às variáveis não bá- 


sicas x, e X, 

A solução da desigualdade a seguir: 
(9 38)+a(-1 466)20 

é dada por: 

(9 38)+a(-1 466)20 

(9 38)+(-a 4,660) >0 


que resulta nas seguintes desigualdades: 


9-a20 
e 
38 44,660.20 


Resolvendo a inequação 9 - a > 0, temos: 


Da inequação 38 + 4,66a. > 0, segue que: 
38 + 4,660, > 0 


4,660. > —38 
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o2 298 
— 4,66 


a 2 -8,15 


Como a <9 e a >-8,15, podemos escrever que -8,15 <a < 9, ou seja, o parâmetro a está en- 
tre -8,15 e 9. 


Assim, como o lucro referente a cada miniatura do modelo esportivo é de R$ 15,00 e -8,15 < 
a € 9, temos que: 


e seo lucro estiver entre 15 — 8,15 e 15 + 9, ou seja, entre R$ 6,85 e R$ 24,00 (6,85 < c, < 
24,00), a solução que obtivemos inicialmente para o problema de PL é ótima; 


e mas, se o lucro for inferior a R$ 6,85 ou superior a R$ 24,00, a solução pode ser melho- 
rada e certamente mais iterações do método simplex serão necessárias para resolver- 
mos o problema de otimização da fábrica de brinquedos. 


E, finalmente, a análise de sensibilidade em relação à variável básica x, será feita como segue: 
(9 38)+o(-2 2,66) >0 


Note que 9 e 38 são os coeficientes das variáveis x, e x, respectivamente, e —2 e 2,66 são os 
coeficientes da linha referente à variável básica x, nas colunas referentes às variáveis não bá- 


sicas x, e x,. 

A solução da desigualdade é apresentada a seguir: 
(9 38)+a(-2 2,66)>0 

(9 38)+(-2a 2,660) 20 

que resulta em: 

9-2az0 


-2a 2-9 


38 + 2,660. 2 0 


2,660 2 —38 


5-98 
— 2,66 


a 2 1429 
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Nesse caso, a < 4,5 e a > -14,29. Podemos escrever que -14,29 < a < 4,5, ou seja, a está com- 
preendido entre -14,29 e 4,5. 


Portanto, a solução original é ótima se o estiver entre —14,29 e 4,5. Caso contrário, a solução 
pode ser melhorada. 


Mas o que significa isso? Relembrando, a variável x, se refere ao número de caminhonetes 
excedentes. Logo, a solução obtida é ótima para qualquer lucro entre R$ -14,29 e R$ 4,50 re- 
ferente às caminhonetes não vendidas. Qualquer valor abaixo de R$ -14,29 ou acima de R$ 
4,50 faz com que a solução não seja ótima. 


É importante ressaltar que, na análise de sensibilidade, quando dizemos que a solução original 
se mantém ótima, significa dizer que as variáveis básicas continuam básicas e que as variáveis 
não básicas são exatamente as mesmas. No entanto, é fácil perceber que, se o lucro referente 
a um determinado item aumentar e a produção for a mesma, o valor da função objetivo tam- 
bém aumentará. Se diminuirmos o lucro unitário de um produto, o valor da função objetivo 
também diminuirá. No caso de problemas de minimização, se reduzirmos o custo unitário de 
um produto, o custo total também diminuirá e, se aumentarmos o custo unitário de um item, 
o custo total aumentará. 


Considere o exemplo da indústria de brinquedos. Se a indústria desejar aumentar o preço do 
modelo esportivo de modo que o lucro unitário passe de R$ 15,00 para R$ 19,50, a solução 
original obtida pelo método simplex ainda será ótima? Nesse caso, o que irá acontecer com 


o lucro máximo? 
Resolução 


Como sabemos, se o lucro do modelo esportivo estiver entre R$ 6,85 e R$ 24,00, a solu- 
ção original será ótima. Como o lucro será de R$ 19,50, então, não haverá modificações na 
solução. No entanto, o valor da função objetivo, dado por z = 12x, + 15x,, passará de R$ 
31.800,00 (conforme vimos no Capítulo 3) para R$ 39.900,00, pois a nova função objetivo é 
z = 12x, + 19,50x,. Assim: 


z = 12 - (400) + 19,50 - (1.800) 
z = 4.800 + 35.100 


z = 39.900 


Os valores 400 e 1.800 se referem à produção de caminhonetes e modelos esportivos, respec- 
tivamente, conforme vimos no Capítulo 3. 
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Alterações nas constantes das restrições 


A análise de sensibilidade para as constantes das restrições é feita de modo análogo à análise 
feita para os coeficientes da função objetivo. Utilizando ainda o exemplo da fábrica de brin- 
quedos, vamos considerar, por exemplo, a restrição referente ao plástico (primeira restrição). 
A quantidade disponível de plástico passará a ser 1.000 + a. Assim, para obtermos as desi- 
gualdades referentes à primeira restrição, basta considerarmos o vetor dos termos indepen- 
dentes da solução ótima somado à coluna referente à variável de x, multiplicada por o, pois 
essa variável é a folga da primeira restrição, todos da solução ótima, como segue: 


400 2 

1.800 |+| -1 |-a € 0 
300 -2 

400 20, 

1.800 |+| -œ |>0 
300 -20 

400 +20 
1.800-0./>0 

300 — 201 


que resulta em: 
400 +20. >0 
1.800- 0.2 0 

300 +20. 2 0 
Logo, temos que: 
a 2 -200 

a € 1.800 

a < 150 


Dessas desigualdades, temos que a > -200 e que a < 150. Logo, -200 < a < 150. Isso signifi- 
ca que a solução obtida é ótima para a quantidade disponível de plástico entre 1.000 — 200 e 
1.000 + 150, ou seja, entre 800 kg e 1.150 kg. Se a quantidade de plástico que a indústria tem 
disponível for inferior a 800 kg ou superior a 1.150 kg, a solução não será ótima. 


É importante destacar que, quando dizemos que a solução original é ótima, queremos dizer 
que as variáveis básicas da solução original continuam na base. Mas é fácil perceber que, se 
alterarmos a quantidade disponível de uma das matérias-primas, também iremos alterar a 
quantidade a ser produzida. Se for possível produzir 400 caminhonetes e 1.800 modelos es- 
portivos com uma tonelada de plástico, se reduzirmos essa matéria-prima para 800 kg, por 
exemplo, a produção também será reduzida. 
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Para as demais restrições, o procedimento é sempre o mesmo: consideramos o vetor dos ter- 
mos independentes somado à coluna referente à variável de folga da restrição em questão, 
multiplicada por a. 


Alterações nos coeficientes das restrições 


A análise de sensibilidade nos coeficientes das restrições é importante quando há a possibili- 
dade de terem ocorrido erros nos dados do problema. Na prática, essa análise não é tão sim- 
ples, tanto sob o ponto de vista teórico quanto, principalmente, sob o ponto de vista compu- 
tacional. Neste livro, não iremos aprofundar o estudo sobre alterações nos coeficientes das 


restrições. Mas você poderá encontrar algumas orientações interessantes na abordagem feita 


por Zionts (1974, p. 161) sobre o assunto. 


Embora a simulação possa ser definida como “uma técnica que envolve a construção de um 
modelo de uma situação real para sua posterior experimentação” (Loesch; Hein, 2009, p. 207), 
dada a abrangência de tal conceito, queremos deixar bem claro que neste estudo, quando fa- 
lamos em simulação, estamos nos referindo a modelos matemáticos e lógicos na busca por 
uma solução ótima nos processos de otimização. Podemos, ainda, pensar em simulação como 
imitação de um processo real, realizado em um determinado momento. 


Por exemplo, podemos simular a linha de produção de uma indústria, para avaliar os efeitos 
de mudanças em um processo produtivo ou, ainda, simular sistemas de distribuição ou de 
controle de estoques, para melhorar o fluxo desses sistemas. 


Essa imitação teria como objetivo apresentar a ideia de um cenário representativo de uma si- 


tuação real, prevendo, assim, o comportamento, ou seja, como funcionaria. 


Muitas vezes usamos recursos computacionais para simular, em decorrência da grande quan- 
tidade de cálculos que são feitos, o que torna a simulação manual quase rara, por causa de 
sua complexidade. Em particular, simulações são utilizadas também para inferir distribui- 


ções amostrais de quantidades, valores ou de qualquer outro interesse, investigando-se, as- 
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sim, uma série de situações. 


A simulação é empregada como ferramenta de análise, quando queremos prever mudanças 
em situações já existentes ou como ferramenta para construção de novos cenários ou situa- 


ções, como já mencionamos. 


Você poderia perguntar: mas por que fazer simulação se posso utilizar a análise de 
sensibilidade? 


Na simulação, podemos criar novos procedimentos operacionais, situações de controle e de 
fluxos de comunicação sem alterar a situação real. Além disso, podemos testar novos equipa- 


mentos e sistemas de transportes, testar a aplicabilidade de novos procedimentos e possibilitar 
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a aceleração ou o retardo do tempo da situação desenhada. Resumindo, podemos entender 
melhor todo um sistema. 


Mas vale lembrar que algumas dificuldades podem aparecer na construção de uma simula- 
ção. Nesse sentido, o parecer de Santos (1999) é que podemos ter dificuldades para interpre- 
tar os resultados e perder mais tempo e consequentemente dinheiro, além de os resultados 
apresentados, na maioria das vezes, não serem exatos. O fato é que as simulações devem ser 


feitas por profissionais competentes e com muita experiência nas áreas em que atuam. 


Etapas e tipos de modelos de simulação 


Em estatística, existem vários resultados que podem ser ilustrados via simulação, o que aju- 
da na compreensão e visualização dos conceitos e dos resultados. Em nosso estudo, os pro- 
cedimentos matemáticos de simulação são desenvolvidos no âmbito das variáveis aleatórias. 
Nesse sentido, Loesch e Hein (2009, p. 208) apresentam quatro fontes para a obtenção de nú- 
meros aleatórios: os métodos manuais, as tabelas, os métodos para computadores analógicos 


e os métodos para computadores digitais. 
Os modelos de simulação podem ser: 
e estáticos - em que a passagem do tempo é irrelevante; 
e dinâmicos- cujo sistema varia com o passar do tempo; podem ser discretos ou contínuos; 
e  determinísticos — que não contêm variável aleatória; as variáveis são conhecidas; 
e estocásticos - que apresentam variáveis estimadas, por médias. 


Vamos, agora, ao aspecto prático das etapas dos procedimentos de simulação. 


Etapas da simulação 


As etapas básicas de um processo de simulação correspondem às mesmas fases de um estudo 
de PO ou de uma pesquisa estatística científica. Compreendem: 


e formulação do problema; 

e determinação dos objetivos; 
e planejamento; 

e construção do modelo; 

e coleta dos dados; 

e codificação; 


* testes; 
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e validação; 
e produção; 
e avaliação; 
e documentação; 


e implementação. 


Cabe ressaltar a importância destas duas últimas (documentação e implementação), pois é fun- 


damental documentar de forma clara e concisa os resultados de todo o processo de simulação. 


Veja, na sequência, um exemplo de simulação apresentado por Santos (1999). 


Exemplo de modelo de simulação 


Uma grande máquina industrial tem 3 rolamentos diferentes que quebram de tempos em 


tempos. A probabilidade da vida útil (em horas de operação) de um rolamento está dada 


na tabela [...] [a seguir]: 


Vida do rolamento (horas) Probabilidade 
1.000 0.10 
1.100 0.13 
1.200 0.25 
1.300 0.13 
1.400 0.09 
1.500 0.12 
1.600 0.02 
1.700 0.06 
1.800 0.05 
1.900 0.05 


Quando um rolamento quebra, a máquina para e um mecânico é chamado para instalar 


um novo rolamento no lugar do que quebrou. 


O tempo que o mecânico demora para chegar ao rolamento quebrado também é uma va- 


riável aleatória, com a distribuição dada na tabela a seguir: 


Tempo de espera (minutos) Probabilidade 
5 0.60 
10 0.30 
15 0.10 
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Cada minuto que a máquina fica parada custa $ 5 e o custo mecánico é de $ 1 por minuto 
trabalhado substituindo o rolamento. O mecánico demora 20 minutos para trocar 1 rola- 
mento, 30 minutos para trocar 2 e 40 minutos para trocar 3. Cada rolamento novo custa 
$ 20. Alguém sugeriu que, ao quebrar um dos rolamentos, se fizesse logo a troca dos 3. 
Deseja-se avaliar a situacáo do ponto de vista económico [grifo nosso]. 

Solução 

Temos de comparar o custo de alternativa atual e da alternativa proposta. Precisamos es- 
tabelecer um horizonte de tempo para fazer esta comparação. Considerando que a menor 
vida útil de um rolamento é 1.000 horas (mais de 1 mês), vamos estabelecer um horizonte 
de 20.000 horas (um pouco mais de 2 anos) para fazer a comparação. 

Como a vida útil dos rolamentos e a espera pelo mecânico são variáveis 
aleatórias que seguem as distribuições vistas anteriormente, temos de relacionar àquelas dis- 
tribuições com uma tabela de números aleatórios. 

Assim sendo, vamos imaginar que temos um gerador de números aleatórios capaz de 
criar qualquer inteiro entre 0 e 99, ou seja, 100 números. Vamos atribuir a cada duração 
de vida útil uma faixa destes números que me garanta que a distribuição probabilística 
seja mantida. 

Como a 1º vida útil (1.000 horas) tem 10% de probabilidade de ocorrer, vamos atribuir a 
esta duração a faixa de 0 a 9 inclusive. Ou seja, 10% dos 100 números. Para a 2º duração 
provável (1.100 horas), com 13% de probabilidade de ocorrência, vamos atribuir a faixa 
de 10 a 22 inclusive, ou seja, 13 números. Podemos continuar para as demais durações 
prováveis dos rolamentos como pode ser visto na tabela a seguir, ressaltando que a proba- 


bilidade acumulada dá o limite das faixas escolhidas. 


ned Pribina — » Eoü—— 

m.) Acumulada Atribuído 
1.000 0.10 0.10 0-9 
1.100 0.13 023 10-22 
1.200 025 0.48 23-47 
1.300 013 0.61 48-60 
1.400 0.09 0.70 61-69 
1.500 0.12 0.82 70-81 
1.600 0.02 0.84 82-83 
1.700 0.06 0.90 84-89 
1.800 0.05 0.95 90-94 
1.900 0.05 1.00 95-99 


Tabela semelhante pode ser construída para a espera pela chegada do mecânico. 
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Tempo án m AE 
cure [omes Nao | esa. 
5 0.60 0.60 00 — 59 
10 0.30 0.90 60 — 89 
its 0.10 1.00 90 — 99 


Com os dados das tabelas [...], podemos executar a simulação que, neste caso, foi realizada 


na planilha [do Microsoft? Office] EXCEL, apresentando os seguintes resultados para o 


rolamento 1: 
ROLAMENTO 1 
Sequência D o Vida Id D . Espe E 
Aleatório (horas) Aleatório (Min) 
(horas) 
1 62 1.400 1.400 61 10 
2 85 1.700 3.100 10 5 
3 89 1.700 4.800 46 5 
4 24 1.200 6.000 28 5 
5 99 1.900 7.900 55 5 
6 27 1.200 9.100 64 10 
y 89 1.700 10.800 63 10 
8 2 1.100 11.900 75 10 
9 2 1.000 12.900 54 5 
10 34 1.200 14.100 67 10 
11 7 1.000 15.100 90 115) 
12 75 1.500 16.600 14 5 
i3 22 1.100 17.700 80 10 
14 97 1.900 19.600 84 10 
15 97 1.200 20.800 9 5 
$2120 


Podemos observar na planilha que para cada sequéncia, ou seja, rolamento novo, é gera- 


do um número aleatório que indica qual a vida útil daquele rolamento. Tendo quebrado, 


após esta vida útil, o mecânico é chamado e um 2º número aleatório é gerado para definir 


o tempo de espera até a troca do rolamento ser iniciada. 


Quando a vida acumulada ultrapassa 20.000 horas, ou seja, a duração da simulação, para- 


mos a execução do processo. 


Processos semelhantes foram executados para os outros 2 rolamentos, como visto a seguir. 
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ROLAMENTO 2 


Sequéncia N : M A ont D, ; Espe ix 
Aleatório (horas) Aleatório (Min) 
(horas) 
1 89 1.700 1.700 58 5 
2 47 1.200 2.900 88 10 
3 60 1.300 4.200 20 5 
4 3 1.000 5.200 98 15 
5 40 1.200 6.400 26 5 
6 64 1.400 7.800 97 15 
7 9 1.000 8.800 41 5 
8 30 1.200 10.000 79 10 
9 32 1.200 11.200 0 5 E 
10 8 1.000 12.200 3 5 $ 
M 94 1.800 14.000 58 5 E 
12 66 1.400 15.400 84 10 š 
13 53 1.300 16.700 61 10 8 
14 17 1.100 17.800 43 5 E 
15 72 1.500 19.300 15 5 E 
16 0 1.000 20.300 97 15 5 
2-130 E 
ROLAMENTO 3 : 
Sequéncia ps : M ge cm I . Espera 
Aleatório (horas) horas) Aleatório (Min) 
1 49 1.300 1.300 44 5 
2 26 1.200 2.500 45 5 
3 2 1.000 3.500 72 10 
4 83 1.600 5.100 87 10 
5 21 1.100 6.200 19 5 
6 20 1.100 7.300 81 10 
7 60 1.300 8.600 56 5 : 
8 34 1.200 9.800 74 10 E 
9 63 1.400 11.200 93 15 5 
10 69 1.400 11.200 93 15 i 
11 44 1.200 13.800 71 10 i 
12 76 1.500 15.300 97 15 : 
13 55 1.300 16.600 59 5 i 
14 85 1.700 18.300 81 10 i 
15 21 1.100 19.400 21 5 E 
16 5 1.000 20.400 1 5 E 
L=130 


Com os dados obtidos na simulação, podemos calcular o custo da situação atual: 


E 
2 


Custo dos rolamentos = (15 + 16 + 16) - $ 20 = $ 940 
Custo da máquina parada esperando pelo mecânico = 
(120 + 130 + 130) - $ 5 = $ 1.900 

Custo da máquina parada trocando rolamento = 

(15 + 16 + 16) - 20 - $ 5 = $ 4.700 

Custo do mecânico = (15 + 16 + 16) - 20 - $ 1 = $ 940 
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Custo Total = 940 + 1.900 + 4.700 + 940 = $ 8.480 


A simulação da situação proposta apresentou os seguintes resultados: 


ROL. 1 ROL. 2 ROL. 3 
i E : Vida ESPERA 
Seq.| NA Vida Hr | NA |Vida Hr | NA | Vida Hr 

Quebra | Acum. NA | (Min) 
1 96 | 1.900 2 1.000 34 | 1.200 1.000 1.000 | 21 5 
2 pá) | iso) 7 1.000 | 47 | 1.200 1.200 2.000 | 36 5 
3 | 96 | 1.900 | 46 | 1.200 | 49 | 1.300 1.200 eom | 2 5 
4 | 48 | 1.300 17 | 1.100 | 42 | 1.200 1.100 4.300 7 5 
S sel 1200 | 99] dc | 20 353099 1.100 5.400 | 58 5 
6 | 36 | 1.200 | 94 | 1.800 | 98 | 1.900 1.200 6.600 | 83 10 
7 Jj 4 p uem i | uen | ee p 9X0 1.100 7.700 | 14 5 
e pi | 1590 2 1.000 | 20 | 1.100 1.000 8.700 | 75 10 
9 4 1.000 | 22 | 1.100 86 | 1.700 1.000 9.700 5 5 
10 | 69 | 1.400 21 | 1.100 0 1.000 1.000 | 10.700 | 65 10 
Tb ome fp ote 89 | 1.700 | 58 | 1.300 1.100 | 11.800 | 15 5 
120567 [12007 | 12 | 100 66 | 1.400 1.100 | 12.900 | 12 5 
ls 75S SN 5728 Boo E 1200 1.200 | 14.100 | 32 5 
14 | 76 | 1.500 | 78 | 1.500 | 95 | 1.900 1.500 | 15.600 | 2 5 
15 | 71 | 1.500 5 1.000 86 | 1.700 1.000 | 16.600 | 31 5 
16 | 98 | 1.900 | 43 | 1200 | 22 | 1.100 1.100 | 17.700 | 51 5 
17 | 98 | 1.900 | 47 | 1.200 | 60 | 1.300 1.200 | 18.900 | 20 5 
18 | 68 | 1.400 61 | 1.400 5722 5915900 1.800 | 20.200 | 35 5 
105 


NA = Nº Aleatório 


Feita a simulação da situação proposta, podemos calcular os custos: 
Custo dos rolamentos = (18 - 3) - $ 20 = $ 1.080 

Custo da máquina parada esperando pelo mecânico = 105 -$5=$525 
Custo da máquina parada trocando rolamento = 18 - 40 - $ 5 = $ 3.600 
Custo do mecânico = 18 - 40 - $ 1 = $ 720 

Custo Total = 1.080 + 525 + 3.600 + 720 = $ 5.925 


Assim a simulação nos mostrou que a situação proposta é bem melhor em termos econô- 


micos. 


Fonte: Santos, 1999, p. 6-10. 


A simulação é uma técnica para estudar o comportamento e as reações de situações reais em 
modelos matemáticos, a qual permite dimensionar suas possíveis soluções. Ela possibilita uma 
análise probabilística com base em dados razoavelmente confiáveis. Nesse encadeamento de 
processos, trabalhamos com os coeficientes da função objetivo e as constantes de restrições, 
seguindo os passos necessários para manter a confiabilidade do processo. Portanto, sempre 


devemos levar em conta aspectos como: formulação do problema, determinação dos objetivos, 
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planejamento, construção do modelo, coleta dos dados, codificação, testes, validação, produ- 
ção, avaliação, documentação e implementação. Assim, estaremos aptos a fazer uma análise de 
sensibilidade, considerando, como vimos, a função desta, que é tornar as decisões mais preci- 
sas, facilitar o processo de refinamento dos dados estimados, os quais irão servir de base a tais 
decisões, bem como centrar o foco nos pontos essenciais ou críticos quando da implementação. 


1. Considere o exercício resolvido referente à empresa de artefatos de madeira apresentado 
no Capítulo 2. Faça a análise de sensibilidade para os coeficientes da função objetivo (de ca- 
da uma das variáveis originais, exceto da variável de folga) e para a constante da restrição. 
2. Considere o problema de PL: 


max Z = 4x, + 2x,+ X, 


s.a. x, + 2x,* 0,5x, € 10 
2x, *X,* X,X 8 


x,20,x,20 


A tabela final é dada a seguir: 


x X X x X 


1 2 3 4 5 
z=16 0 0 1 0 2 

Xa 6 0 1,5) 0 jl -0,5 

zen 4 1 0,5 0,5 0 0,5 


Determine o que irá acontecer com a solução ótima se o coeficiente de x, na função objetivo 
passar de 4 para 3. Caso haja alteração na base, determine a nova solução ótima. 


3. Faça a análise de sensibilidade para os coeficientes das variáveis originais do problema (x, 
e x,) e para a constante da restrição referente ao tempo disponível de molho (segunda restri- 
ção) do problema 3 do Capítulo 3: 


Um fabricante produz dois tipos de aço: normal e especial. Uma tonelada de aço normal re- 
quer 2 horas no forno de soleira aberta e 5 horas de molho; 1 tonelada de aço especial requer 
2 horas no forno de soleira aberta e 3 horas de molho. O forno de soleira aberta está disponí- 
vel 8 horas por dia, e o molho está disponível 15 horas por dia. O lucro em 1 tonelada de aço 
normal é de $ 120,00 e de $ 100,00 para 1 tonelada de aço especial. A empresa precisa produ- 
zir diariamente no mínimo 2 toneladas de aço normal e 1 tonelada de aço especial. Considere 


que a tabela final do problema é dada a seguir: 
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X 


1 2 3 5 6 7 8 

z = 406,66 0 0 0 333 | 467 |-4666*M 0 M 
x. 0,66 0 0 1 |-0,66 | -133 1,33 0 0 
X 0,66 0 0 0 0,33 | 1,66 -1,66 1 -1 
X; 2,0 1 0 0 0 -1 1 0 0 
X; 1,66 0 1 0 0,33 | 1,66 -1,66 0 0 


4. Qual é a importância da análise de sensibilidade em problemas de PL? 


5. O que é simulação e quais são as suas principais aplicações? 
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Conteüdos do capítulo 
e Passo a passo do desenvolvimento do conceito de matriz. 
e Formas de representação de dados correspondente às matrizes. 


e Possibilidades de cálculos aritméticos com as matrizes. 


Após o estudo deste capítulo, você será capaz de: 
1. identificar e relembrar conceitualmente as matrizes; 
2. classificar uma matriz quanto ao seu tipo e formato; 


3. recordar o cálculo das operações fundamentais de matrizes em seus diversos tipos. 


A pesquisa operacional (PO) trabalha muito com a linguagem matemática, com a programa- 
ção linear (PL), cuja base está centrada no estudo das matrizes. Por isso, neste capítulo aborda- 
remos alguns tópicos relacionados a matrizes. Esse é um conteúdo que você precisa conhe- 
cer para utilizar a PO em tomadas de decisão, sendo esta a temática ou o foco que esta obra 
privilegia — o uso da PO em processos decisórios. 


8.1 Conceituando matriz 


Matriz é um conjunto retangular de números que pode ser representado entre colchetes ou 
entre parênteses. Você já deve ter se deparado com essa definição, ou similares, durante seus 
estudos de ensino básico. Mas como isso funciona? A que se aplica? 


Para entender melhor esse tema, vamos considerar o exemplo da situação de um fabricante 
de automóveis que produz três modelos de carros nas versões X, LX e GLX. 


Vamos, então, supor que cada modelo é manufaturado parcialmente na fábrica F,, localizada 
em Londrina - PR, e, depois, terminado na fábrica F,, em São José dos Pinhais — PR. 


Nessa situação, para realizarmos os cálculos relacionados aos custos de produção dos automóveis, va- 
mos lembrar que o custo total de cada produto é a soma do custo de produção com o custo de transporte. 
Os custos de cada fábrica são apresentados a seguir. 


Custos de produção e transporte referentes à fábrica F, 


Custo de produção Custo de transpor- 
(R$) te (R$) 
X 12.000,00 100,00 
LX 11.000,00 90,00 
GLX 14.400,00 95,00 


ão dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 
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Custos de produção e transporte referentes à fábrica F, 


Custo de producáo Custo de transpor- 
(R$) te (R$) 
X 9.500,00 97,00 
LX 14.300,00 102,00 
GLX 15.700,00 96,00 


Essas tabelas podem conceituar matematicamente o que chamamos de 


matrizes e podem ser escritas como no esquema apresentado a seguir, entre parénteses: 


12.000 100 9.500 97 
E|11.000 90 |eE, 14.3800 102 
14.400 95 15.700 96 


Isso porque toda a disposição numérica no formato da primeira tabela pode ser distribuída 


entre linhas e colunas, precedidas por parênteses ou colchetes (como fizemos agora, em F, e F,). 


Genericamente, essa representação (entre parênteses) pode ser entendida como uma matriz 
qualquer por meio de um arranjo retangular composto por mXn números reais! (ou comple- 
xos), organizados em m (linhas horizontais) e n (colunas verticais), que representam a ordem 
mXn (lemos "m por n" ou “mn”. 


É importante ressaltar que: 


e usualmente sempre se dá um nome para a matriz, por meio de uma letra maiúscula 
latina. Por exemplo, A; 


e os subscritos (geralmente acompanham a mesma letra utilizada para a matriz, agora 
minúscula) significam respectivamente a linha e a coluna ocupadas pelos elementos. 
Por exemplo, a, ou aj 

e mindica a quantidade de linhas que uma matriz apresenta; 


* nindicaa quantidade de colunas que uma matriz apresenta. 


Em outras palavras, é a ordem mXn, que representa o tamanho da matriz e que entendemos 
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como a quantidade de linhas (m) e de colunas (n), como já ressaltamos. 


Como você pode verificar, os nomes das matrizes mostradas são F, e F,, isto é, foram nome- 
adas com uma letra maiúscula latina (F). Além disso, apresentam, respectivamente, a ordem 
3X2. Isso porque a ordem, entendida como o tamanho, é definida por três linhas (horizontal) 
e duas colunas (vertical). 


1 Números reais são os números do conjunto R, cujos elementos são números decimais, fracionários, negativos e inteiros. 
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Em virtude de nào podermos exemplificar todos os formatos de matrizes, generalizamos o 


conceito de matriz pela representação a seguir: 


Essa matriz, nomeada como A, apresenta ordem mXn (nào sabemos quantas colunas e linhas 


ela terá, mas podem ser muitas, por isso a generalização mXn). 
Matematicamente falando, em particular, a i-ésima linha de A é dada por: 
A -(a4a,..a,), (1€iX m) 


E a j-ésima coluna de A é escrita como: 


Vamos explicar essa representação. Cada número a, que está dentro da matriz é chamado 
de elemento da matriz, pois esses índices (i e j) indicam a posição do elemento dentro da ma- 
triz. Vejamos: 


a do ee Ain 


e Anotação a, indica que esse elemento da matriz está na primeira linha (i = 1) e na pri- 
meira coluna (j = 1). 


e Anotação a, indica que esse elemento da matriz está na segunda linha (i = 2) e na se- 


gunda coluna (j = 2). 


e A notação a,, que você nào encontra grafada na representação (imediatamente ante- 


327 
rior) da matriz, embora possa ser presumida pela disposição gráfica, indica que esse 
elemento da matriz está na terceira linha (i = 3) e na segunda coluna (j = 2) e assim su- 
cessivamente, como foi possível observar na representação gráfica que utilizamos para 


a demonstração da matriz com seus respectivos elementos. 
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* A leitura da composição que representa um elemento da matriz (expresso geralmen- 
te pela letra representativa da matriz em fonte minúscula mais os índices), como a,, 
é feita da seguinte forma: “a três dois". Exemplificando com outras composições: 


a, (a um um); a,, (a dois dois). 


Dessa maneira, os subscritos significam respectivamente a linha e a coluna ocupadas pelos 
elementos. Logo, vamos repetir, se a matriz for A (letra maiúscula), o elemento da matriz se- 
rá a (letra minúscula) com os numerais representativos da linha e da coluna (subscritos), co- 


mo em a,,. 
Como fazemos, então, a indicação de uma matriz? 


A notação para a matriz pode ser realizada da seguinte maneira: podemos, por exemplo, in- 
dicar uma matriz A pelos seus elementos, na forma A, -[ajl. Nela, mXn indica a quanti- 
dade de linhas e colunas, e a, indica quais sào os elementos da matriz, ou seja, quais sào os 
elementos a,, a, ..., à, sendo estes expressos por expressões matemáticas, como veremos 
nos exercícios resolvidos. 


Hyer 


| 
L/C d 


Sabendo que a matriz A, é definida por a, = 2i*j, construa a matriz A. 
Resolução 


Se a matriz A é de ordem 3X2, sabemos que ela tem 3 linhas e 2 colunas e cada elemento a, 
será descrito pela expressão matemática 2i+j, já que não sabemos quais são os elementos des- 


sa matriz. Assim, temos: 
a = 2*jf 


a=2:1+1=2+1=3>a,53; a,=2:1+2=2+2=4>a,„=4; 


a,=2:2+1=4+1=55a,=5; a,=2:2+1=4+2=65a,,=6; 


2 or 1l=6+1=7 »8477; a,=2:3+1=6+2=8 >a, = 8; 


Logo, a matriz A será: 


34 
A=|5 6 
7 8 
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Sào inümeros os tipos de matrizes, entre os quais podemos citar as matrizes linha, coluna, 
adjunta etc.; porém, chamaremos a atenção somente para quatro casos, por serem o objetivo 
e o foco do nosso interesse. Sào as matrizes classificadas como quadrada, diagonal, identi- 
dade e transposta. 


Matriz quadrada 


Como vocé pode imaginar, a matriz assim denominada apresenta-se sempre em um formato 


quadrado, ou seja, nesse caso, a quantidade de linhas é igual à quantidade de colunas. 


Conceituando matematicamente, dizemos que, se tivermos uma matriz Ae se m - n, a matriz 


é dita matriz quadrada de ordem n, e os elementos a,, a a, formam a diagonal princi- 


LH or vt 


pal da matriz A. 


Você pode observar essa formulação de matriz tanto na matriz B como na M que apresenta- 


mos a seguir: 


10 À 
2 a 
B=| 3 2 0,4|,M= 
412 
6 0 -4 


Você pode notar que a matriz B é de ordem 3X3, pois apresenta três linhas e três colunas; já 
a matriz M é de ordem 2X2. Não haveria nenhum problema em dizer que as matrizes B e 
M são de ordem 3 e 2, respectivamente (lembrando que m = n; logo, temos n =n e, portanto, 
ordem n). 


Podemos observar também que os elementos da diagonal principal das respectivas matrizes 
são os valores 1, 2 e -4 para B e -1 e 12 para M. 


Dessa maneira, se dissermos que temos uma matriz de ordem 10, imediatamente lembrare- 
mos que essa matriz é quadrada e apresenta 10 linhas e 10 colunas. 


Matriz diagonal 


Uma matriz é dita diagonal sempre que todos os elementos, acima ou abaixo, ou ambos, da 


diagonal principal sào iguais a zero, ou seja, os valores a 


ar Apr eee A (sendo i = j) são todos di- 


ferentes de zero, e os termos em que a, sendo i +j, apresentam valores iguais a zero. 


Observe as matrizes a seguir. Podemos notar algumas diferenças e semelhanças entre elas, como: 
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e todas apresentam uma diagonal (pois todo formato quadrado tem uma diagonal, e es- 


sa diagonal está na posição em que se encontram os termos a,, a,,, ..., a, (sendo i =j); 


e algumas têm valores abaixo dessa diagonal igual a zero, outras acima, podendo ser 


acima e abaixo, como na matriz T. 
Vejamos: 


Dadas as matrizes T, G e C: 


-1 2 0 1 
1 0 0 
20 0. -3 -3 2 
T=|0 2 0,G- ,Cz 
-11 0 0 13 
0 0 3 
0 0 0 4 


podemos interpretar que: 
e amatriz T é diagonal de ordem 3X3; 


e a matriz G é diagonal inferior de ordem 2X2, já que os valores acima da diagonal são 
iguais a zero; 


e  amatriz C é de ordem 4X4, e ela é chamada de matriz diagonal superior. 


Atenção: perceba que toda matriz diagonal deve ser quadrada. 


Matriz identidade 


É toda matriz diagonal que apresenta todos os elementos da diagonal principal iguais a 1 (po- 
sitivo) e os demais elementos iguais a zero. Veja os exemplos das matrizes T, Ge C a seguir: 


Sono 
oh o o 
= ooo 


1 
0 
0 
0 


Importantíssimo reforçar: Observe, também, que ela é quadrada, e os elementos da diago- 
nal principal são iguais a 1. 
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FK 


» ] e - 
V EXercicio 


e aplicando o que você aprendeu, qual das afirmações a seguir é verdadeira? 
a) A matriz I é matriz identidade de ordem 2. 

b) A matriz G é conhecida como matriz identidade de ordem 3, 

c) A matriz F é uma matriz diagonal. 

d) A matriz G é conhecida como matriz diagonal. 

Resolução 

Para simplificar, apenas a alternativa “b” está correta. 


A alternativa “a” é falsa, pois, para ser considerada identidade, a matriz deve ter em sua dia- 
gonal elementos iguais a 1, e nào -1. A alternativa "c" está incorreta porque apresenta um 
elemento diferente de zero fora da diagonal principal (a, = —1). A última alternativa também 


é falsa, pois a matriz G é, na verdade, a matriz identidade. 


Matriz transposta 


A representação que fazemos logo a seguir, já vista no início deste capítulo, apresenta as ma- 
trizes "custos" das tabelas “Produção e transporte F,” e “Produção e transporte F,” dispos- 
tas de outra forma se comparadas à configuração das tabelas. Essas tabelas você encontra na 


sequência. 
12.000 100 9.500 97 
F,| 11.000 90 |eF,| 14.300 102 
14.400 95 15.700 96 


Custos de produção e transporte referentes à fábrica F, dispostos em linhas 


X LX GLX 


Custos de produção 12.000,00 11.000,00 14.400,00 


Custos de transporte 100,00 90,00 95,00 
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Custos de produção referentes à fábrica F, dispostos em linhas 


X LX GLX 


Custos de produção 9.500,00 14.300,00 15.700,00 


Custos de transporte 97,00 102,00 96,00 


No entanto, em vez de os custos de produção e de transportes estarem escritos em colunas, 
como nas tabelas que acabamos de ver, podemos apresentá-los em linhas, conforme as repre- 


sentacóes que faremos na sequéncia, resultando nas seguintes matrizes: 


1 


2 


pe 12.000 411.000 14.400 E! 9.500 14300 15.700 
100 90 95 97 102 96 


A matriz F’ é chamada de transposta de F, e a matriz F’ é a transposta de F,. Podemos dizer, 


então, que, se A é uma matriz, então a sua transposta é a matriz AT, em que: 


a; -a, com 1<i<n,1<j<m 


Assim, a transposta de uma matriz A é obtida trocando-se as linhas de A por suas colunas e 


vice-versa. 


11 -5 11 
0 0 EB 14 1 14 

zb ;|6- i AU E 7 E mL x B E j 
= 0 0 1 1 3 4 14 14 14 
14 -14 14 

28 28 28 | 


vocé pode afirmar: 

a) A matriz F é diagonal. 

b) A matriz G é de ordem 4x4. 

C) A matriz H é de ordem 3x3. 

d) A matriz Ié a matriz de 4x4. 
Resolução 


Analisando uma a uma cada alternativa, podemos verificar que: 
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Ma! 


e naalternativa “a”, a matriz F não é diagonal, pois, para ser diagonal, os elementos aci- 


ma ou abaixo, ou ambos, deveriam ser iguais a zero; 


e na alternativa "b", a matriz G tem três linhas e três colunas; portanto, ela é de ordem 
3x3, ou simplesmente de ordem 3; 


e naalternativa “c”, a matriz H é de ordem 3, ou seja, 3x3. Então, essa é a alternativa correta; 


* na alternativa "d", a matriz I nào pode ser de ordem igual a 4x4, pois ela apresenta 
apenas três linhas e três colunas. 


Nosso objetivo, ao descrevermos brevemente como as matrizes se apresentam, é que você 
esteja atualizado para continuar os caminhos da PO, pois nessa área trabalhamos com dados 
quantitativos modelados com a linguagem matricial. Portanto, nesse contexto, a compreen- 
são dos conceitos de matrizes é fundamental. 


Cn 
Qo 


A 


com ma 
Com as matrizes podemos também efetuar algumas operações, como já é de seu conheci- 
mento por meio do estudo da aritmética. Lembre-se de que as operações de soma, subtração, 
multiplicação, divisão e potência foram algumas das diversas operações aritméticas que você 


aprendeu durante sua jornada no ensino básico. 


Demonstraremos alguns cálculos matriciais, apenas aqueles necessários para que você se si- 
tue no contexto da PO: igualdade de matrizes, multiplicação por um escalar, soma de matri- 


zes e multiplicação de matrizes. 


Igualdade de matrizes 


Dizemos que duas matrizes A = aje Bs b, são iguais se a, = bi, isto é, se todos os elemen- 


tos correspondentes forem iguais. Vejamos o exemplo a seguir. 


Sendo as matrizes A e G de ordem 3: 
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123 1 b 3 
A=|0 1 1|eG=|0 c 1 : 
0 2 1 a 2d : 


A =G, se e somente se a= 0, b=2,c=1ed=1. 


Observe a sentença, pois ela diz que A = G, se e somente se a =0, b =2, c=1 e d =1. Visualizando 
essa troca na matriz G, ficará claro o processo de igualdade. 
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Multiplicação por um escalar 


Nessa operação aritmética, a matriz B dada pelo produto de uma matriz A por um escalar k 


(um nümero qualquer) é obtida multiplicando-se cada elemento de A por k. Vejamos: 


BA E up dusdekuic ^ OO d 
e À= 15 ek-5,logo, o produto de 754 55 , gerando, 
10 20 
assim, a matriz B= : 
5 25 


Na linguagem matemática, podemos definir a multiplicação de um escalar pela expressão: 
B=k-A<ob,=ka,emquel<ism, 1<j<n 
Vejamos outro exemplo: 


A empresa Trabalhando Muito Ltda., em 2002, tinha os custos de produção e de transporte 
dos produtos P1 e P2 fornecidos pelas tabelas que apresentamos logo a seguir, cotados em 
dólar. Se quisermos transpor a tabela para valores em reais, da época, basta multiplicarmos 
cada elemento da tabela por 1,994, como apresentado a seguir: 


Custos dos produtos P1 e P2 


Custos P1 P2 


Custos de produção 2,33 | 15 


Custos de transportes 134 | 20 


Se transformarmos em matriz a tabela mostrada, teremos: 


2,33 LB 
134 2,0 


, sendo C a matriz custo. 


Para fazermos essa conversão de dólar para real, basta multiplicarmos o escalar 1,994 pela 
matriz C, ficando: 


[2,331,994 1,5.1,994| [4,65 2,99 
1,34-1,994 2,0-1,994 


i 2,67 3,99 


Esse resultado apresenta os custos de transporte e de produção, em reais, dos produtos P1 e P2. 


Soma de matrizes 


Se quisermos efetuar a soma de matrizes, basta somarmos os elementos correspondentes de 
cada matriz. É importante ressaltar que as matrizes devem apresentar o mesmo nümero de 
linhas e o mesmo nümero de colunas. Em outras palavras, se A e B sào matrizes, a soma de 
A e B resulta na matriz C definida por: 


c=a+b. 
1 1j E) 


com 1<i<m,1<j<n 
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Por exemplo, os custos totais de produção e de transporte de cada um dos modelos de auto- 
móveis, se utilizarmos os dados do início deste capítulo, sào fornecidos pela soma dos custos 
de produção e de transporte dos automóveis nas fábricas F, e F,. Sendo F, a matriz custo de 


produção e F, a matriz custo de transporte, temos, portanto: 


12.000 100 9.500 97 
F| 11.000 90 |eF,| 14.300 102 
14.400 95 15.700 96 


Os custos totais de produção são definidos pela soma dos custos de produção e de transpor- 
te. Logo: 


12.000 100 9.500 97 
C=F +F, > C=| 11.000 90 |+| 14.300 102 
14.400 95 15.700 96 


Se somarmos cada elemento da matriz, teremos: 


12.000 + 9.500 100+97 
C=| 11.000+14.300 904102 
14.400 + 15.700 95+96 


Portanto, a matriz custo total, chamada de C, será: 


21.500 197 
C=| 25.300 192 
30.100 191 


Assim, os custos totais de produção e de transporte dos três modelos de automóveis são es- 
pecificados na tabela a seguir. 


Custos totais de produção e transporte dos modelos de veículos X, LXe GLX 


Custos de produção | Custos de transporte 
X 21.500,00 197,00 
LX 25.300,00 192,00 
GLX 30.100,00 191,00 


Multiplicação de matrizes 


Para multiplicarmos duas matrizes quaisquer, devemos prestar muita atenção, pois essa opera- 
ção não se faz por meio do mesmo processo pelo qual estamos acostumados a multiplicar um 


número escalar pela matriz. Vamos exemplificar essa situação partindo do pressuposto de que, 


Nenhuma parte desta publicação poderá ser reproduzida por qualquer meio ou forma sem a prévia autorização da Editora InterSaberes. A violação dos direitos autorais é crime estabelecido na Lei nº 9.610/1998 e punido pelo art. 184 do Código Penal 


dada uma matriz A do tipo mXn e sendo B uma matriz do tipo nXo, define-se como produto 
da matriz A pela matriz B a matriz C do tipo mXo, tal que cada elemento de C (c,) satisfaça: 


c; = ab, + apb, + (..) + a,b 


inn 


Em outras palavras, cada elemento da matriz C é calculado pela multiplicação ordenada dos 
elementos da linha i da matriz A pelos elementos correspondentes da coluna j da matriz B e, 


a seguir, somamos os produtos obtidos. 


[[2.3+3-2] 241494 
AsBse| 13400 1-1+0-4 
| 43492 4.14 5.4 


12 14 


l| AB. AB 
| 22 24 RANA 


Observação importantíssima: Só podemos multiplicar matrizes quando o número de colu- 


nas da primeira matriz for igual ao número de linhas da segunda matriz. Isso você pode ob- 


servar na representação que apresentamos anteriormente. 


Agora vamos realizar alguns exercícios para que você possa acompanhar a aplicação desses 
cálculos. 


Se você acompanhar a resolução desses exercícios, com certeza poderá entender de forma 
mais rápida o que abordamos neste capítulo. 
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1. Na confecção de três tamanhos de camisetas de futebol (P, Me G) são usados distintivos 
dos clubes nas cores verde (V) e vermelha (v). O número de distintivos usados, por tamanho, 
é dado pela tabela a seguir: 


Camiseta P Camiseta M Camiseta G 
v 3 1 3 
V 6 5 5 
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O nümero de camisetas fabricadas, de cada tamanho, nos meses de maio e junho, é dado pe- 
la seguinte tabela: 


Maio junho 
Camiseta P 100 50 
Camiseta M 50 100 
Camiseta G 50 50 


Nessas condições, obtenha, por meio do uso de matrizes, a tabela que dá o total de distinti- 
vos usados em maio e junho. 


Resolução 


A solução do problema se resume na multiplicação das matrizes: 


100 50 

313 500 400 
[50 100|= 

655 1.100 1.050 
50 50 


Logo, a tabela obtida foi: 


Maio junho 
v 500 400 
V 1.100 1.050 
2. Dadas as matrizes B, D e C: 
-1 3 2 » Og. 
» Ji 
B=|3 -2 “dA au w, 
6 3 5 -2 2 3 


calcule as expressóes a seguir: 
a)3-C 

b)B-D 

c)D+B 


d)C-D 
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Resolução 


a) Multiplicamos cada elemento da matriz C por 3, pois se trata de multiplicação de um es- 


calar por uma matriz. 


Š 2 1| |32 31) |6 3 
1 4X 131 Ga] 13 3 
b) Temos um caso de subtração de matrizes. Logo, devemos subtrair cada elemento da ma- 


triz B pelo respectivo elemento da matriz D. 
=]. 3 2 -2 0 1| |-1-(-2) 1 3 1 


3-0 2-1 
B-D=|3 -2 0-|1 0 2| 3-(-1 -2-0 0-2 |=|4 -2 -2 
6 3 5 2283 6-(-2) 3-2 5-3 8 1 2 


Observação: Lembre-se da regra de sinais nas operações com números negativos. 


c) Somando as matrizes diretamente, temos: 


-2 0 1 -13 2 -2+(-1) 0+3 2-1 -3 3 3 
D+B=|-1 0 24|3 -2 0|=|-1+3 0+(-2) 2+0|=]2 -2 2 
-2 2 3| |6 Emo -246 243 3+5 4 5 8 


d) Como a matriz C apresenta duas colunas e a matriz D apresenta três linhas (2 7 3), o pro- 


duto C - D é impossível. 


e) Para a multiplicação da matriz B pela D, temos: 
SE PM D V 


B-D=| 3-2 0||-1 0 2 
6 3E? 9 


5 


(79: C2)* (3) (D+(2) (2) (799: (0) (3) (0) +(2)-(2) C 0) (3) (2)+(2)-(3) 
(3) (2) + (22) (C + (0): (—2) (3) (0)*(72) (0) + (0)-(2) (3) (1) +(-2) - (2) + (0) (3) 
(6): (-2)* (3) CD (5: C2). (6) - (0)+ (3) (0) € (5) (2). (6) -+ (3) -(2)+ (5): () 


2-3-4 040-4 —1+6+6 
-6+2+0 00-40-40 3-4+0|= 
|-12-3-10 0+0+10 +6+6+15 


-5 4 1 
-4 0 -1 
-25 10 27 
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| Sintese 


O contexto matemático da PO cria a necessidade de domínio do conceito de matrizes, da ela- 
boração estrutural destas com conhecimento de seus elementos básicos. Esses conhecimentos 
são prerrequisitos. Por esse motivo, agregamos ao nosso estudo de PO essa abordagem. Nos 
capítulos anteriores, você, com certeza, já deve ter percebido a necessidade do conhecimento 
das operações matriciais para a operacionalização da PO. Queremos, no entanto, deixar bas- 
tante claro que o conteúdo apresentado aqui constitui apenas um material de apoio para re- 


lembrar conceitos e operações que você deve ter aprendido na educação básica. 


At Questões para revisão 


Estes exercícios são fundamentais para que você perceba como o seu aprendizado está se de- 
senvolvendo. E muito importante que você faça todos, de forma a reforçar a base teórica pa- 


ra os demais capítulos. 


1. Construa a matriz A = (a), de modo que a, -3i^ - j. 
1 2 -3 D d 
2. Dadas as matrizes A = qe eame ie B- 0 |, determine A +2B7. 
4 -3 
2 PT — 
3. Determine x e y sabendo que TAE 2 l , sendo x > 0. 
4 0 2x -y 0 


4. Defina matriz identidade. 


5. Quando podemos chamar uma matriz de quadrada? 
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Ao chegarmos ao final desta obra, queremos destacar o quanto a pesquisa operacional (PO) 
é importante quando tratamos de problemas que envolvem otimização de processos, seja para 


maximizar o lucro de uma empresa, seja para minimizar os custos de transporte. 


No entanto, vimos que a PO não se limita apenas a esses dois tipos de problemas. Temos uma 
grande variedade de processos que se enquadram nos temas por nós abordados. Podemos 
otimizar a utilização dos recursos de modo a diminuir o desperdício de matéria-prima e me- 
lhorar a programação e o controle da produção. Vimos também que a PO está baseada na 
construção de modelos matemáticos determinísticos (problemas estáticos) e probabilísticos 
(problemas dinâmicos). Isso explica o fato de, atualmente, muitas corporações recorrerem a 
profissionais dotados de conhecimento em PO para otimizar seus processos produtivos ou 
organizacionais. Cada vez mais esses profissionais estão presentes, executando suas ativida- 
des e, na maioria dos casos, contribuindo para que as empresas obtenham melhorias signifi- 
cativas nas atividades desenvolvidas. 


Muitos problemas apresentam uma grande quantidade de variáveis e restrições. Nesse caso, 
é importante a utilização de um software específico para a resolução de problemas. Destacam- 
-se o Solver (Microsoft? Office Excel), o Linear Interactive and Discret Optmizer (Lindo), 
o Language for Interactive General Optimizer (Lingo), o Quantitative System for Business 
(WinQSB), entre outros. 


É importante ainda ressaltar que no Brasil existe a Sociedade Brasileira de Pesquisa Operacional 
(Sobrapo!), que reúne profissionais de PO. A Sobrapo organiza simpósios anuais e publica sua 
própria revista científica, denominada Pesquisa Operacional, que é indexada nos International 
Abstracts in Operations Research, da International Federation of Operational Research Societies 
(Ifors) e, desde 2002, no SciELO. 


Este livro finaliza aqui, mas você pode ampliar seus conhecimentos na área dando continui- 
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dade aos seus estudos nas várias vertentes apresentadas nesta obra de iniciação à PO. Isso é 


possível por meio das leituras recomendadas na seção “Para saber mais” e de inúmeras ou- 
tras que você encontra na bibliografia sobre o assunto. Bons estudos! 
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V Respostas 


m d 


Capítulo 1 

1. V, V, V, B, V, P, V, V. 
2. d 

3. a 

4. Resposta pessoal. 


5. Os problemas são situações que a orga- 
nização precisa resolver para dar segui- 
mento aos seus propósitos ou para atin- 
gir seus objetivos. 


Capítulo 2 


1. Vocé deve ter encontrado o se- 
guinte resultado: x, - 10, x, - 0, 
z=30. 


Solução 
ótima 


2. Você deve ter encontrado a se- 
guinte solução: x, = 4,5, x, = 3,5, 
Z=28,5. 


24 


a 


2 


Solucào 
ótima 


"di 


> 


2x,<9 Aa 
e xtxi9 


3. Você deve ter encontrado a seguinte 
solução: x = 0, x,= 384, z = 11.520. 


Em que: x, é a quantidade de implemen- 
tos do modelo A e x, é a quantidade de 
implementos do modelo B. 


XA 
Ax, * 5x,* 1920 Solução 
otima 


3x,+ 3,5 x,º 1440 ⁄ 


x,+1,5x,º 576 


4. A função objetivo em PO direciona a 


construção, reconstrução, solução e va- 
lidação do modelo, sempre consideran- 
do o lucro máximo ou o custo mínimo. 
Nesse processo, atribuir valores aleató- 
rios para z é uma técnica que possibilita a 
representação gráfica da função objetivo. 


5. A forma-padrão de um problema de PL 


é caracterizada pela padronização com o 
objetivo de facilitar o entendimento. 
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Capítulo 3 


1. 


Vocé deve ter encontrado o seguinte re- 
sultado: x, - 65.040, x, = 100, z = 1.830.320. 


Em que: x, é a quantidade de carteiras e 
x, é a quantidade de bolsas. 


. Vocé deve ter encontrado o seguinte re- 


sultado: x, = 5,33, x, = 19,33, z = 2186,67. 


Em que: x, é a quantidade de biombos 
do modelo Ie x, é a quantidade de biom- 
bos do modelo II. 


. Vocé deve ter encontrado o se- 


guinte resultado: x, - 2, x, - 1,67, 
z = 406,67. 


Em que: x, é a quantidade, em toneladas, 
de aço normal e x, é a quantidade, em 
toneladas, de aço especial. 


Foram aplicadas operações de 


maximização. 


. Esse método baseia-se no acréscimo de 


uma variável artificial, em um processo 
semelhante à “folga”, visto no método 
simplex. A diferença é que nesse cálculo 
forçarmos essa variável a anular-se no 
fim. O “M” representa um número ar- 
bitrariamente grande, que podemos re- 
presentar por 10.000 ou mais. Por isso 
Big-M, ou seja, Grande M, é o coeficiente 
de penalização atribuído às variáveis. 


Capítulo 4 


1, 


2. 


3. 


O resultado é: x, = 30, x, = 100, z = 86.500. 


O resultado é: x, = 110.000, x, = 10.000, z 
= 81.500. 


O resultado é: x, = 6, x, = 2, z = 360. 


4. O resultado é: x, = 6, x, = 3, z = 177.000. 


5. O resultado é: x, = 2.000.000, x, = 4.722.222, 
x= 277.777,8, z = 15.726.110. 


Capítulo 5 


1. O resultado é: x, = 25, x, = 15, x 
X4, = 20, x4, = 30, z = 535. 


2. O resultado é: x,, = 50, x,, = 30, x,, = 40, 
X, = 40, X» 30, X, 20, z = 820. 


31 


3. O resultado é: x, = 4.000, x, = 3.000, x,, 
= 5.000, x, = 8.000, x = 4.000, 
z = 1.272.000 


1 DESTINO FICTÍCIO 


4. Os fatores que permitem identificar se o 
sistema está em equilíbrio ou desequilíbrio 


são a origem, o destino e a quantidade. 


5. Usar modelos matemáticos para a reso- 
lução de problemas de logística possi- 
bilita encontrarmos uma solução ótima 
para problemas que envolvem quan- 
tias a serem transportadas para outro(s) 


destino(s), além de otimizar as rotas. 


6. Os nós a serem conectados são: A-B, B-C, 
C-D, D-E e D-F. O custo total é: 167. 


7. Os nós a serem conectados são: A-B, A-C, 
F-D, C-E, E-F, F-G e G-H. O custo total 
é: 1.322. 
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Capítulo 6 


1. 


De Curitiba para Paranaguá: 16 mil. 


De Curitiba para Itajaí: 8 mil. 


Estoque nào utilizado de Curitiba: 1 mil. 


De São Paulo para Santos: 12 mil. 


Estoque nào utilizado de Sào Paulo: 8 


mil. 


Custo total de transporte: R$ 5.960.000,00. 


De Colombo para Curitiba: 700. 


De Colombo para Almirante Tamandaré: 
300. 


De Colombo para Pinhais: 300. 
De Campo Largo para Pinhais: 10. 


De Campo Largo para Campo Largo: 
540. 


De Campo Largo para Araucária: 310. 
Produtos nào recebidos em Pinhais: 90. 


Custo total de transporte: R$ 21.330,00. 


De A para R: 300. 
De B para P: 250. 
De B para Q: 100. 
De B para R: 100. 
Estoque nào utilizado de B: 50. 


Custo total de transporte: R$ 17.400,00. 


. Osnósaserem conectados são: D-B, A-C, 


C-D, C-E, D-F, D-G, I-H e G-I. O custo 
total é: R$ 107,00. 


. Osnós a serem conectados são: A-B, B-C, 


A-D, D-E, E-F, D-G, G-H, H-I, I-J, J-K, 
K-L, I-M, M-N e M-O. O custo total é: 
R$ 210,00. 


Obs.: A solução obtida não é única. Há 
outras possibilidades que resultam em 
custo total igual a R$ 210,00. 


Capítulo 7 


1. Coeficiente de x, menor ou igual a 5,00, 


a solução original é ótima (- o « c, € 5,0). 


Coeficiente de x, na funcáo objetivo maior 
do que 2,80, a solução original é ótima 
(2,80 < c, < 00). 


Constante da restrição maior ou igual a 
zero, a solução original é ótima (0 « b, < 


oo). 


. A solução ótima do problema original é: x, 


-4, x, - 0, z - 16. Pela análise de sensibili- 
dade referente à variável x, temos que 4 < 
c, «ov. Com c, passando de 4 para 3, temos 
que a solução original não é mais ótima. 
Resolvendo o novo problema pelo método 
simplex, temos que a nova solução ótima é: 
x=2,x=4,2=14. 


. Solução ótima para 


-œ< c, € 166,66, 72 « c, « o e 13 <b, < 16. 


. A análise de sensibilidade verifica quais 


são os efeitos causados pela variação de 
um determinado elemento na solução fi- 
nal de um problema de PL. Utilizamos a 
análise de sensibilidade para verificar o 
impacto de uma alteração no problema 
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original sem termos de resolver nova- 


mente esse problema. 


5. A simulação consiste no uso de técnicas 
matemáticas implementadas computa- 
cionalmente para imitar um determina- 
do processo do mundo real. Utilizamos 
a simulação para descrever o compor- 
tamento de um determinado processo, 
prever comportamentos futuros ou cons- 
truir teorias com base em determinadas 


observações. 
Capítulo 8 


1. Você deve ter chegado à seguinte matriz: 


2 10 
11 10 9 


2. Vocé deve ter chegado à seguinte matriz: 


3 8 5 
0 5-6 


3. O resultado para essa questão é: x = 3e 
y=2. 


4. É toda matriz diagonal que apresenta 
todos os elementos da diagonal princi- 
paliguais a 1 (positivo) e os demais ele- 


mentos iguais a zero. 


5. Conceituando matematicamente, dize- 
mos que, se tivermos uma matriz A e se 


m - n, então a matriz é dita quadrada. 
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Na prática, a PO fornece às empresas, aos estudantes e 
aos profissionais em geral uma série de métodos para 
diminuir tar lucros e, principalmente, 
encontrar a > para uma situacáo que 


envolva cálculos dedsvosn na vida de uma organização. 


Ideal para quem atua nas áreas de administração 
de empresas e engenharia, este livro é uma ferra- 
menta para quem deseja aprender, de maneira 
dinâmica e direta, a ler, entender e formular 
problemas de forma sistemática e racional, para, 
dessa forma, aplicar os conhecimentos adquiri- 
dos em sua prática profissional. 
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